Ocean Modeling Discussion

that is the problem of multiple forcing files. I encountered the same problem. You need the modify the code in ROMS/Utility/inquiry.F when it CALL netcdf_close (ng, model, ncid). In the case of multiples forcing files, the default value for Lupdate is true. The code will try to update global ...

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 ...

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 ...

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" ;
float SST(sst_time ...