[MITgcm-support] Initial theta field in tutorial (Global Ocean Simulation at 4 deg)

[Martin Losch]

Hi Krishnakumar,

you are right in that lev_t.bin contains 12 monthly averages of temperature  (from an old  Levitus data base, I believe the 1994-version). In ini_theta.F, the model reads only the first 90x40x15 values, corresponding to the January mean, and uses them as the initial conditions. The other 11 months are not use. I agree that this is a little misleading. You could also provide a file that contains just one month (with size 90x40x15), and the initialisation will be fine.

Martin

On Nov 15, 2011, at 5:42 AM, Krishnakumar Rajagopalan wrote:

> Hi All,
>  
> In the tutorial, tutorial_global_oce_latlon (Global Ocean Simulation at 4 deg), the initial theta field file, lev_t.bin has size 90x40x15x12. It appears that after reading the complete field - theta(90,40,15,12), the code assigns theta(:,:,:,1) to the initial theta field and the other 11 fields of theta - theta(:,:,:,2:12) are discarded. Is this correct? Also, is theta(:,:,:,1) the January theta field from Levitas database?
>  
> Thanks a lot
>  
> Best regards
>  
> Krishnakumar
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support