Changeset 218 for LMDZ.3.3/trunk


Ignore:
Timestamp:
May 14, 2001, 5:01:36 PM (24 years ago)
Author:
lmdz
Message:

Generalisation des if NC_DOUBLE pour passer sur VPP MAF
LF

Location:
LMDZ.3.3/trunk/libf/dyn3d
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/trunk/libf/dyn3d/grilles_gcm_netcdf.F

    r207 r218  
    168168c rajoute l'ecriture de la grille
    169169
    170         status=NF_PUT_VARA_REAL(ncid_out,out_lonuid,1,iim+1,rlonudeg)
    171         status=NF_PUT_VARA_REAL(ncid_out,out_lonvid,1,iim+1,rlonvdeg)
    172         status=NF_PUT_VARA_REAL(ncid_out,out_latuid,1,jjm+1,rlatudeg)
    173         status=NF_PUT_VARA_REAL(ncid_out,out_latvid,1,jjm,rlatvdeg)
     170#ifdef NC_DOUBLE
     171      status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonuid,1,iim+1,rlonudeg)
     172      status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonvid,1,iim+1,rlonvdeg)
     173      status=NF_PUT_VARA_DOUBLE(ncid_out,out_latuid,1,jjm+1,rlatudeg)
     174      status=NF_PUT_VARA_DOUBLE(ncid_out,out_latvid,1,jjm,rlatvdeg)
     175#else
     176      status=NF_PUT_VARA_REAL(ncid_out,out_lonuid,1,iim+1,rlonudeg)
     177      status=NF_PUT_VARA_REAL(ncid_out,out_lonvid,1,iim+1,rlonvdeg)
     178      status=NF_PUT_VARA_REAL(ncid_out,out_latuid,1,jjm+1,rlatudeg)
     179      status=NF_PUT_VARA_REAL(ncid_out,out_latvid,1,jjm,rlatvdeg)
     180#endif
    174181
    175182        start(1)=1
     
    189196        enddo
    190197
     198#ifdef NC_DOUBLE
     199        status=NF_PUT_VARA_DOUBLE(ncid_out,out_varid,start,
     200     s  count,temp)
     201#else
    191202        status=NF_PUT_VARA_REAL(ncid_out,out_varid,start,
    192203     s  count,temp)
     204#endif
     205
    193206
    194207c fermeture du fichier netcdf
  • LMDZ.3.3/trunk/libf/dyn3d/read_dstoke.F

    r207 r218  
    4242        debut(2) = 1
    4343
     44#ifdef NC_DOUBLE
     45        status=NF_GET_VARA_DOUBLE(ncidfd,variddt,debut,epais,dtv)
     46#else
    4447        status=NF_GET_VARA_REAL(ncidfd,variddt,debut,epais,dtv)
     48#endif
    4549        zdtvr=dtv(1,1,1)
    4650
     51#ifdef NC_DOUBLE
     52        status=NF_GET_VARA_DOUBLE(ncidfd,varididvt,debut,epais,adv)
     53#else
    4754        status=NF_GET_VARA_REAL(ncidfd,varididvt,debut,epais,adv)
     55#endif
    4856        ziadvtrac= adv(1,1,1)
    4957       
     58#ifdef NC_DOUBLE
     59        status=NF_GET_VARA_DOUBLE(ncidfd,varididvp,debut,epais,adv2)
     60#else
    5061        status=NF_GET_VARA_REAL(ncidfd,varididvp,debut,epais,adv2)
     62#endif
    5163        ziadvtrac2= adv2(1,1,1)
    5264
  • LMDZ.3.3/trunk/libf/dyn3d/read_fstoke.F

    r207 r218  
    2828     . rlatu_dy(imo1,jmo1),rlatv_dy(imo1,jmo)
    2929      integer*4 ncrec,ncim,ncjm,nclm
    30       real zrec,zim,zjm,zlm
     30      integer*4 zrec,zim,zjm,zlm
    3131      integer*4 xid,yid,zid,tid
    3232      real*4 zdtvr,ziadvtrac
     
    126126c niveaux de pression
    127127       
     128#ifdef NC_DOUBLE
     129       status=NF_GET_VARA_DOUBLE(ncidf,varidpl,1,zlm,pl)
     130#else
    128131       status=NF_GET_VARA_REAL(ncidf,varidpl,1,zlm,pl)
     132#endif
    129133
    130134c rlonu,rlonv,rlatu,rlatv
     
    140144      count(4)=0
    141145       
    142         status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy)
    143         status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy)
     146#ifdef NC_DOUBLE
     147      status=NF_GET_VARA_DOUBLE(ncidf,varidnlo,start,count,rlonu_dy)
     148      status=NF_GET_VARA_DOUBLE(ncidf,varidnla,start,count,rlatu_dy)
     149#else
     150      status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy)
     151      status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy)
     152#endif
    144153       
    145154      start(1)=1
     
    153162      count(4)=0
    154163
    155         status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy)
    156         status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy)
     164#ifdef NC_DOUBLE
     165      status=NF_GET_VARA_DOUBLE(ncidfv,varidnlov,start,count,rlonv_dy)
     166      status=NF_GET_VARA_DOUBLE(ncidfv,varidnlav,start,count,rlatv_dy)
     167#else
     168      status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy)
     169      status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy)
     170#endif
    157171
    158172c Lecture de phis et aire
     
    169183
    170184c phis
     185#ifdef NC_DOUBLE
     186      status=NF_GET_VARA_DOUBLE(ncidf,varidps,start,count,phis)
     187#else
    171188      status=NF_GET_VARA_REAL(ncidf,varidps,start,count,phis)
    172 c      print*,'WARNING!!! Correction bidon pour palier a un '
    173 c      print*,'probleme dans la creation des fichiers nc'
    174 c      call correctbid(iim,jmo1*1,phis)
    175        call dump2d(imo1,jmo1,phis,'PHISNC')
     189#endif
    176190
    177191c aire
     192#ifdef NC_DOUBLE
     193       status=NF_GET_VARA_DOUBLE(ncidf,varidai,start,count,airedy)
     194#else
    178195       status=NF_GET_VARA_REAL(ncidf,varidai,start,count,airedy)
    179 c      call correctbid(iim,jmo1*1,airedy)
    180 c       call dump2d(imo1,jmo1,airedy,'AIRENC')
    181 
    182 c        status= NF_CLOSE(ncidf)
     196#endif
    183197
    184198       else
    185199
    186200      print*,'ok1'
    187 c       stop'TEST lectfluxnc'
    188201
    189202c ---------------------
     
    205218
    206219c masse
     220#ifdef NC_DOUBLE
     221      status=NF_GET_VARA_DOUBLE(ncidf,varidm,start,count,masse)
     222#else
    207223      status=NF_GET_VARA_REAL(ncidf,varidm,start,count,masse)
     224#endif
    208225c      print*,'WARNING!!! Correction bidon pour palier a un '
    209226c      print*,'probleme dans la creation des fichiers nc'
     
    212229       call dump2d(imo1,jmo1,masse,'MASSE COUCHE 1')
    213230c pbaru
     231#ifdef NC_DOUBLE
     232      status=NF_GET_VARA_DOUBLE(ncidf,varidpu,start,count,pbaru)
     233#else
    214234      status=NF_GET_VARA_REAL(ncidf,varidpu,start,count,pbaru)
     235#endif
    215236c      print*,'WARNING!!! Correction bidon pour palier a un '
    216237c      print*,'probleme dans la creation des fichiers nc'
     
    219240
    220241c w
     242#ifdef NC_DOUBLE
     243      status=NF_GET_VARA_DOUBLE(ncidf,varidw,start,count,w)
     244#else
    221245      status=NF_GET_VARA_REAL(ncidf,varidw,start,count,w)
     246#endif
    222247c      print*,'WARNING!!! Correction bidon pour palier a un '
    223248c      print*,'probleme dans la creation des fichiers nc'
     
    226251
    227252c teta
     253#ifdef NC_DOUBLE
     254      status=NF_GET_VARA_DOUBLE(ncidf,varidt,start,count,teta)
     255#else
    228256      status=NF_GET_VARA_REAL(ncidf,varidt,start,count,teta)
     257#endif
    229258c      print*,'WARNING!!! Correction bidon pour palier a un '
    230259c      print*,'probleme dans la creation des fichiers nc'
     
    233262
    234263c phi
     264#ifdef NC_DOUBLE
     265      status=NF_GET_VARA_DOUBLE(ncidf,varidph,start,count,phi)
     266#else
    235267      status=NF_GET_VARA_REAL(ncidf,varidph,start,count,phi)
     268#endif
    236269c      print*,'WARNING!!! Correction bidon pour palier a un '
    237270c      print*,'probleme dans la creation des fichiers nc'
     
    243276c  pbarv
    244277        print*,'Lecture de pbarv irec=',irec
     278#ifdef NC_DOUBLE
     279      status=NF_GET_VARA_DOUBLE(ncidfv,varidpv,start,count,pbarv)
     280#else
    245281      status=NF_GET_VARA_REAL(ncidfv,varidpv,start,count,pbarv)
     282#endif
    246283c      call correctbid(iim,jjm*nlevnc,pbarv)
    247284        PRINT*,'*******LECTURE DE PBARV******'
  • LMDZ.3.3/trunk/libf/dyn3d/read_fstoke0.F

    r207 r218  
    2424     . rlatu_dy(iip1,jjp1),rlatv_dy(iip1,jjm)
    2525      integer*4 ncrec,ncim,ncjm,nclm
    26       real zrec,zim,zjm,zlm
     26      integer*4 zrec,zim,zjm,zlm
    2727      integer*4 xid,yid,zid,tid
    2828      real*4 zdtvr,ziadvtrac
     
    122122c niveaux de pression
    123123       
     124#ifdef NC_DOUBLE
     125       status=NF_GET_VARA_DOUBLE(ncidf,varidpl,1,zlm,pl)
     126#else
    124127       status=NF_GET_VARA_REAL(ncidf,varidpl,1,zlm,pl)
     128#endif
    125129
    126130c rlonu,rlonv,rlatu,rlatv
     
    136140      count(4)=0
    137141       
    138         status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy)
    139         status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy)
    140        
     142#ifdef NC_DOUBLE
     143      status=NF_GET_VARA_DOUBLE(ncidf,varidnlo,start,count,rlonu_dy)
     144      status=NF_GET_VARA_DOUBLE(ncidf,varidnla,start,count,rlatu_dy)
     145#else
     146      status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy)
     147      status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy)
     148#endif
     149
    141150      start(1)=1
    142151      start(2)=1
     
    149158      count(4)=0
    150159
    151         status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy)
    152         status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy)
     160#ifdef NC_DOUBLE
     161      status=NF_GET_VARA_DOUBLE(ncidfv,varidnlov,start,count,rlonv_dy)
     162      status=NF_GET_VARA_DOUBLE(ncidfv,varidnlav,start,count,rlatv_dy)
     163#else
     164      status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy)
     165      status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy)
     166#endif
    153167
    154168c Lecture de phis et aire
     
    165179
    166180c phis
     181#ifdef NC_DOUBLE
     182      status=NF_GET_VARA_DOUBLE(ncidf,varidps,start,count,phis)
     183#else
    167184      status=NF_GET_VARA_REAL(ncidf,varidps,start,count,phis)
     185#endif
    168186
    169187c aire
     188#ifdef NC_DOUBLE
     189       status=NF_GET_VARA_DOUBLE(ncidf,varidai,start,count,airedy)
     190#else
    170191       status=NF_GET_VARA_REAL(ncidf,varidai,start,count,airedy)
     192#endif
    171193
    172194       else
     
    193215c masse
    194216        PRINT*,'LECTURE DE masse'
     217#ifdef NC_DOUBLE
     218      status=NF_GET_VARA_DOUBLE(ncidf,varidm,start,count,masse)
     219#else
    195220      status=NF_GET_VARA_REAL(ncidf,varidm,start,count,masse)
     221#endif
    196222
    197223c pbaru
     224#ifdef NC_DOUBLE
     225      status=NF_GET_VARA_DOUBLE(ncidf,varidpu,start,count,pbaru)
     226#else
    198227      status=NF_GET_VARA_REAL(ncidf,varidpu,start,count,pbaru)
     228#endif
    199229
    200230c w
     231#ifdef NC_DOUBLE
     232      status=NF_GET_VARA_DOUBLE(ncidf,varidw,start,count,w)
     233#else
    201234      status=NF_GET_VARA_REAL(ncidf,varidw,start,count,w)
     235#endif
    202236
    203237c teta
     238#ifdef NC_DOUBLE
     239      status=NF_GET_VARA_DOUBLE(ncidf,varidt,start,count,teta)
     240#else
    204241      status=NF_GET_VARA_REAL(ncidf,varidt,start,count,teta)
     242#endif
    205243
    206244c phi
     245#ifdef NC_DOUBLE
     246      status=NF_GET_VARA_DOUBLE(ncidf,varidph,start,count,phi)
     247#else
    207248      status=NF_GET_VARA_REAL(ncidf,varidph,start,count,phi)
     249#endif
    208250       
    209251        count(2) = zjm-1
    210252
    211253c  pbarv
     254#ifdef NC_DOUBLE
     255      status=NF_GET_VARA_DOUBLE(ncidfv,varidpv,start,count,pbarv)
     256#else
    212257      status=NF_GET_VARA_REAL(ncidfv,varidpv,start,count,pbarv)
     258#endif
    213259
    214260      start(3)=irec
     
    222268      return
    223269      end
    224 
  • LMDZ.3.3/trunk/libf/dyn3d/read_pstoke.F

    r185 r218  
    196196c niveaux de pression
    197197
     198#ifdef NC_DOUBLE
     199            status=NF_GET_VARA_DOUBLE(ncidp,varidpl,1,kev,pl)
     200#else
    198201            status=NF_GET_VARA_REAL(ncidp,varidpl,1,kev,pl)
     202#endif
    199203
    200204c lecture de aire et phis
     
    212216c
    213217c phis
     218#ifdef NC_DOUBLE
     219      status=NF_GET_VARA_DOUBLE(ncidp,varidps,start,count,phisfi2)
     220#else
    214221      status=NF_GET_VARA_REAL(ncidp,varidps,start,count,phisfi2)
     222#endif
     223
     224
    215225      call gr_ecrit_fi(1,kon,iim,jjm+1,phisfi2,phisfi)
    216226
    217227c aire
     228#ifdef NC_DOUBLE
     229      status=NF_GET_VARA_DOUBLE(ncidp,varidai,start,count,airefi2)
     230#else
    218231      status=NF_GET_VARA_REAL(ncidp,varidai,start,count,airefi2)
     232#endif
     233
    219234      call gr_ecrit_fi(1,kon,iim,jjm+1,airefi2,airefi)
    220235      else
     
    241256c frac_impa
    242257
     258#ifdef NC_DOUBLE
     259      status=NF_GET_VARA_DOUBLE(ncidp,varidfi,start,count,frac_impa2)
     260#else
    243261      status=NF_GET_VARA_REAL(ncidp,varidfi,start,count,frac_impa2)
     262#endif
    244263      call gr_ecrit_fi(kev,kon,iim,jjm+1,frac_impa2,frac_impa)
    245264
    246265c frac_nucl
    247266
     267#ifdef NC_DOUBLE
     268      status=NF_GET_VARA_DOUBLE(ncidp,varidfn,start,count,frac_nucl2)
     269#else
    248270      status=NF_GET_VARA_REAL(ncidp,varidfn,start,count,frac_nucl2)
     271#endif
    249272      call gr_ecrit_fi(kev,kon,iim,jjm+1,frac_nucl2,frac_nucl)
    250273
    251274c abder t
     275
     276#ifdef NC_DOUBLE
     277      status=NF_GET_VARA_DOUBLE(ncidp,varidt,start,count,t2)
     278#else
    252279      status=NF_GET_VARA_REAL(ncidp,varidt,start,count,t2)
     280#endif
     281
    253282      call gr_ecrit_fi(kev,kon,iim,jjm+1,t2,t)
    254283
    255284c mfu
     285
     286#ifdef NC_DOUBLE
     287      status=NF_GET_VARA_DOUBLE(ncidp,varidmfu,start,count,mfu2)
     288#else
    256289      status=NF_GET_VARA_REAL(ncidp,varidmfu,start,count,mfu2)
     290#endif
     291
    257292      call gr_ecrit_fi(kev,kon,iim,jjm+1,mfu2,mfu)
    258293
    259294c mfd
     295#ifdef NC_DOUBLE
     296      status=NF_GET_VARA_DOUBLE(ncidp,varidmfd,start,count,mfd2)
     297#else
    260298      status=NF_GET_VARA_REAL(ncidp,varidmfd,start,count,mfd2)
     299#endif
    261300      call gr_ecrit_fi(kev,kon,iim,jjm+1,mfd2,mfd)
    262301
    263302c en_u
     303#ifdef NC_DOUBLE
     304      status=NF_GET_VARA_DOUBLE(ncidp,varidenu,start,count,en_u2)
     305#else
    264306      status=NF_GET_VARA_REAL(ncidp,varidenu,start,count,en_u2)
     307#endif
    265308      call gr_ecrit_fi(kev,kon,iim,jjm+1,en_u2,en_u)
    266309
    267310c de_u
     311#ifdef NC_DOUBLE
     312      status=NF_GET_VARA_DOUBLE(ncidp,variddeu,start,count,de_u2)
     313#else
    268314      status=NF_GET_VARA_REAL(ncidp,variddeu,start,count,de_u2)
     315#endif
    269316      call gr_ecrit_fi(kev,kon,iim,jjm+1,de_u2,de_u)
    270317
    271318c en_d
     319#ifdef NC_DOUBLE
     320      status=NF_GET_VARA_DOUBLE(ncidp,varidend,start,count,en_d2)
     321#else
    272322      status=NF_GET_VARA_REAL(ncidp,varidend,start,count,en_d2)
     323#endif
    273324      call gr_ecrit_fi(kev,kon,iim,jjm+1,en_d2,en_d)
    274325
    275326c de_d
     327#ifdef NC_DOUBLE
     328      status=NF_GET_VARA_DOUBLE(ncidp,varidded,start,count,de_d2)
     329#else
    276330      status=NF_GET_VARA_REAL(ncidp,varidded,start,count,de_d2)
     331#endif
    277332      call gr_ecrit_fi(kev,kon,iim,jjm+1,de_d2,de_d)
    278333
    279334c coefh
    280335        print*,'LECTURE de coefh a irec =',irec
     336#ifdef NC_DOUBLE
     337       status=NF_GET_VARA_DOUBLE(ncidp,varidch,start,count,coefh2)
     338#else
    281339       status=NF_GET_VARA_REAL(ncidp,varidch,start,count,coefh2)
     340#endif
     341
    282342       call gr_ecrit_fi(kev,kon,iim,jjm+1,coefh2,coefh)
    283343
     
    289349c pyu1
    290350        print*,'LECTURE de yu1 a irec =',irec
     351#ifdef NC_DOUBLE
     352      status=NF_GET_VARA_DOUBLE(ncidp,varidyu1,start,count,pyu12)
     353#else
    291354      status=NF_GET_VARA_REAL(ncidp,varidyu1,start,count,pyu12)
     355#endif
     356
    292357      call gr_ecrit_fi(1,kon,iim,jjm+1,pyu12,pyu1)
    293358
    294359c pyv1
    295360        print*,'LECTURE de yv1 a irec =',irec
     361#ifdef NC_DOUBLE
     362      status=NF_GET_VARA_DOUBLE(ncidp,varidyv1,start,count,pyv12)
     363#else
    296364      status=NF_GET_VARA_REAL(ncidp,varidyv1,start,count,pyv12)
     365#endif
    297366      call gr_ecrit_fi(1,kon,iim,jjm+1,pyv12,pyv1)
    298367
    299368c ftsol1
    300369        print*,'LECTURE de ftsol1 a irec =',irec
     370#ifdef NC_DOUBLE
     371      status=NF_GET_VARA_DOUBLE(ncidp,varidfts1,start,count,ftsol12)
     372#else
    301373      status=NF_GET_VARA_REAL(ncidp,varidfts1,start,count,ftsol12)
     374#endif
    302375       call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol12,ftsol1)
    303376
    304377c ftsol2
    305378        print*,'LECTURE de ftsol2 a irec =',irec
     379#ifdef NC_DOUBLE
     380      status=NF_GET_VARA_DOUBLE(ncidp,varidfts2,start,count,ftsol22)
     381#else
    306382      status=NF_GET_VARA_REAL(ncidp,varidfts2,start,count,ftsol22)
     383#endif
    307384      call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol22,ftsol2)
    308385
    309386c ftsol3
    310387         print*,'LECTURE de ftsol3 a irec =',irec
     388#ifdef NC_DOUBLE
     389      status=NF_GET_VARA_DOUBLE(ncidp,varidfts3,start,count,ftsol32)
     390#else
    311391      status=NF_GET_VARA_REAL(ncidp,varidfts3,start,count,ftsol32)
     392#endif
    312393      call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol32,ftsol3)
    313394
    314395c ftsol4
     396#ifdef NC_DOUBLE
     397      status=NF_GET_VARA_DOUBLE(ncidp,varidfts4,start,count,ftsol42)
     398#else
    315399      status=NF_GET_VARA_REAL(ncidp,varidfts4,start,count,ftsol42)
     400#endif
    316401      call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol42,ftsol4)
    317402
    318403c psrf1
     404#ifdef NC_DOUBLE
     405      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr1,start,count,psrf12)
     406#else
    319407      status=NF_GET_VARA_REAL(ncidp,varidpsr1,start,count,psrf12)
     408#endif
    320409c      call dump2d(iip1-1,jjm+1,psrf12,'PSRF1NC')
    321410      call gr_ecrit_fi(1,kon,iim,jjm+1,psrf12,psrf1)
    322411
    323412c psrf2
     413#ifdef NC_DOUBLE
     414      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr2,start,count,psrf22)
     415#else
    324416      status=NF_GET_VARA_REAL(ncidp,varidpsr2,start,count,psrf22)
     417#endif
    325418c      call dump2d(iip1-1,jjm+1,psrf22,'PSRF2NC')
    326419      call gr_ecrit_fi(1,kon,iim,jjm+1,psrf22,psrf2)
    327420
    328421c psrf3
     422#ifdef NC_DOUBLE
     423      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr3,start,count,psrf32)
     424#else
    329425      status=NF_GET_VARA_REAL(ncidp,varidpsr3,start,count,psrf32)
     426#endif
    330427      call gr_ecrit_fi(1,kon,iim,jjm+1,psrf32,psrf3)
    331428
    332429c psrf4
     430#ifdef NC_DOUBLE
     431      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr4,start,count,psrf42)
     432#else
    333433      status=NF_GET_VARA_REAL(ncidp,varidpsr4,start,count,psrf42)
     434#endif
    334435      call gr_ecrit_fi(1,kon,iim,jjm+1,psrf42,psrf4)
    335436       
Note: See TracChangeset for help on using the changeset viewer.