Changeset 256 for LMDZ.3.3/trunk


Ignore:
Timestamp:
Jul 17, 2001, 5:34:41 PM (23 years ago)
Author:
lmdz
Message:

Recuperation de la version qui se trouvait dans dyn3d et remise des NC_DOUBLE
pour VPP. MAF
LF

Location:
LMDZ.3.3/trunk/libf/phylmd
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/trunk/libf/phylmd/read_pstoke.F

    r207 r256  
    199199c niveaux de pression
    200200
     201#ifdef NC_DOUBLE
     202            status=NF_GET_VARA_DOUBLE(ncidp,varidpl,1,zklevo,pl)
     203#else
    201204            status=NF_GET_VARA_REAL(ncidp,varidpl,1,zklevo,pl)
     205#endif
     206
    202207
    203208c lecture de aire et phis
    … …  
    214219
    215220c phis
     221#ifdef NC_DOUBLE
     222      status=NF_GET_VARA_DOUBLE(ncidp,varidps,start,count,phisfi2)
     223#else
    216224      status=NF_GET_VARA_REAL(ncidp,varidps,start,count,phisfi2)
     225#endif
    217226c      print*,'WARNING!!! Correction bidon pour palier a un '
    218227c      print*,'probleme dans la creation des fichiers nc'
    … …  
    222231
    223232c aire
     233#ifdef NC_DOUBLE
     234      status=NF_GET_VARA_DOUBLE(ncidp,varidai,start,count,airefi2)
     235#else
    224236      status=NF_GET_VARA_REAL(ncidp,varidai,start,count,airefi2)
     237#endif
     238
    225239c      call correctbid(iim,jjp1*1,airefi2)
    226240c       call dump2d(iip1-1,jjp1,airefi2,'AIRENC')
    … …  
    249263c frac_impa
    250264
     265#ifdef NC_DOUBLE
     266      status=NF_GET_VARA_DOUBLE(ncidp,varidfi,start,count,frac_impa2)
     267#else
    251268      status=NF_GET_VARA_REAL(ncidp,varidfi,start,count,frac_impa2)
     269#endif
    252270c      print*,'WARNING!!! Correction bidon pour palier a un '
    253271c      print*,'probleme dans la creation des fichiers nc'
    … …  
    258276c frac_nucl
    259277
     278#ifdef NC_DOUBLE
     279      status=NF_GET_VARA_DOUBLE(ncidp,varidfn,start,count,frac_nucl2)
     280#else
    260281      status=NF_GET_VARA_REAL(ncidp,varidfn,start,count,frac_nucl2)
     282#endif
     283
    261284c      print*,'WARNING!!! Correction bidon pour palier a un '
    262285c      print*,'probleme dans la creation des fichiers nc'
    … …  
    266289
    267290c abder t
     291#ifdef NC_DOUBLE
     292      status=NF_GET_VARA_DOUBLE(ncidp,varidt,start,count,t2)
     293#else
    268294      status=NF_GET_VARA_REAL(ncidp,varidt,start,count,t2)
     295#endif
    269296      call gr_ecrit_fi(klevo,klono,imo,jmo+1,t2,t)
    270297
    271298c mfu
     299#ifdef NC_DOUBLE
     300      status=NF_GET_VARA_DOUBLE(ncidp,varidmfu,start,count,mfu2)
     301#else
    272302      status=NF_GET_VARA_REAL(ncidp,varidmfu,start,count,mfu2)
     303#endif
     304
    273305c      print*,'WARNING!!! Correction bidon pour palier a un '
    274306c      print*,'probleme dans la creation des fichiers nc'
    … …  
    278310
    279311c mfd
     312#ifdef NC_DOUBLE
     313      status=NF_GET_VARA_DOUBLE(ncidp,varidmfd,start,count,mfd2)
     314#else
    280315      status=NF_GET_VARA_REAL(ncidp,varidmfd,start,count,mfd2)
     316#endif
    281317c      print*,'WARNING!!! Correction bidon pour palier a un '
    282318c      print*,'probleme dans la creation des fichiers nc'
    … …  
    286322
    287323c en_u
     324#ifdef NC_DOUBLE
     325      status=NF_GET_VARA_DOUBLE(ncidp,varidenu,start,count,en_u2)
     326#else
    288327      status=NF_GET_VARA_REAL(ncidp,varidenu,start,count,en_u2)
     328#endif
    289329c      print*,'WARNING!!! Correction bidon pour palier a un '
    290330c      print*,'probleme dans la creation des fichiers nc'
    … …  
    294334
    295335c de_u
     336#ifdef NC_DOUBLE
     337      status=NF_GET_VARA_DOUBLE(ncidp,variddeu,start,count,de_u2)
     338#else
    296339      status=NF_GET_VARA_REAL(ncidp,variddeu,start,count,de_u2)
     340#endif
     341
    297342c      print*,'WARNING!!! Correction bidon pour palier a un '
    298343c      print*,'probleme dans la creation des fichiers nc'
    … …  
    302347
    303348c en_d
     349#ifdef NC_DOUBLE
     350      status=NF_GET_VARA_DOUBLE(ncidp,varidend,start,count,en_d2)
     351#else
    304352      status=NF_GET_VARA_REAL(ncidp,varidend,start,count,en_d2)
     353#endif
    305354c      print*,'WARNING!!! Correction bidon pour palier a un '
    306355c      print*,'probleme dans la creation des fichiers nc'
    … …  
    310359
    311360c de_d
     361#ifdef NC_DOUBLE
     362      status=NF_GET_VARA_DOUBLE(ncidp,varidded,start,count,de_d2)
     363#else
    312364      status=NF_GET_VARA_REAL(ncidp,varidded,start,count,de_d2)
     365#endif
    313366c      print*,'WARNING!!! Correction bidon pour palier a un '
    314367c      print*,'probleme dans la creation des fichiers nc'
    … …  
    318371
    319372c coefh
    320       status=NF_GET_VARA_REAL(ncidp,varidch,start,count,coefh2)
     373        print*,'LECTURE de coefh a irec =',irec
     374#ifdef NC_DOUBLE
     375       status=NF_GET_VARA_DOUBLE(ncidp,varidch,start,count,coefh2)
     376#else
     377       status=NF_GET_VARA_REAL(ncidp,varidch,start,count,coefh2)
     378#endif
    321379c      print*,'WARNING!!! Correction bidon pour palier a un '
    322380c      print*,'probleme dans la creation des fichiers nc'
    … …  
    339397
    340398c pyv1
     399#ifdef NC_DOUBLE
     400      status=NF_GET_VARA_DOUBLE(ncidp,varidyv1,start,count,pyv12)
     401#else
    341402      status=NF_GET_VARA_REAL(ncidp,varidyv1,start,count,pyv12)
     403#endif
    342404c      print*,'WARNING!!! Correction bidon pour palier a un '
    343405c      print*,'probleme dans la creation des fichiers nc'
    … …  
    347409
    348410c ftsol1
     411#ifdef NC_DOUBLE
     412      status=NF_GET_VARA_DOUBLE(ncidp,varidfts1,start,count,ftsol12)
     413#else
    349414      status=NF_GET_VARA_REAL(ncidp,varidfts1,start,count,ftsol12)
     415#endif
    350416c      print*,'WARNING!!! Correction bidon pour palier a un '
    351417c      print*,'probleme dans la creation des fichiers nc'
    … …  
    355421
    356422c ftsol2
     423#ifdef NC_DOUBLE
     424      status=NF_GET_VARA_DOUBLE(ncidp,varidfts2,start,count,ftsol22)
     425#else
    357426      status=NF_GET_VARA_REAL(ncidp,varidfts2,start,count,ftsol22)
     427#endif
    358428c      print*,'WARNING!!! Correction bidon pour palier a un '
    359429c      print*,'probleme dans la creation des fichiers nc'
    … …  
    363433
    364434c ftsol3
     435#ifdef NC_DOUBLE
     436      status=NF_GET_VARA_DOUBLE(ncidp,varidfts3,start,count,ftsol32)
     437#else
    365438      status=NF_GET_VARA_REAL(ncidp,varidfts3,start,count,ftsol32)
     439#endif
    366440c      print*,'WARNING!!! Correction bidon pour palier a un '
    367441c      print*,'probleme dans la creation des fichiers nc'
    … …  
    371445
    372446c ftsol4
     447#ifdef NC_DOUBLE
     448      status=NF_GET_VARA_DOUBLE(ncidp,varidfts4,start,count,ftsol42)
     449#else
    373450      status=NF_GET_VARA_REAL(ncidp,varidfts4,start,count,ftsol42)
     451#endif
    374452c      print*,'WARNING!!! Correction bidon pour palier a un '
    375453c      print*,'probleme dans la creation des fichiers nc'
    … …  
    379457
    380458c psrf1
    381       status=NF_GET_VARA_REAL(ncidp,varidpsr1,start,count,psrf12)
     459#ifdef NC_DOUBLE
     460      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr2,start,count,psrf22)
     461#else
     462      status=NF_GET_VARA_REAL(ncidp,varidpsr2,start,count,psrf22)
     463#endif
    382464c      print*,'WARNING!!! Correction bidon pour palier a un '
    383465c      print*,'probleme dans la creation des fichiers nc'
    … …  
    387469
    388470c psrf2
     471#ifdef NC_DOUBLE
     472      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr2,start,count,psrf22)
     473#else
    389474      status=NF_GET_VARA_REAL(ncidp,varidpsr2,start,count,psrf22)
     475#endif
    390476c      print*,'WARNING!!! Correction bidon pour palier a un '
    391477c      print*,'probleme dans la creation des fichiers nc'
    … …  
    395481
    396482c psrf3
     483#ifdef NC_DOUBLE
     484      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr3,start,count,psrf32)
     485#else
    397486      status=NF_GET_VARA_REAL(ncidp,varidpsr3,start,count,psrf32)
     487#endif
    398488c      print*,'WARNING!!! Correction bidon pour palier a un '
    399489c      print*,'probleme dans la creation des fichiers nc'
    … …  
    403493
    404494c psrf4
     495#ifdef NC_DOUBLE
     496      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr4,start,count,psrf42)
     497#else
    405498      status=NF_GET_VARA_REAL(ncidp,varidpsr4,start,count,psrf42)
     499#endif
    406500c      print*,'WARNING!!! Correction bidon pour palier a un '
    407501c      print*,'probleme dans la creation des fichiers nc'

  • LMDZ.3.3/trunk/libf/phylmd/read_pstoke0.F

    r207 r256  
    1 c
    2 c $Header$
    3 c
    41        subroutine read_pstoke0(irec,
    52     .   zrec,zkon,zkev,airefi,phisfi,
    … …  
    199196c niveaux de pression
    200197
     198#ifdef NC_DOUBLE
     199            status=NF_GET_VARA_DOUBLE(ncidp,varidpl,1,kev,pl)
     200#else
    201201            status=NF_GET_VARA_REAL(ncidp,varidpl,1,kev,pl)
     202#endif
    202203
    203204c lecture de aire et phis
    … …  
    215216c
    216217c phis
     218#ifdef NC_DOUBLE
     219      status=NF_GET_VARA_DOUBLE(ncidp,varidps,start,count,phisfi2)
     220#else
    217221      status=NF_GET_VARA_REAL(ncidp,varidps,start,count,phisfi2)
     222#endif
     223
     224
    218225      call gr_ecrit_fi(1,kon,iim,jjm+1,phisfi2,phisfi)
    219226
    220227c aire
     228#ifdef NC_DOUBLE
     229      status=NF_GET_VARA_DOUBLE(ncidp,varidai,start,count,airefi2)
     230#else
    221231      status=NF_GET_VARA_REAL(ncidp,varidai,start,count,airefi2)
     232#endif
     233
    222234      call gr_ecrit_fi(1,kon,iim,jjm+1,airefi2,airefi)
    223235      else
    … …  
    244256c frac_impa
    245257
     258#ifdef NC_DOUBLE
     259      status=NF_GET_VARA_DOUBLE(ncidp,varidfi,start,count,frac_impa2)
     260#else
    246261      status=NF_GET_VARA_REAL(ncidp,varidfi,start,count,frac_impa2)
     262#endif
    247263      call gr_ecrit_fi(kev,kon,iim,jjm+1,frac_impa2,frac_impa)
    248264
    249265c frac_nucl
    250266
     267#ifdef NC_DOUBLE
     268      status=NF_GET_VARA_DOUBLE(ncidp,varidfn,start,count,frac_nucl2)
     269#else
    251270      status=NF_GET_VARA_REAL(ncidp,varidfn,start,count,frac_nucl2)
     271#endif
    252272      call gr_ecrit_fi(kev,kon,iim,jjm+1,frac_nucl2,frac_nucl)
    253273
    254274c abder t
     275
     276#ifdef NC_DOUBLE
     277      status=NF_GET_VARA_DOUBLE(ncidp,varidt,start,count,t2)
     278#else
    255279      status=NF_GET_VARA_REAL(ncidp,varidt,start,count,t2)
     280#endif
     281
    256282      call gr_ecrit_fi(kev,kon,iim,jjm+1,t2,t)
    257283
    258284c mfu
     285
     286#ifdef NC_DOUBLE
     287      status=NF_GET_VARA_DOUBLE(ncidp,varidmfu,start,count,mfu2)
     288#else
    259289      status=NF_GET_VARA_REAL(ncidp,varidmfu,start,count,mfu2)
     290#endif
     291
    260292      call gr_ecrit_fi(kev,kon,iim,jjm+1,mfu2,mfu)
    261293
    262294c mfd
     295#ifdef NC_DOUBLE
     296      status=NF_GET_VARA_DOUBLE(ncidp,varidmfd,start,count,mfd2)
     297#else
    263298      status=NF_GET_VARA_REAL(ncidp,varidmfd,start,count,mfd2)
     299#endif
    264300      call gr_ecrit_fi(kev,kon,iim,jjm+1,mfd2,mfd)
    265301
    266302c en_u
     303#ifdef NC_DOUBLE
     304      status=NF_GET_VARA_DOUBLE(ncidp,varidenu,start,count,en_u2)
     305#else
    267306      status=NF_GET_VARA_REAL(ncidp,varidenu,start,count,en_u2)
     307#endif
    268308      call gr_ecrit_fi(kev,kon,iim,jjm+1,en_u2,en_u)
    269309
    270310c de_u
     311#ifdef NC_DOUBLE
     312      status=NF_GET_VARA_DOUBLE(ncidp,variddeu,start,count,de_u2)
     313#else
    271314      status=NF_GET_VARA_REAL(ncidp,variddeu,start,count,de_u2)
     315#endif
    272316      call gr_ecrit_fi(kev,kon,iim,jjm+1,de_u2,de_u)
    273317
    274318c en_d
     319#ifdef NC_DOUBLE
     320      status=NF_GET_VARA_DOUBLE(ncidp,varidend,start,count,en_d2)
     321#else
    275322      status=NF_GET_VARA_REAL(ncidp,varidend,start,count,en_d2)
     323#endif
    276324      call gr_ecrit_fi(kev,kon,iim,jjm+1,en_d2,en_d)
    277325
    278326c de_d
     327#ifdef NC_DOUBLE
     328      status=NF_GET_VARA_DOUBLE(ncidp,varidded,start,count,de_d2)
     329#else
    279330      status=NF_GET_VARA_REAL(ncidp,varidded,start,count,de_d2)
     331#endif
    280332      call gr_ecrit_fi(kev,kon,iim,jjm+1,de_d2,de_d)
    281333
    282334c coefh
    283335        print*,'LECTURE de coefh a irec =',irec
     336#ifdef NC_DOUBLE
     337       status=NF_GET_VARA_DOUBLE(ncidp,varidch,start,count,coefh2)
     338#else
    284339       status=NF_GET_VARA_REAL(ncidp,varidch,start,count,coefh2)
     340#endif
     341
    285342       call gr_ecrit_fi(kev,kon,iim,jjm+1,coefh2,coefh)
    286343
    … …  
    292349c pyu1
    293350        print*,'LECTURE de yu1 a irec =',irec
     351#ifdef NC_DOUBLE
     352      status=NF_GET_VARA_DOUBLE(ncidp,varidyu1,start,count,pyu12)
     353#else
    294354      status=NF_GET_VARA_REAL(ncidp,varidyu1,start,count,pyu12)
     355#endif
     356
    295357      call gr_ecrit_fi(1,kon,iim,jjm+1,pyu12,pyu1)
    296358
    297359c pyv1
    298360        print*,'LECTURE de yv1 a irec =',irec
     361#ifdef NC_DOUBLE
     362      status=NF_GET_VARA_DOUBLE(ncidp,varidyv1,start,count,pyv12)
     363#else
    299364      status=NF_GET_VARA_REAL(ncidp,varidyv1,start,count,pyv12)
     365#endif
    300366      call gr_ecrit_fi(1,kon,iim,jjm+1,pyv12,pyv1)
    301367
    302368c ftsol1
    303369        print*,'LECTURE de ftsol1 a irec =',irec
     370#ifdef NC_DOUBLE
     371      status=NF_GET_VARA_DOUBLE(ncidp,varidfts1,start,count,ftsol12)
     372#else
    304373      status=NF_GET_VARA_REAL(ncidp,varidfts1,start,count,ftsol12)
     374#endif
    305375       call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol12,ftsol1)
    306376
    307377c ftsol2
    308378        print*,'LECTURE de ftsol2 a irec =',irec
     379#ifdef NC_DOUBLE
     380      status=NF_GET_VARA_DOUBLE(ncidp,varidfts2,start,count,ftsol22)
     381#else
    309382      status=NF_GET_VARA_REAL(ncidp,varidfts2,start,count,ftsol22)
     383#endif
    310384      call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol22,ftsol2)
    311385
    312386c ftsol3
    313387         print*,'LECTURE de ftsol3 a irec =',irec
     388#ifdef NC_DOUBLE
     389      status=NF_GET_VARA_DOUBLE(ncidp,varidfts3,start,count,ftsol32)
     390#else
    314391      status=NF_GET_VARA_REAL(ncidp,varidfts3,start,count,ftsol32)
     392#endif
    315393      call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol32,ftsol3)
    316394
    317395c ftsol4
     396#ifdef NC_DOUBLE
     397      status=NF_GET_VARA_DOUBLE(ncidp,varidfts4,start,count,ftsol42)
     398#else
    318399      status=NF_GET_VARA_REAL(ncidp,varidfts4,start,count,ftsol42)
     400#endif
    319401      call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol42,ftsol4)
    320402
    321403c psrf1
     404#ifdef NC_DOUBLE
     405      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr1,start,count,psrf12)
     406#else
    322407      status=NF_GET_VARA_REAL(ncidp,varidpsr1,start,count,psrf12)
     408#endif
    323409c      call dump2d(iip1-1,jjm+1,psrf12,'PSRF1NC')
    324410      call gr_ecrit_fi(1,kon,iim,jjm+1,psrf12,psrf1)
    325411
    326412c psrf2
     413#ifdef NC_DOUBLE
     414      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr2,start,count,psrf22)
     415#else
    327416      status=NF_GET_VARA_REAL(ncidp,varidpsr2,start,count,psrf22)
     417#endif
    328418c      call dump2d(iip1-1,jjm+1,psrf22,'PSRF2NC')
    329419      call gr_ecrit_fi(1,kon,iim,jjm+1,psrf22,psrf2)
    330420
    331421c psrf3
     422#ifdef NC_DOUBLE
     423      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr3,start,count,psrf32)
     424#else
    332425      status=NF_GET_VARA_REAL(ncidp,varidpsr3,start,count,psrf32)
     426#endif
    333427      call gr_ecrit_fi(1,kon,iim,jjm+1,psrf32,psrf3)
    334428
    335429c psrf4
     430#ifdef NC_DOUBLE
     431      status=NF_GET_VARA_DOUBLE(ncidp,varidpsr4,start,count,psrf42)
     432#else
    336433      status=NF_GET_VARA_REAL(ncidp,varidpsr4,start,count,psrf42)
     434#endif
    337435      call gr_ecrit_fi(1,kon,iim,jjm+1,psrf42,psrf4)
    338436       
Note: See TracChangeset for help on using the changeset viewer.