Ocean Modeling Discussion

In the script designed to configure the nudging coefficients, the loop iterates from 1 to NTCLM. However, the array Tnudg contains NT elements. I think there is a mismatch of the indices IF (ANY(LnudgeTCLM(:,ng))) THEN DO itrc=1,NTCLM(ng) DO k=1,N(ng) DO j=JstrT,JendT DO i=IstrT,IendT CLIMA(ng)%Tnud...

Hi Tao, Thank you very much. I will read your paper and contact you. Wentao Ma Hi, As for the isopycnal flux, you could divide one section into a few density classes and calculate the flux across each class. As for the diapycnal flux, there are two ways: one is through the control volume method, and...

Hi folks
Do you have any idea to calculate isopycnal and diapycnal fluxes from the ROMS output? I know how to calculate the lateral and vertical flexes. But I need to decompose the fluxes to isopycnal and diapycnal.

Hi folks
Is there a way to calculate diapycnal velocity from the ROMS output?

I checked the forcing data in the dir. of the case test. The attribute should be as variables: double sst_time(sst_time) ; sst_time:long_name = "time for sea surface temperature" ; sst_time:units = "days since 1968-05-23 00:00:00 GMT" ; sst_time:calendar = "gregorian" ;...