Changeset 533 for trunk/LMDZ.GENERIC


Ignore:
Timestamp:
Feb 16, 2012, 5:52:18 PM (13 years ago)
Author:
jleconte
Message:
  • Mesh area now present in 1d diag files.
  • Spectral Bandwidth included in the diagspec files
  • Correction of a initialzation bug on day_step and ecritphy in rcm1d
Location:
trunk/LMDZ.GENERIC
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/README

    r527 r533  
    583583- corrects a small bug when no clouds are present in aeropacity
    584584
    585 
    586 
    587 
    588 
    589 
     585== 16/02/2012 == JL + EM
     586- Mesh area now present in 1d diag files.
     587- Spectral Bandwidth included in the diagspec files
     588- Correction of a initialzation bug on day_step and ecritphy in rcm1d
     589
     590
     591
     592
     593
  • trunk/LMDZ.GENERIC/libf/phystd/iniwrite_specIR.F

    r526 r533  
    22
    33      use radinc_h, only: L_NSPECTI
    4       use radcommon_h, only: WNOI
     4      use radcommon_h, only: WNOI,DWNI
    55
    66      implicit none
     
    107107
    108108      ierr = NF_DEF_DIM (nid, "index", length, idim_index)
    109 !      ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
    110109      ierr = NF_DEF_DIM (nid, "latitude", jjp1, idim_rlatu)
    111110      ierr = NF_DEF_DIM (nid, "longitude", iip1, idim_rlonv)
    112 !      ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)
    113 !      ierr = NF_DEF_DIM (nid, "interlayer", (llm+1), idim_llmp1)
    114       ierr = NF_DEF_DIM (nid, "altitude", llm, idim_llm)
    115 !      ierr = NF_DEF_DIM (nid, "subsurface_layers",nsoilmx,idim_nsoilmx)
    116111      ierr = NF_DEF_DIM (nid, "IR Wavenumber",L_NSPECTI,idim_bandsIR)
    117112
     
    140135
    141136
    142 c$$$      ierr = NF_REDEF (nid)
    143 c$$$#ifdef NC_DOUBLE
    144 c$$$      ierr = NF_DEF_VAR (nid, "rlonu", NF_DOUBLE, 1, idim_rlonu,nvarid)
    145 c$$$#else
    146 c$$$      ierr = NF_DEF_VAR (nid, "rlonu", NF_FLOAT, 1, idim_rlonu,nvarid)
    147 c$$$#endif
    148 c$$$      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21,
    149 c$$$     .                       "Longitudes at u nodes")
    150 c$$$      ierr = NF_ENDDEF(nid)
    151 c$$$#ifdef NC_DOUBLE
    152 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonu/pi*180)
    153 c$$$#else
    154 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonu/pi*180)
    155 c$$$#endif
    156 
    157 c --------------------------
    158137      ierr = NF_REDEF (nid)
    159138#ifdef NC_DOUBLE
     
    189168#endif
    190169c
    191 c --------------------------
    192       ierr = NF_REDEF (nid)
    193 #ifdef NC_DOUBLE
    194       ierr = NF_DEF_VAR (nid, "altitude", NF_DOUBLE, 1,
    195      .       idim_llm,nvarid)
    196 #else
    197       ierr = NF_DEF_VAR (nid, "altitude", NF_FLOAT, 1,
    198      .       idim_llm,nvarid)
    199 #endif
    200       ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",10,"pseudo-alt")
    201       ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"km")
    202       ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',2,"up")
    203 
    204       ierr = NF_ENDDEF(nid)
    205 #ifdef NC_DOUBLE
    206       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pseudoalt)
    207 #else
    208       ierr = NF_PUT_VAR_REAL (nid,nvarid,pseudoalt)
    209 #endif
    210 
    211 
    212 c$$$c --------------------------
    213 c$$$      ierr = NF_REDEF (nid)
    214 c$$$#ifdef NC_DOUBLE
    215 c$$$      ierr = NF_DEF_VAR (nid, "rlatv", NF_DOUBLE, 1, idim_rlatv,nvarid)
    216 c$$$#else
    217 c$$$      ierr = NF_DEF_VAR (nid, "rlatv", NF_FLOAT, 1, idim_rlatv,nvarid)
    218 c$$$#endif
    219 c$$$      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20,
    220 c$$$     .                       "Latitudes at v nodes")
    221 c$$$      ierr = NF_ENDDEF(nid)
    222 c$$$#ifdef NC_DOUBLE
    223 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatv/pi*180)
    224 c$$$#else
    225 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatv/pi*180)
    226 c$$$#endif
    227 
    228 
    229 c$$$c --------------------------
    230 c$$$c  Vertical levels
    231 c$$$      call def_var(nid,"aps","hybrid pressure at midlayers ","Pa",
    232 c$$$     .            1,idim_llm,nvarid,ierr)
    233 c$$$#ifdef NC_DOUBLE
    234 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aps)
    235 c$$$#else
    236 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,aps)
    237 c$$$#endif
    238 c$$$
    239 c$$$      call def_var(nid,"bps","hybrid sigma at midlayers"," ",
    240 c$$$     .            1,idim_llm,nvarid,ierr)
    241 c$$$#ifdef NC_DOUBLE
    242 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bps)
    243 c$$$#else
    244 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,bps)
    245 c$$$#endif
    246 c$$$
    247 c$$$      call def_var(nid,"ap","hybrid pressure at interlayers","Pa",
    248 c$$$     .            1,idim_llmp1,nvarid,ierr)
    249 c$$$#ifdef NC_DOUBLE
    250 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ap)
    251 c$$$#else
    252 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,ap)
    253 c$$$#endif
    254 c$$$
    255 c$$$      call def_var(nid,"bp","hybrid sigma at interlayers"," ",
    256 c$$$     .            1,idim_llmp1,nvarid,ierr)
    257 c$$$#ifdef NC_DOUBLE
    258 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bp)
    259 c$$$#else
    260 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,bp)
    261 c$$$#endif
    262 
    263 c$$$!-------------------------------
    264 c$$$! (soil) depth variable mlayer() (known from comsoil.h)
    265 c$$$!-------------------------------
    266 c$$$      ierr=NF_REDEF (nid) ! Enter NetCDF (re-)define mode
    267 c$$$      ! define variable
    268 c$$$#ifdef NC_DOUBLE
    269 c$$$      ierr=NF_DEF_VAR(nid,"soildepth",NF_DOUBLE,1,idim_nsoilmx,nvarid)
    270 c$$$#else
    271 c$$$      ierr=NF_DEF_VAR(nid,"soildepth",NF_FLOAT,1,idim_nsoilmx,nvarid)
    272 c$$$#endif
    273 c$$$      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 20,
    274 c$$$     .                        "Soil mid-layer depth")
    275 c$$$      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"units",1,"m")
    276 c$$$      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"positive",4,"down")
    277 c$$$      ierr=NF_ENDDEF(nid) ! Leave NetCDF define mode
    278 c$$$      ! write variable
    279 c$$$#ifdef NC_DOUBLE
    280 c$$$      ierr=NF_PUT_VAR_DOUBLE (nid,nvarid,mlayer)
    281 c$$$#else
    282 c$$$      ierr=NF_PUT_VAR_REAL (nid,nvarid,mlayer)
    283 c$$$#endif
    284170
    285171!-------------------------------
     
    307193#endif
    308194
     195!-------------------------------
     196! Width of bands in the IR
     197!-------------------------------
     198
     199      ierr=NF_REDEF (nid) ! Enter NetCDF (re-)define mode
     200      ! define variable
     201#ifdef NC_DOUBLE
     202      ierr=NF_DEF_VAR(nid,"IR Bandwidth",NF_DOUBLE,1,
     203     .                          idim_bandsIR,nvarid)
     204#else
     205      ierr=NF_DEF_VAR(nid,"IR Bandwidth",NF_FLOAT,1,
     206     .                          idim_bandsIR,nvarid)
     207#endif
     208      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 25,
     209     .                        "Bandwidth in the infrared")
     210      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"units",5,"cm^-1")
     211      ierr=NF_ENDDEF(nid) ! Leave NetCDF define mode
     212      ! write variable
     213#ifdef NC_DOUBLE
     214      ierr=NF_PUT_VAR_DOUBLE (nid,nvarid,dble(DWNI))
     215#else
     216      ierr=NF_PUT_VAR_REAL (nid,nvarid,real(DWNI))
     217#endif
    309218c
    310219c --------------------------
    311220c  Mesh area and conversion coefficients cov. <-> contra. <--> natural
    312221
    313       id(1)=idim_rlonu
    314       id(2)=idim_rlatu
    315 c
    316       ierr = NF_REDEF (nid)
    317 #ifdef NC_DOUBLE
    318       ierr = NF_DEF_VAR (nid, "cu", NF_DOUBLE, 2, id,nvarid)
    319 #else
    320       ierr = NF_DEF_VAR (nid, "cu", NF_FLOAT, 2, id,nvarid)
    321 #endif
    322       ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 40,
    323      .             "Conversion coefficients cov <--> natural")
    324       ierr = NF_ENDDEF(nid)
    325 #ifdef NC_DOUBLE
    326       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cu)
    327 #else
    328       ierr = NF_PUT_VAR_REAL (nid,nvarid,cu)
    329 #endif
    330 c
    331       id(1)=idim_rlonv
    332       id(2)=idim_rlatv
    333 c
    334 c --------------------------
    335       ierr = NF_REDEF (nid)
    336 #ifdef NC_DOUBLE
    337       ierr = NF_DEF_VAR (nid, "cv", NF_DOUBLE, 2, id,nvarid)
    338 #else
    339       ierr = NF_DEF_VAR (nid, "cv", NF_FLOAT, 2, id,nvarid)
    340 #endif
    341       ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 40,
    342      .             "Conversion coefficients cov <--> natural")
    343       ierr = NF_ENDDEF(nid)
    344 #ifdef NC_DOUBLE
    345       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cv)
    346 #else
    347       ierr = NF_PUT_VAR_REAL (nid,nvarid,cv)
    348 #endif
    349 c
     222
    350223      id(1)=idim_rlonv
    351224      id(2)=idim_rlatu
     
    367240#endif
    368241
    369 c  Surface geopotential
    370       id(1)=idim_rlonv
    371       id(2)=idim_rlatu
    372 
    373       ierr = NF_REDEF (nid)
    374 #ifdef NC_DOUBLE
    375       ierr = NF_DEF_VAR (nid, "phisinit", NF_DOUBLE, 2, id,nvarid)
    376 #else
    377       ierr = NF_DEF_VAR (nid, "phisinit", NF_FLOAT, 2, id,nvarid)
    378 #endif
    379       ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 27,
    380      .                       "Geopotential at the surface")
    381       ierr = NF_ENDDEF(nid)
    382 #ifdef NC_DOUBLE
    383       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,phis)
    384 #else
    385       ierr = NF_PUT_VAR_REAL (nid,nvarid,phis)
    386 #endif,
    387 
    388 
    389 !      write(*,*)'iniwrite: iim,jjm,llm,idayref',iim,jjm,llm,idayref
    390 !      write(*,*)'iniwrite: rad,omeg,g,mugaz,rcp',rad,omeg,g,mugaz,rcp
    391 !      write(*,*)'iniwrite: daysec,dtphys',daysec,dtphys
    392242
    393243      END
  • trunk/LMDZ.GENERIC/libf/phystd/iniwrite_specVI.F

    r526 r533  
    22
    33      use radinc_h, only: L_NSPECTV
    4       use radcommon_h, only: WNOV
     4      use radcommon_h, only: WNOV,DWNV
    55
    66      implicit none
     
    107107
    108108      ierr = NF_DEF_DIM (nid, "index", length, idim_index)
    109 !      ierr = NF_DEF_DIM (nid, "rlonu", iip1, idim_rlonu)
    110109      ierr = NF_DEF_DIM (nid, "latitude", jjp1, idim_rlatu)
    111110      ierr = NF_DEF_DIM (nid, "longitude", iip1, idim_rlonv)
    112 !      ierr = NF_DEF_DIM (nid, "rlatv", jjm, idim_rlatv)
    113 !      ierr = NF_DEF_DIM (nid, "interlayer", (llm+1), idim_llmp1)
    114       ierr = NF_DEF_DIM (nid, "altitude", llm, idim_llm)
    115 !      ierr = NF_DEF_DIM (nid, "subsurface_layers",nsoilmx,idim_nsoilmx)
    116111      ierr = NF_DEF_DIM (nid, "VI Wavenumber",L_NSPECTV,idim_bandsVI)
    117112
     
    139134c  longitudes and latitudes
    140135
    141 
    142 c$$$      ierr = NF_REDEF (nid)
    143 c$$$#ifdef NC_DOUBLE
    144 c$$$      ierr = NF_DEF_VAR (nid, "rlonu", NF_DOUBLE, 1, idim_rlonu,nvarid)
    145 c$$$#else
    146 c$$$      ierr = NF_DEF_VAR (nid, "rlonu", NF_FLOAT, 1, idim_rlonu,nvarid)
    147 c$$$#endif
    148 c$$$      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 21,
    149 c$$$     .                       "Longitudes at u nodes")
    150 c$$$      ierr = NF_ENDDEF(nid)
    151 c$$$#ifdef NC_DOUBLE
    152 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlonu/pi*180)
    153 c$$$#else
    154 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlonu/pi*180)
    155 c$$$#endif
    156 
    157 c --------------------------
    158136      ierr = NF_REDEF (nid)
    159137#ifdef NC_DOUBLE
     
    189167#endif
    190168c
    191 c --------------------------
    192       ierr = NF_REDEF (nid)
    193 #ifdef NC_DOUBLE
    194       ierr = NF_DEF_VAR (nid, "altitude", NF_DOUBLE, 1,
    195      .       idim_llm,nvarid)
    196 #else
    197       ierr = NF_DEF_VAR (nid, "altitude", NF_FLOAT, 1,
    198      .       idim_llm,nvarid)
    199 #endif
    200       ierr = NF_PUT_ATT_TEXT (nid,nvarid,"long_name",10,"pseudo-alt")
    201       ierr = NF_PUT_ATT_TEXT (nid,nvarid,'units',2,"km")
    202       ierr = NF_PUT_ATT_TEXT (nid,nvarid,'positive',2,"up")
    203 
    204       ierr = NF_ENDDEF(nid)
    205 #ifdef NC_DOUBLE
    206       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,pseudoalt)
    207 #else
    208       ierr = NF_PUT_VAR_REAL (nid,nvarid,pseudoalt)
    209 #endif
    210 
    211 
    212 c$$$c --------------------------
    213 c$$$      ierr = NF_REDEF (nid)
    214 c$$$#ifdef NC_DOUBLE
    215 c$$$      ierr = NF_DEF_VAR (nid, "rlatv", NF_DOUBLE, 1, idim_rlatv,nvarid)
    216 c$$$#else
    217 c$$$      ierr = NF_DEF_VAR (nid, "rlatv", NF_FLOAT, 1, idim_rlatv,nvarid)
    218 c$$$#endif
    219 c$$$      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 20,
    220 c$$$     .                       "Latitudes at v nodes")
    221 c$$$      ierr = NF_ENDDEF(nid)
    222 c$$$#ifdef NC_DOUBLE
    223 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,rlatv/pi*180)
    224 c$$$#else
    225 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,rlatv/pi*180)
    226 c$$$#endif
    227 
    228 
    229 c$$$c --------------------------
    230 c$$$c  Vertical levels
    231 c$$$      call def_var(nid,"aps","hybrid pressure at midlayers ","Pa",
    232 c$$$     .            1,idim_llm,nvarid,ierr)
    233 c$$$#ifdef NC_DOUBLE
    234 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,aps)
    235 c$$$#else
    236 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,aps)
    237 c$$$#endif
    238 c$$$
    239 c$$$      call def_var(nid,"bps","hybrid sigma at midlayers"," ",
    240 c$$$     .            1,idim_llm,nvarid,ierr)
    241 c$$$#ifdef NC_DOUBLE
    242 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bps)
    243 c$$$#else
    244 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,bps)
    245 c$$$#endif
    246 c$$$
    247 c$$$      call def_var(nid,"ap","hybrid pressure at interlayers","Pa",
    248 c$$$     .            1,idim_llmp1,nvarid,ierr)
    249 c$$$#ifdef NC_DOUBLE
    250 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,ap)
    251 c$$$#else
    252 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,ap)
    253 c$$$#endif
    254 c$$$
    255 c$$$      call def_var(nid,"bp","hybrid sigma at interlayers"," ",
    256 c$$$     .            1,idim_llmp1,nvarid,ierr)
    257 c$$$#ifdef NC_DOUBLE
    258 c$$$      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,bp)
    259 c$$$#else
    260 c$$$      ierr = NF_PUT_VAR_REAL (nid,nvarid,bp)
    261 c$$$#endif
    262 
    263 c$$$!-------------------------------
    264 c$$$! (soil) depth variable mlayer() (known from comsoil.h)
    265 c$$$!-------------------------------
    266 c$$$      ierr=NF_REDEF (nid) ! Enter NetCDF (re-)define mode
    267 c$$$      ! define variable
    268 c$$$#ifdef NC_DOUBLE
    269 c$$$      ierr=NF_DEF_VAR(nid,"soildepth",NF_DOUBLE,1,idim_nsoilmx,nvarid)
    270 c$$$#else
    271 c$$$      ierr=NF_DEF_VAR(nid,"soildepth",NF_FLOAT,1,idim_nsoilmx,nvarid)
    272 c$$$#endif
    273 c$$$      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 20,
    274 c$$$     .                        "Soil mid-layer depth")
    275 c$$$      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"units",1,"m")
    276 c$$$      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"positive",4,"down")
    277 c$$$      ierr=NF_ENDDEF(nid) ! Leave NetCDF define mode
    278 c$$$      ! write variable
    279 c$$$#ifdef NC_DOUBLE
    280 c$$$      ierr=NF_PUT_VAR_DOUBLE (nid,nvarid,mlayer)
    281 c$$$#else
    282 c$$$      ierr=NF_PUT_VAR_REAL (nid,nvarid,mlayer)
    283 c$$$#endif
    284169
    285170!-------------------------------
     
    308193
    309194c
    310 c --------------------------
    311 c  Mesh area and conversion coefficients cov. <-> contra. <--> natural
    312 
    313       id(1)=idim_rlonu
    314       id(2)=idim_rlatu
    315 c
    316       ierr = NF_REDEF (nid)
    317 #ifdef NC_DOUBLE
    318       ierr = NF_DEF_VAR (nid, "cu", NF_DOUBLE, 2, id,nvarid)
    319 #else
    320       ierr = NF_DEF_VAR (nid, "cu", NF_FLOAT, 2, id,nvarid)
    321 #endif
    322       ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 40,
    323      .             "Conversion coefficients cov <--> natural")
    324       ierr = NF_ENDDEF(nid)
    325 #ifdef NC_DOUBLE
    326       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cu)
    327 #else
    328       ierr = NF_PUT_VAR_REAL (nid,nvarid,cu)
    329 #endif
    330 c
    331       id(1)=idim_rlonv
    332       id(2)=idim_rlatv
    333 c
    334 c --------------------------
    335       ierr = NF_REDEF (nid)
    336 #ifdef NC_DOUBLE
    337       ierr = NF_DEF_VAR (nid, "cv", NF_DOUBLE, 2, id,nvarid)
    338 #else
    339       ierr = NF_DEF_VAR (nid, "cv", NF_FLOAT, 2, id,nvarid)
    340 #endif
    341       ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 40,
    342      .             "Conversion coefficients cov <--> natural")
    343       ierr = NF_ENDDEF(nid)
    344 #ifdef NC_DOUBLE
    345       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,cv)
    346 #else
    347       ierr = NF_PUT_VAR_REAL (nid,nvarid,cv)
    348 #endif
    349 c
     195!-------------------------------
     196! Width of bands in the Visible
     197!-------------------------------
     198
     199      ierr=NF_REDEF (nid) ! Enter NetCDF (re-)define mode
     200      ! define variable
     201#ifdef NC_DOUBLE
     202      ierr=NF_DEF_VAR(nid,"VI Bandwidth",NF_DOUBLE,1,
     203     .                          idim_bandsVI,nvarid)
     204#else
     205      ierr=NF_DEF_VAR(nid,"VI Bandwidth",NF_FLOAT,1,
     206     .                          idim_bandsVI,nvarid)
     207#endif
     208      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 24,
     209     .                        "Bandwidth in the visible")
     210      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"units",5,"cm^-1")
     211      ierr=NF_ENDDEF(nid) ! Leave NetCDF define mode
     212      ! write variable
     213#ifdef NC_DOUBLE
     214      ierr=NF_PUT_VAR_DOUBLE (nid,nvarid,dble(DWNV))
     215#else
     216      ierr=NF_PUT_VAR_REAL (nid,nvarid,real(DWNV))
     217#endif
     218c
     219c --------------------------
     220c   Mesh area
     221c --------------------------
     222
    350223      id(1)=idim_rlonv
    351224      id(2)=idim_rlatu
    352 c
    353 c --------------------------
     225
     226
    354227      ierr = NF_REDEF (nid)
    355228#ifdef NC_DOUBLE
     
    367240#endif
    368241
    369 c  Surface geopotential
    370       id(1)=idim_rlonv
    371       id(2)=idim_rlatu
    372 
    373       ierr = NF_REDEF (nid)
    374 #ifdef NC_DOUBLE
    375       ierr = NF_DEF_VAR (nid, "phisinit", NF_DOUBLE, 2, id,nvarid)
    376 #else
    377       ierr = NF_DEF_VAR (nid, "phisinit", NF_FLOAT, 2, id,nvarid)
    378 #endif
    379       ierr = NF_PUT_ATT_TEXT (nid, nvarid, "title", 27,
    380      .                       "Geopotential at the surface")
    381       ierr = NF_ENDDEF(nid)
    382 #ifdef NC_DOUBLE
    383       ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,phis)
    384 #else
    385       ierr = NF_PUT_VAR_REAL (nid,nvarid,phis)
    386 #endif,
    387 
    388 
    389 !      write(*,*)'iniwrite: iim,jjm,llm,idayref',iim,jjm,llm,idayref
    390 !      write(*,*)'iniwrite: rad,omeg,g,mugaz,rcp',rad,omeg,g,mugaz,rcp
    391 !      write(*,*)'iniwrite: daysec,dtphys',daysec,dtphys
     242
    392243
    393244      END
  • trunk/LMDZ.GENERIC/libf/phystd/writediagspecIR.F

    r526 r533  
    5757#include "temps.h"
    5858#include "surfdat.h"
     59#include "callkeys.h"
    5960
    6061! Arguments on input:
     
    102103!Sortie des variables au rythme voulu
    103104
    104       irythme = int(ecritphy) ! sortie au rythme de ecritphy
    105 !     irythme = iconser  ! sortie au rythme des variables de controle
     105      irythme = int(ecritphy)*iradia ! sortie au rythme de ecritphy*iradia
     106!EM+JL if the spetra need to be output more frequently, need to define a ecritSpec...
    106107!     irythme = iphysiq  ! sortie a tous les pas physique
    107 !     irythme = iecri*day_step ! sortie au rythme des fichiers histoires
    108 !     irythme = periodav*day_step ! sortie au rythme des fichiers histmoy
     108
    109109
    110110!***************************************************************
  • trunk/LMDZ.GENERIC/libf/phystd/writediagspecVI.F

    r526 r533  
    5757#include "temps.h"
    5858#include "surfdat.h"
     59#include "callkeys.h"
    5960
    6061! Arguments on input:
     
    102103!Sortie des variables au rythme voulu
    103104
    104       irythme = int(ecritphy) ! sortie au rythme de ecritphy
    105 !     irythme = iconser  ! sortie au rythme des variables de controle
     105      irythme = int(ecritphy)*iradia ! sortie au rythme de ecritphy
     106!EM+JL if the spetra need to be output more frequently, need to define a ecritSpec...
    106107!     irythme = iphysiq  ! sortie a tous les pas physique
    107 !     irythme = iecri*day_step ! sortie au rythme des fichiers histoires
    108 !     irythme = periodav*day_step ! sortie au rythme des fichiers histmoy
    109108
    110109!***************************************************************
Note: See TracChangeset for help on using the changeset viewer.