[MITgcm-support] OBCS checkpoint

[Gianmaria Sannino]

Hello,

I'm using version checkpoint52b_post of the MITgcm in a quite simple
configuration (115x62x31) on a single processor. The model has
activated only a sponge obc on the western side.

In this configuration there is a problem with the obcs checkpoint. If you
look at the obcs_checkpoint.F file you can see:

#ifdef ALLOW_OBCS_EAST
      CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBEu,1,myIt,myThid)
      CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBEv,2,myIt,myThid)
      CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBEt,3,myIt,myThid)
      CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBEs,4,myIt,myThid)
# ifdef ALLOW_NONHYDROSTATIC
      CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBEw,9,myIt,myThid)
# endif
#endif


#ifdef ALLOW_OBCS_WEST
       CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBWu,5,myIt,myThid)
       CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBWv,6,myIt,myThid)
       CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBWt,7,myIt,myThid)
       CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBWs,8,myIt,myThid)
# ifdef ALLOW_NONHYDROSTATIC
      CALL MDSWRITEFIELDYZ(fn,prec,lgf,'RS',Nr,OBWw,10,myIt,myThid)
# endif


This means that in my case the first record of the pickup_obEW file is 5,
but this is a problem for the Subrotuine MDSWRITEFIELDYZ that is waiting a
record number 1 to open the pickup_obEW file as a new file. Probably this
issue can be easly resolved using four pichup_ob files, i.e. one for each
open side.


Thanks
gianmaria




******************************************************
Gianmaria Sannino
Climate Project - Ocean Modelling Unit 
SP 91 - via Anguillarese 301
S.M. di Galeria, 00060 ROMA, ITALY
Voice: (+39) 06 3048 6799   Fax: (+39) 06 3048 4264
e-mail: gianmaria.sannino at casaccia.enea.it
******************************************************