Changeset 352


Ignore:
Timestamp:
Mar 21, 2002, 4:09:30 PM (23 years ago)
Author:
lmdzadmin
Message:

2 changements pour les fichiers histoire:

  • utilisation de l'entree "rectilineaire" de IOIPSL pour ne plus avoir a lancer ncregular a chaque fois
  • le calendrier des fichiers histoire est maintenant base sur la date d'initialisation de la simulation plutot que sur la date de depart du job en cours

LF

Location:
LMDZ.3.3/branches/rel-LF/libf
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/etat0_netcdf.F

    r332 r352  
    2323#include "dimphy.h"
    2424#include "dimsoil.h"
     25#include "temps.h"
    2526      !
    2627      LOGICAL interbar
    … …  
    3839      REAL :: tsoil(klon,nsoilmx,nbsrf)
    3940      REAL :: radsol(klon),rain_fall(klon), snow_fall(klon)
    40       REAL :: solsw(klon), sollw(klon)
    41       REAL :: deltat(klon), frugs(klon,nbsrf), agesno(klon),rugmer(klon)
     41      REAL :: solsw(klon), sollw(klon), fder(klon)
     42      REAL :: deltat(klon), frugs(klon,nbsrf), agesno(klon,nbsrf)
     43      REAL :: rugmer(klon)
    4244      REAL :: zmea(iip1*jjp1), zstd(iip1*jjp1)
    4345      REAL :: zsig(iip1*jjp1), zgam(iip1*jjp1), zthe(iip1*jjp1)
    … …  
    5052      INTEGER :: iml_lic, jml_lic, llm_tmp, ttm_tmp, iret
    5153      INTEGER :: itaul(1), fid
    52       REAL :: lev(1), date, dt
     54      REAL :: lev(1), date
    5355      REAL, ALLOCATABLE, DIMENSION(:,:) :: lon_lic, lat_lic
    5456      REAL, ALLOCATABLE, DIMENSION(:)  :: dlon_lic, dlat_lic
    … …  
    193195      pctsrf=0.
    194196      !cree le masque a partir du fichier relief
    195       varname = 'zmasq'
    196       zmasq(:) = 0.
    197       CALL startget(varname,iip1,jjp1,rlonv,rlatu,klon,zmasq,0.0,
    198      , jjm ,rlonu,rlatv , interbar )
     197      call gr_dyn_fi(1, iip1, jjp1, klon, masque, zmasq)
    199198      WHERE (zmasq(1 : klon) .LT. EPSFRA)
    200199          zmasq(1 : klon) = 0.
    … …  
    323322     .     jjm, rlonu, rlatv , interbar )
    324323      !
    325       varname = 'agsno'
    326       agesno(:) = 0.0
    327       CALL startget(varname,iip1,jjp1,rlonv,rlatu,klon,agesno,0.0,
    328      .     jjm, rlonu, rlatv , interbar )
     324!      varname = 'agesno'
     325!      agesno(:) = 0.0
     326!      CALL startget(varname,iip1,jjp1,rlonv,rlatu,klon,agesno,0.0,
     327!     .     jjm, rlonu, rlatv , interbar )
    329328
    330329      varname = 'zmea'
    … …  
    511510      ENDIF
    512511
     512!      where (pctsrf(1:klon, is_ter) >= .5)
     513!        pctsrf(1:klon, is_ter) = 1.
     514!        pctsrf(1:klon, is_oce) = 0.
     515!        pctsrf(1:klon, is_sic) = 0.
     516!        pctsrf(1:klon, is_lic) = 0.
     517!        zmasq = 1.
     518!      endwhere
     519!      where (pctsrf(1:klon, is_lic) >= .5)
     520!        pctsrf(1:klon, is_ter) = 0.
     521!        pctsrf(1:klon, is_oce) = 0.
     522!        pctsrf(1:klon, is_sic) = 0.
     523!        pctsrf(1:klon, is_lic) = 1.
     524!        zmasq = 1.
     525!      endwhere
     526!      where (pctsrf(1:klon, is_oce) >= .5)
     527!        pctsrf(1:klon, is_ter) = 0.
     528!        pctsrf(1:klon, is_oce) = 1.
     529!        pctsrf(1:klon, is_sic) = 0.
     530!        pctsrf(1:klon, is_lic) = 0.
     531!        zmasq = 0.
     532!      endwhere
     533!      where (pctsrf(1:klon, is_sic) >= .5)
     534!        pctsrf(1:klon, is_ter) = 0.
     535!        pctsrf(1:klon, is_oce) = 0.
     536!        pctsrf(1:klon, is_sic) = 1.
     537!        pctsrf(1:klon, is_lic) = 0.
     538!        zmasq = 0.
     539!      endwhere
     540!      call gr_fi_dyn(1, klon, iip1, jjp1, zmasq, masque)
     541C
     542C verif que somme des sous surface = 1
     543C
     544!      ji=count( (abs( sum(pctsrf(1 : klon, 1 : nbsrf), dim = 2)) - 1.0 )
     545!     $    .GT. EPSFRA)
     546!      IF (ji .NE. 0) THEN
     547!          WRITE(*,*) 'pb repartition sous maille pour ',ji,' points'
     548!     ENDIF
    513549
    514550
    … …  
    544580      print*,'sortie inidissip'
    545581      itau = 0
     582      itau_dyn = 0
     583      itau_phy = 0
    546584      iday = dayref +itau/day_step
    547585      time = FLOAT(itau-(iday-dayref)*day_step)/day_step
    … …  
    551589       iday = iday + 1
    552590      ENDIF
     591      day_ref = dayref
     592      annee_ref = anneeref
     593
    553594      CALL geopot  ( ip1jmp1, tpot  , pk , pks,  phis  , phi   )
    554595      print*,'sortie geopot'
    … …  
    557598     *                phi,w, pbaru,pbarv,time+iday-dayref   )
    558599       print*,'sortie caldyn0'     
    559       CALL dynredem0("start.nc",dayref,anneeref,phis,nqmx)
     600      CALL dynredem0("start.nc",dayref,phis,nqmx)
    560601      print*,'sortie dynredem0'
    561602      CALL dynredem1("start.nc",0.0,vvent,uvent,tpot,q3d,nqmx,masse ,
    … …  
    569610      co2_ppm   = 330.0
    570611      solaire   = 1370.0
    571 
    572 c      call physdem(lonfi, latfi, phystep,radpas,co2_ppm,
    573 c     .                   solaire,tsol, qsol,
    574 c     .                   sn, radsol, deltat, rugmer,
    575 c     .                   agesno, zmea, zstd, zsig,
    576 c     .                   zgam, zthe, zpic, zval,
    577 c     .                   rugsrel)
    578612
    579613c
    … …  
    614648      frugs(1:klon,is_lic) = MAX(1.0e-05, zstd(1:klon)*zsig(1:klon)/2.0)
    615649      frugs(1:klon,is_sic) = 0.001
    616 
    617       call physdem("startphy.nc",phystep,radpas, co2_ppm, solaire,
     650      fder = 0.0
     651
     652      call phyredem("startphy.nc",phystep,radpas, co2_ppm, solaire,
    618653     $    latfi, lonfi, pctsrf, tsolsrf, tsoil, deltat, qsolsrf, snsrf,
    619      $    albe, evap, rain_fall, snow_fall, solsw, sollw,
     654     $    albe, evap, rain_fall, snow_fall, solsw, sollw, fder,
    620655     $    radsol, frugs,  agesno,
    621656     $    zmea, zstd, zsig, zgam, zthe, zpic, zval, rugsrel,

  • LMDZ.3.3/branches/rel-LF/libf/phylmd/phyredem.F

    r258 r352  
    9999      tab_cntrl(14) = anneeref
    100100      tab_cntrl(13) = day_end
    101       tab_cntrl(14) = anne_ini
    102 c
    103       ierr = NF_REDEF (nid)
     101      tab_cntrl(14) = annee_ref
     102      tab_cntrl(15) = itau_phy
     103c
     104      ierr = NF_REDEF (nid)
     105#ifdef NC_DOUBLE
     106      ierr = NF_DEF_VAR (nid, "controle", NF_DOUBLE, 1, idim1,nvarid)
     107#else
    104108      ierr = NF_DEF_VAR (nid, "controle", NF_FLOAT, 1, idim1,nvarid)
     109#endif
    105110      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 22,
    106111     .                        "Parametres de controle")
    … …  
    113118c
    114119      ierr = NF_REDEF (nid)
     120#ifdef NC_DOUBLE
     121      ierr = NF_DEF_VAR (nid, "longitude", NF_DOUBLE, 1, idim2,nvarid)
     122#else
    115123      ierr = NF_DEF_VAR (nid, "longitude", NF_FLOAT, 1, idim2,nvarid)
     124#endif
    116125      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 32,
    117126     .                        "Longitudes de la grille physique")
    … …  
    124133c
    125134      ierr = NF_REDEF (nid)
     135#ifdef NC_DOUBLE
     136      ierr = NF_DEF_VAR (nid, "latitude", NF_DOUBLE, 1, idim2,nvarid)
     137#else
    126138      ierr = NF_DEF_VAR (nid, "latitude", NF_FLOAT, 1, idim2,nvarid)
     139#endif
    127140      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 31,
    128141     .                        "Latitudes de la grille physique")
    … …  
    137150C
    138151      ierr = NF_REDEF (nid)
     152#ifdef NC_DOUBLE
     153      ierr = NF_DEF_VAR (nid, "masque", NF_DOUBLE, 1, idim2,nvarid)
     154#else
    139155      ierr = NF_DEF_VAR (nid, "masque", NF_FLOAT, 1, idim2,nvarid)
     156#endif
    140157      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 16,
    141158     .                        "masque terre mer")
    … …  
    151168C
    152169      ierr = NF_REDEF (nid)
     170#ifdef NC_DOUBLE
     171      ierr = NF_DEF_VAR (nid, "FTER", NF_DOUBLE, 1, idim2,nvarid)
     172#else
    153173      ierr = NF_DEF_VAR (nid, "FTER", NF_FLOAT, 1, idim2,nvarid)
     174#endif
    154175      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21,
    155176     .                        "fraction de continent")
    … …  
    164185C
    165186      ierr = NF_REDEF (nid)
     187#ifdef NC_DOUBLE
     188      ierr = NF_DEF_VAR (nid, "FLIC", NF_DOUBLE, 1, idim2,nvarid)
     189#else
    166190      ierr = NF_DEF_VAR (nid, "FLIC", NF_FLOAT, 1, idim2,nvarid)
     191#endif
    167192      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 24,
    168193     .                        "fraction glace de terre")
    … …  
    177202C
    178203      ierr = NF_REDEF (nid)
     204#ifdef NC_DOUBLE
     205      ierr = NF_DEF_VAR (nid, "FOCE", NF_DOUBLE, 1, idim2,nvarid)
     206#else
    179207      ierr = NF_DEF_VAR (nid, "FOCE", NF_FLOAT, 1, idim2,nvarid)
     208#endif
    180209      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 14,
    181210     .                        "fraction ocean")
    … …  
    190219C
    191220      ierr = NF_REDEF (nid)
     221#ifdef NC_DOUBLE
     222      ierr = NF_DEF_VAR (nid, "FSIC", NF_DOUBLE, 1, idim2,nvarid)
     223#else
    192224      ierr = NF_DEF_VAR (nid, "FSIC", NF_FLOAT, 1, idim2,nvarid)
     225#endif
    193226      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 18,
    194227     .                        "fraction glace mer")
    … …  
    206239        WRITE(str2,'(i2.2)') nsrf
    207240        ierr = NF_REDEF (nid)
     241#ifdef NC_DOUBLE
     242        ierr = NF_DEF_VAR (nid, "TS"//str2, NF_DOUBLE, 1, idim2,nvarid)
     243#else
    208244        ierr = NF_DEF_VAR (nid, "TS"//str2, NF_FLOAT, 1, idim2,nvarid)
     245#endif
    209246        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
    210247     .                        "Temperature de surface No."//str2)
    … …  
    226263        WRITE(str7,'(i2.2,"srf",i2.2)') isoil,nsrf
    227264        ierr = NF_REDEF (nid)
     265#ifdef NC_DOUBLE
     266        ierr = NF_DEF_VAR (nid, "Tsoil"//str7,NF_DOUBLE,1,idim2,nvarid)
     267#else
    228268        ierr = NF_DEF_VAR (nid, "Tsoil"//str7,NF_FLOAT,1,idim2,nvarid)
     269#endif
    229270        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 29,
    230271     .                        "Temperature du sol No."//str7)
    … …  
    244285c
    245286      ierr = NF_REDEF (nid)
     287#ifdef NC_DOUBLE
     288      ierr = NF_DEF_VAR (nid, "DELTAT", NF_DOUBLE, 1, idim2,nvarid)
     289#else
    246290      ierr = NF_DEF_VAR (nid, "DELTAT", NF_FLOAT, 1, idim2,nvarid)
     291#endif
    247292      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 33,
    248293     .                        "Ecart de la SST (pour slab-ocean)")
    … …  
    258303        WRITE(str2,'(i2.2)') nsrf
    259304        ierr = NF_REDEF (nid)
     305#ifdef NC_DOUBLE
     306        ierr = NF_DEF_VAR (nid,"QS"//str2,NF_DOUBLE,1,idim2,nvarid)
     307#else
    260308        ierr = NF_DEF_VAR (nid,"QS"//str2,NF_FLOAT,1,idim2,nvarid)
     309#endif
    261310        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 25,
    262311     .                        "Humidite de surface No."//str2)
    … …  
    277326        WRITE(str2,'(i2.2)') nsrf
    278327        ierr = NF_REDEF (nid)
     328#ifdef NC_DOUBLE
     329        ierr = NF_DEF_VAR (nid,"ALBE"//str2,NF_DOUBLE,1,idim2,nvarid)
     330#else
    279331        ierr = NF_DEF_VAR (nid,"ALBE"//str2,NF_FLOAT,1,idim2,nvarid)
     332#endif
    280333        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 23,
    281334     .                        "albedo de surface No."//str2)
    … …  
    296349        WRITE(str2,'(i2.2)') nsrf
    297350        ierr = NF_REDEF (nid)
     351#ifdef NC_DOUBLE
     352        ierr = NF_DEF_VAR (nid,"EVAP"//str2,NF_DOUBLE,1,idim2,nvarid)
     353#else
    298354        ierr = NF_DEF_VAR (nid,"EVAP"//str2,NF_FLOAT,1,idim2,nvarid)
     355#endif
    299356        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
    300357     .                        "Evaporation de surface No."//str2)
    … …  
    316373        WRITE(str2,'(i2.2)') nsrf
    317374        ierr = NF_REDEF (nid)
     375#ifdef NC_DOUBLE
     376        ierr = NF_DEF_VAR (nid,"SNOW"//str2,NF_DOUBLE,1,idim2,nvarid)
     377#else
    318378        ierr = NF_DEF_VAR (nid,"SNOW"//str2,NF_FLOAT,1,idim2,nvarid)
     379#endif
    319380        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 22,
    320381     .                        "Neige de surface No."//str2)
    … …  
    333394c
    334395      ierr = NF_REDEF (nid)
     396#ifdef NC_DOUBLE
     397      ierr = NF_DEF_VAR (nid, "RADS", NF_DOUBLE, 1, idim2,nvarid)
     398#else
    335399      ierr = NF_DEF_VAR (nid, "RADS", NF_FLOAT, 1, idim2,nvarid)
     400#endif
    336401      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
    337402     .                        "Rayonnement net a la surface")
    … …  
    344409c
    345410      ierr = NF_REDEF (nid)
     411#ifdef NC_DOUBLE
     412      ierr = NF_DEF_VAR (nid, "solsw", NF_DOUBLE, 1, idim2,nvarid)
     413#else
    346414      ierr = NF_DEF_VAR (nid, "solsw", NF_FLOAT, 1, idim2,nvarid)
     415#endif
    347416      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 32,
    348417     .                        "Rayonnement solaire a la surface")
    … …  
    355424c
    356425      ierr = NF_REDEF (nid)
     426#ifdef NC_DOUBLE
     427      ierr = NF_DEF_VAR (nid, "sollw", NF_DOUBLE, 1, idim2,nvarid)
     428#else
    357429      ierr = NF_DEF_VAR (nid, "sollw", NF_FLOAT, 1, idim2,nvarid)
     430#endif
    358431      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 27,
    359432     .                        "Rayonnement IF a la surface")
    … …  
    366439c
    367440      ierr = NF_REDEF (nid)
     441#ifdef NC_DOUBLE
     442      ierr = NF_DEF_VAR (nid, "fder", NF_DOUBLE, 1, idim2,nvarid)
     443#else
    368444      ierr = NF_DEF_VAR (nid, "fder", NF_FLOAT, 1, idim2,nvarid)
     445#endif
    369446      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 14,
    370447     .                        "Derive de flux")
    … …  
    377454c
    378455      ierr = NF_REDEF (nid)
     456#ifdef NC_DOUBLE
     457      ierr = NF_DEF_VAR (nid, "rain_f", NF_DOUBLE, 1, idim2,nvarid)
     458#else
    379459      ierr = NF_DEF_VAR (nid, "rain_f", NF_FLOAT, 1, idim2,nvarid)
     460#endif
    380461      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21,
    381462     .                        "precipitation liquide")
    … …  
    388469c
    389470      ierr = NF_REDEF (nid)
     471#ifdef NC_DOUBLE
     472      ierr = NF_DEF_VAR (nid, "snow_f", NF_DOUBLE, 1, idim2,nvarid)
     473#else
    390474      ierr = NF_DEF_VAR (nid, "snow_f", NF_FLOAT, 1, idim2,nvarid)
     475#endif
    391476      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20,
    392477     .                        "precipitation solide")
    … …  
    402487        WRITE(str2,'(i2.2)') nsrf
    403488        ierr = NF_REDEF (nid)
     489#ifdef NC_DOUBLE
     490        ierr = NF_DEF_VAR (nid,"RUG"//str2,NF_DOUBLE,1,idim2,nvarid)
     491#else
    404492        ierr = NF_DEF_VAR (nid,"RUG"//str2,NF_FLOAT,1,idim2,nvarid)
     493#endif
    405494        ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 23,
    406495     .                        "rugosite de surface No."//str2)
    … …  
    421510            WRITE(str2,'(i2.2)') nsrf
    422511            ierr = NF_REDEF (nid)
     512#ifdef NC_DOUBLE
     513            ierr = NF_DEF_VAR (nid,"AGESNO"//str2,NF_DOUBLE,1,idim2
     514     $          ,nvarid)
     515#else
    423516            ierr = NF_DEF_VAR (nid,"AGESNO"//str2,NF_FLOAT,1,idim2
    424517     $          ,nvarid)
     518#endif
    425519            ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 15,
    426520     .                        "Age de la neige surface No."//str2)
    … …  
    438532c
    439533      ierr = NF_REDEF (nid)
     534#ifdef NC_DOUBLE
     535      ierr = NF_DEF_VAR (nid, "ZMEA", NF_DOUBLE, 1, idim2,nvarid)
     536#else
    440537      ierr = NF_DEF_VAR (nid, "ZMEA", NF_FLOAT, 1, idim2,nvarid)
     538#endif
    441539      ierr = NF_ENDDEF(nid)
    442540#ifdef NC_DOUBLE
    … …  
    447545c
    448546      ierr = NF_REDEF (nid)
     547#ifdef NC_DOUBLE
     548      ierr = NF_DEF_VAR (nid, "ZSTD", NF_DOUBLE, 1, idim2,nvarid)
     549#else
    449550      ierr = NF_DEF_VAR (nid, "ZSTD", NF_FLOAT, 1, idim2,nvarid)
     551#endif
    450552      ierr = NF_ENDDEF(nid)
    451553#ifdef NC_DOUBLE
    … …  
    454556      ierr = NF_PUT_VAR_REAL (nid,nvarid,zstd)
    455557#endif
    456 c
    457       ierr = NF_REDEF (nid)
     558      ierr = NF_REDEF (nid)
     559#ifdef NC_DOUBLE
     560      ierr = NF_DEF_VAR (nid, "ZSIG", NF_DOUBLE, 1, idim2,nvarid)
     561#else
    458562      ierr = NF_DEF_VAR (nid, "ZSIG", NF_FLOAT, 1, idim2,nvarid)
     563#endif
    459564      ierr = NF_ENDDEF(nid)
    460565#ifdef NC_DOUBLE
    … …  
    463568      ierr = NF_PUT_VAR_REAL (nid,nvarid,zsig)
    464569#endif
    465 c
    466       ierr = NF_REDEF (nid)
     570      ierr = NF_REDEF (nid)
     571#ifdef NC_DOUBLE
     572      ierr = NF_DEF_VAR (nid, "ZGAM", NF_DOUBLE, 1, idim2,nvarid)
     573#else
    467574      ierr = NF_DEF_VAR (nid, "ZGAM", NF_FLOAT, 1, idim2,nvarid)
     575#endif
    468576      ierr = NF_ENDDEF(nid)
    469577#ifdef NC_DOUBLE
    … …  
    472580      ierr = NF_PUT_VAR_REAL (nid,nvarid,zgam)
    473581#endif
    474 c
    475       ierr = NF_REDEF (nid)
     582      ierr = NF_REDEF (nid)
     583#ifdef NC_DOUBLE
     584      ierr = NF_DEF_VAR (nid, "ZTHE", NF_DOUBLE, 1, idim2,nvarid)
     585#else
    476586      ierr = NF_DEF_VAR (nid, "ZTHE", NF_FLOAT, 1, idim2,nvarid)
     587#endif
    477588      ierr = NF_ENDDEF(nid)
    478589#ifdef NC_DOUBLE
    … …  
    481592      ierr = NF_PUT_VAR_REAL (nid,nvarid,zthe)
    482593#endif
    483 c
    484       ierr = NF_REDEF (nid)
     594      ierr = NF_REDEF (nid)
     595#ifdef NC_DOUBLE
     596      ierr = NF_DEF_VAR (nid, "ZPIC", NF_DOUBLE, 1, idim2,nvarid)
     597#else
    485598      ierr = NF_DEF_VAR (nid, "ZPIC", NF_FLOAT, 1, idim2,nvarid)
     599#endif
    486600      ierr = NF_ENDDEF(nid)
    487601#ifdef NC_DOUBLE
    … …  
    490604      ierr = NF_PUT_VAR_REAL (nid,nvarid,zpic)
    491605#endif
    492 c
    493       ierr = NF_REDEF (nid)
     606      ierr = NF_REDEF (nid)
     607#ifdef NC_DOUBLE
     608      ierr = NF_DEF_VAR (nid, "ZVAL", NF_DOUBLE, 1, idim2,nvarid)
     609#else
    494610      ierr = NF_DEF_VAR (nid, "ZVAL", NF_FLOAT, 1, idim2,nvarid)
     611#endif
    495612      ierr = NF_ENDDEF(nid)
    496613#ifdef NC_DOUBLE
    … …  
    499616      ierr = NF_PUT_VAR_REAL (nid,nvarid,zval)
    500617#endif
    501 c
    502       ierr = NF_REDEF (nid)
     618      ierr = NF_REDEF (nid)
     619#ifdef NC_DOUBLE
     620      ierr = NF_DEF_VAR (nid, "RUGSREL", NF_DOUBLE, 1, idim2,nvarid)
     621#else
    503622      ierr = NF_DEF_VAR (nid, "RUGSREL", NF_FLOAT, 1, idim2,nvarid)
     623#endif
    504624      ierr = NF_ENDDEF(nid)
    505625#ifdef NC_DOUBLE
    … …  
    510630c
    511631      ierr = NF_REDEF (nid)
     632#ifdef NC_DOUBLE
     633      ierr = NF_DEF_VAR (nid, "TANCIEN", NF_DOUBLE, 1, idim3,nvarid)
     634#else
    512635      ierr = NF_DEF_VAR (nid, "TANCIEN", NF_FLOAT, 1, idim3,nvarid)
     636#endif
    513637      ierr = NF_ENDDEF(nid)
    514638#ifdef NC_DOUBLE
    … …  
    519643c
    520644      ierr = NF_REDEF (nid)
     645#ifdef NC_DOUBLE
     646      ierr = NF_DEF_VAR (nid, "QANCIEN", NF_DOUBLE, 1, idim3,nvarid)
     647#else
    521648      ierr = NF_DEF_VAR (nid, "QANCIEN", NF_FLOAT, 1, idim3,nvarid)
     649#endif
    522650      ierr = NF_ENDDEF(nid)
    523651#ifdef NC_DOUBLE
    … …  
    528656c
    529657      ierr = NF_REDEF (nid)
     658#ifdef NC_DOUBLE
     659      ierr = NF_DEF_VAR (nid, "RUGMER", NF_DOUBLE, 1, idim2,nvarid)
     660#else
    530661      ierr = NF_DEF_VAR (nid, "RUGMER", NF_FLOAT, 1, idim2,nvarid)
     662#endif
    531663      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 28,
    532664     .                        "Longueur de rugosite sur mer")
Note: See TracChangeset for help on using the changeset viewer.