Changeset 3725
- Timestamp:
- Apr 15, 2025, 9:34:31 AM (8 weeks ago)
- Location:
- trunk/LMDZ.GENERIC
- Files:
-
- 24 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.GENERIC/changelog.txt
r3723 r3725 2063 2063 == 14/04/2025 == EM+MT 2064 2064 Change default CO2 ice aerosol properties file from old "optprop_co2_vis_ff.dat" 2065 to "optprop_co2ice_vis_n50.dat". Kee possibility for now to still used old file2065 to "optprop_co2ice_vis_n50.dat". Keep possibility for now to still used old file 2066 2066 by seting "forgetit=.true." in suaer_corrk.F90 2067 2067 2068 == 15/04/2025 == EM 2069 Change the way the rate of outputs for diagfi.nc files is specified: 2070 IMPORTANT: Specifying "ecritphy" no longer possible and will trigger an error. 2071 Use "diagfi_output_rate" to specify output rate (in physics steps) instead. 2072 This should makes things (hopefully) clearer for users and also better 2073 enforces a cleaner and clearer separation between dynamics and physics. -
trunk/LMDZ.GENERIC/deftank/callphys.GJ581d
r2972 r3725 3 3 # Run with or without tracer transport ? 4 4 tracer = .true. 5 # rate of outputs for diagfi files (in physics steps, default: 24) 6 diagfi_output_rate = 24 5 7 # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating 6 8 diurnal = .false. -
trunk/LMDZ.GENERIC/deftank/callphys.earlymars
r2972 r3725 3 3 # Run with or without tracer transport ? 4 4 tracer = .true. 5 # rate of outputs for diagfi files (in physics steps, default: 24) 6 diagfi_output_rate = 6 5 7 # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating 6 8 diurnal = .true. -
trunk/LMDZ.GENERIC/deftank/callphys.earthslab
r3100 r3725 3 3 # Run with or without tracer transport ? 4 4 tracer = .true. 5 # rate of outputs for diagfi files (in physics steps, default: 24) 6 diagfi_output_rate = 6 5 7 # Diurnal cycle ? if diurnal=false, diurnally averaged solar heating 6 8 diurnal = .true. -
trunk/LMDZ.GENERIC/deftank/earth1d.rcm1d.def
r3220 r3725 35 35 ndt = 1300 36 36 # Number of steps between each writting in diagfi 37 ecritphy=12037 diagfi_output_rate=120 38 38 39 39 -
trunk/LMDZ.GENERIC/deftank/jupiter1d/rcm1d.def
r1678 r3725 38 38 ndt = 63000 39 39 # ecriture physiq 40 ecritphy= 20040 diagfi_output_rate = 200 41 41 42 42 ## Vertical levels -
trunk/LMDZ.GENERIC/deftank/rcm1d.def.earth
r2110 r3725 15 15 ndt =20 16 16 # Number of steps between each writing in diagfi 17 ecritphy=1217 diagfi_output_rate=12 18 18 19 19 #### Physical parameters -
trunk/LMDZ.GENERIC/deftank/run.def.32x24x25
r2552 r3725 70 70 grireg=.true. 71 71 72 # Output in diagfi file every ecritphy dynamical steps73 ecritphy=12074 75 72 # longitude (degrees) of zoom center 76 73 clon=63. -
trunk/LMDZ.GENERIC/deftank/run.def.64x48x25
r2552 r3725 70 70 grireg=.true. 71 71 72 # Output in diagfi file every ecritphy dynamical steps73 ecritphy=24074 75 72 # longitude (degrees) of zoom center 76 73 clon=63. -
trunk/LMDZ.GENERIC/deftank/run.def.earlymars_hires
r2552 r3725 70 70 grireg=.true. 71 71 72 # Output in diagfi file every ecritphy dynamical steps73 ecritphy=192074 75 72 # longitude (degrees) of zoom center 76 73 clon=63. -
trunk/LMDZ.GENERIC/deftank/run.def.earth_hires
r2552 r3725 67 67 grireg=.true. 68 68 69 # Output in diagfi file every ecritphy dynamical steps70 ecritphy=36071 72 69 # longitude (degrees) of zoom center 73 70 clon=0. -
trunk/LMDZ.GENERIC/deftank/run.def.earthslab.64x48x30
r3100 r3725 65 65 grireg=.true. 66 66 67 # Output in diagfi file every ecritphy dynamical steps68 ecritphy=960069 # ecritphy=170 71 67 # longitude (degrees) of zoom center 72 68 clon=0. -
trunk/LMDZ.GENERIC/deftank/saturn1d/rcm1d.def
r1876 r3725 42 42 ndt = 342020 43 43 # Number of steps between each writing in diagfi 44 ecritphy=500044 diagfi_output_rate=5000 45 45 46 46 # OVERRIDE callphys.def settings -
trunk/LMDZ.GENERIC/libf/aeronostd/calchim_asis.F90
r3312 r3725 10 10 USE comcstfi_mod 11 11 use callkeys_mod 12 use time_phylmdz_mod, only: ecritphy, iphysiq ! output rate set by ecritphy12 use time_phylmdz_mod, only: diagfi_output_rate ! output rate 13 13 use types_asis, only: v_phot_3d, jlabel, nb_phot_hv_max 14 14 use chimiedata_h … … 133 133 integer :: iter(nlayer) ! Number of chemical iterations 134 134 ! within one physical timestep 135 integer :: icount135 integer,intent(in) :: icount ! physics iteration number 136 136 ! for output: 137 137 … … 285 285 ig, zzlay, zzlev, zdens, zycol, ptimestep) 286 286 surface_flux2(ig,:) = surface_flux2(ig,:) + surface_flux(ig,:) 287 if (ngrid==1) then 288 if(mod(icount,ecritphy).eq.0) then 289 surface_flux2(ig,:) = surface_flux2(ig,:)/float(ecritphy) 290 endif 291 else 292 if(mod(icount*iphysiq,ecritphy).eq.0) then 293 surface_flux2(ig,:) = surface_flux2(ig,:)*iphysiq/float(ecritphy) 294 endif 287 if(mod(icount,diagfi_output_rate).eq.0) then 288 !note that surface_flux2(:,:) is reset every diagfi_output_rate 289 surface_flux2(ig,:) = surface_flux2(ig,:)/float(diagfi_output_rate) 295 290 endif 296 end if 291 end if ! of if (depos) 297 292 298 293 !======================================================================= … … 405 400 end do 406 401 ! Restart flux average 407 if (ngrid == 1) then 408 if(mod(icount,ecritphy).eq.0) then 409 surface_flux2(:,:) = 0.0 410 endif 411 else 412 if(mod(icount*iphysiq,ecritphy).eq.0) then 413 surface_flux2(:,:) = 0.0 414 endif 402 if(mod(icount,diagfi_output_rate).eq.0) then 403 surface_flux2(:,:) = 0.0 415 404 endif 416 405 endif ! end depos … … 424 413 425 414 end if ! of if (output) 426 427 return428 415 429 416 -
trunk/LMDZ.GENERIC/libf/aeronostd/writediagspecUV.F
r3309 r3725 30 30 & nbp_lon, nbp_lat, grid_type, 31 31 & unstructured 32 use time_phylmdz_mod, only: ecritphy, iphysiq, day_step, day_ini 32 use time_phylmdz_mod, only: diagfi_output_rate, dtphys, daysec 33 use time_phylmdz_mod, only: day_ini 33 34 34 35 implicit none … … 48 49 ! real dx0 49 50 50 integer i rythme51 integer isample 51 52 integer ierr 52 53 integer iq … … 87 88 !output frequency 88 89 89 i rythme = ecritphy! output with ecritphy frequency90 !EM+JL if the spe tra need to be output more frequently, need to define a ecritSpec...91 ! irythme = iphysiq ! sortie a tous les pas physique90 isample = diagfi_output_rate ! output with ecritphy frequency 91 !EM+JL if the spectra need to be output more frequently, need to define a 92 ! diagspec_output_rate... 92 93 93 94 !*************************************************************** … … 173 174 !------------------------------------------------------------------------ 174 175 if (nom.eq.firstnom) then 175 zitau = zitau + iphysiq176 zitau = zitau + 1 176 177 end if 177 178 … … 180 181 !-------------------------------------------------------- 181 182 182 if ( MOD(zitau+1,i rythme) .eq.0.) then183 if ( MOD(zitau+1,isample) .eq.0.) then 183 184 184 185 ! Compute/write/extend 'Time' coordinate (date given in days) … … 193 194 ntime=ntime+1 ! increment # of stored time steps 194 195 ! compute corresponding date (in days and fractions thereof) 195 date= float (zitau +1)/float (day_step)196 date=(zitau +1)*(dtphys/daysec) 196 197 197 198 if (is_master) then … … 313 314 endif ! of if (dimpx.eq.3) 314 315 315 endif ! of if ( MOD(zitau+1,i rythme) .eq.0.)316 endif ! of if ( MOD(zitau+1,isample) .eq.0.) 316 317 317 318 ! Close the NetCDF file -
trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F
r3716 r3725 21 21 use comcstfi_mod, only: pi, cpp, rad, g, r, 22 22 & mugaz, rcp, omeg 23 use time_phylmdz_mod, only: daysec, dtphys, day_step, ecritphy,24 & nday, iphysiq23 use time_phylmdz_mod, only: daysec, dtphys, nday, 24 & diagfi_output_rate 25 25 use callkeys_mod, only: tracer,rings_shadow, 26 26 & specOLR,water,pceil,ok_slab_ocean,photochem … … 71 71 !================================================================== 72 72 73 #include "dimensions.h" 74 #include "paramet.h" 75 !include "dimphys.h" 76 #include "netcdf.inc" 77 #include "comgeom.h" 73 include "dimensions.h" 74 include "paramet.h" 75 include "comgeom.h" 78 76 79 77 c -------------------------------------------------------------- … … 91 89 c 92 90 INTEGER day0 ! date initial (sol ; =0 a Ls=0) 93 REAL day ! date durant le run 91 REAL day ! date during the run 92 INTEGER day_step ! number of time steps per day 94 93 REAL time ! time (0<time<1 ; time=0.5 a midi) 95 94 REAL play(llm) ! Pressure at the middle of the layers (Pa) … … 596 595 write(*,*) " day_step = ",day_step 597 596 598 iphysiq=1 ! in 1D model physics are called evry time step 599 ecritphy=day_step ! default value for ecritphy 597 diagfi_output_rate=day_step ! default value for diagfi_output_rate 600 598 PRINT *,'Nunber of steps between writediagfi ?' 601 call getin(" ecritphy",ecritphy)602 write(*,*) " ecritphy = ",ecritphy599 call getin("diagfi_output_rate",diagfi_output_rate) 600 write(*,*) " diagfi_output_rate = ",diagfi_output_rate 603 601 604 602 ndt=10 ! default value for ndt -
trunk/LMDZ.GENERIC/libf/phystd/inifis_mod.F90
r3701 r3725 17 17 use comgeomfi_h, only: totarea, totarea_planet 18 18 use comsoil_h, only: ini_comsoil_h, nsoilmx, lay1_soil, alpha_soil 19 use time_phylmdz_mod, only: ecritphy,day_step,iphysiq, &19 use time_phylmdz_mod, only: diagfi_output_rate, & 20 20 init_time, daysec, dtphys 21 21 use comcstfi_mod, only: rad, cpp, g, r, rcp, & … … 76 76 integer,intent(in) :: day_ini 77 77 78 INTEGER ig 78 INTEGER :: ig 79 REAL :: ecritphy ! to check that this obsolete flag is no longer used... 79 80 CHARACTER(len=20) :: rname="inifis" ! routine name, for messages 80 81 … … 103 104 ! read in some parameters from "run.def" for physics, 104 105 ! or shared between dynamics and physics. 106 ecritphy=-666 ! dummy default value 105 107 call getin_p("ecritphy",ecritphy) ! frequency of outputs in physics, 106 108 ! in dynamical steps 107 call getin_p("day_step",day_step) ! number of dynamical steps per day 108 call getin_p("iphysiq",iphysiq) ! call physics every iphysiq dyn step 109 if (ecritphy/=-666) then 110 call abort_physic(rname, & 111 "Error: parameter ecritphy is obsolete! Remove it! "//& 112 "And use diagfi_output_rate instead",1) 113 endif 109 114 110 115 ! do we read a startphy.nc file? (default: .true.) … … 142 147 call getin_p("tracer",tracer) 143 148 if (is_master) write(*,*) trim(rname)//": tracer = ",tracer 149 150 if (is_master) write(*,*) trim(rname)//& 151 ": Output rate for diagfi.nc file (in physics steps) ?" 152 diagfi_output_rate=24 !default value 153 call getin_p("diagfi_output_rate",diagfi_output_rate) 154 if (is_master) write(*,*) trim(rname)//": diagfi_output_rate = ",& 155 diagfi_output_rate 144 156 145 157 if (is_master) write(*,*) trim(rname)//& -
trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90
r3724 r3725 37 37 igcm_h2o_ice, igcm_h2o_vap, igcm_dustbin, & 38 38 igcm_co2_ice, nesp, is_chim, is_condensable,constants_epsi_generic 39 ! igcm_co2_ice, nesp, is_chim, is_condensable,constants_epsi_generic, & 40 ! half_life, top_prod, bot_prod 41 use time_phylmdz_mod, only: ecritphy, iphysiq, nday 39 use time_phylmdz_mod, only: diagfi_output_rate, nday 42 40 use phyetat0_mod, only: phyetat0,tab_cntrl_mod 43 41 use surfini_mod, only: surfini … … 2249 2247 2250 2248 if(meanOLR .and. is_master)then 2251 if((ngrid.gt.1) .or. (mod(icount-1, ecritphy).eq.0))then2249 if((ngrid.gt.1) .or. (mod(icount-1,diagfi_output_rate).eq.0))then 2252 2250 ! to record global radiative balance 2253 2251 open(92,file="rad_bal.out",form='formatted',position='append') … … 2326 2324 2327 2325 if(meanOLR .and. is_master)then 2328 if((ngrid.gt.1) .or. (mod(icount-1, ecritphy).eq.0))then2326 if((ngrid.gt.1) .or. (mod(icount-1,diagfi_output_rate).eq.0))then 2329 2327 ! to record global water balance 2330 2328 open(98,file="h2o_bal.out",form='formatted',position='append') … … 2403 2401 2404 2402 if(lastcall) then 2405 ztime_fin = ptime + ptimestep/(float(iphysiq)*daysec)2403 ztime_fin = ptime + (ptimestep/daysec) 2406 2404 2407 2405 #ifndef MESOSCALE … … 2525 2523 ! OUTPUT in netcdf file "DIAGFI.NC", containing any variable for diagnostic 2526 2524 ! 2527 ! Note 1 : output with period " ecritphy", set in "run.def"2525 ! Note 1 : output with period "diagfi_output_rate", set in "run.def" 2528 2526 ! 2529 2527 ! Note 2 : writediagfi can also be called from any other subroutine for any variable, -
trunk/LMDZ.GENERIC/libf/phystd/time_phylmdz_mod.F90
r1525 r3725 4 4 REAL,SAVE :: dtphys ! physics time step (s) 5 5 !$OMP THREADPRIVATE(dtphys) 6 INTEGER,SAVE :: day_step ! number of dynamical steps per day7 ! (set via inifis)8 !$OMP THREADPRIVATE(day_step)9 6 INTEGER,SAVE :: nday ! number of days to run 10 7 !$OMP THREADPRIVATE(nday) … … 14 11 !$OMP THREADPRIVATE(day_ini) 15 12 16 INTEGER,SAVE :: ecritphy ! for diagfi.nc outputs, write every ecritphy 17 ! dynamical steps (set via inifis) 18 !$OMP THREADPRIVATE(ecritphy) 19 INTEGER,SAVE :: iphysiq ! call physics every iphysiq dynamical step 20 ! (set via inifis) 21 !$OMP THREADPRIVATE(iphysiq) 13 INTEGER,SAVE :: diagfi_output_rate ! for diagfi.nc outputs, write every 14 ! diagfi_output_rate physical step 15 ! (set via inifis) 16 !$OMP THREADPRIVATE(diagfi_output_rate) 22 17 23 18 CONTAINS -
trunk/LMDZ.GENERIC/libf/phystd/volcano.F90
r3436 r3725 338 338 READ(407,*,iostat=ierr) nvoltrac 339 339 write(*,*)trim(rname)//' number of outgased tracers:',nvoltrac 340 !Read in eruption frequency (in multiples of iphysiq) and offset of eruption340 !Read in eruption frequency (in sols) and offset of eruption 341 341 READ(407,'(A)',iostat=ierr) tracline 342 342 READ(tracline,*) eruptfreq, eruptoffset, eruptdur -
trunk/LMDZ.GENERIC/libf/phystd/writediagfi.F
r2957 r3725 41 41 use surfdat_h, only: phisfi 42 42 use geometry_mod, only: cell_area 43 use time_phylmdz_mod, only: ecritphy, day_step, iphysiq, day_ini 43 use time_phylmdz_mod, only: diagfi_output_rate, dtphys, daysec 44 use time_phylmdz_mod, only: day_ini 44 45 USE mod_phys_lmdz_para, only : is_parallel, is_mpi_root, 45 46 & is_master, gather … … 73 74 REAL area((nbp_lon+1),nbp_lat) 74 75 75 integer irythme76 76 integer ierr,ierr2 77 77 integer i,j,l, ig0 … … 119 119 return 120 120 endif 121 122 !***************************************************************123 !Sortie des variables au rythme voulu124 125 irythme = int(ecritphy) ! output rate set by ecritphy126 121 127 122 !*************************************************************** … … 267 262 !------------------------------------------------------------------------ 268 263 if (nom.eq.firstnom) then 269 zitau = zitau + iphysiq264 zitau = zitau + 1 270 265 end if 271 266 … … 274 269 !-------------------------------------------------------- 275 270 276 if ( MOD(zitau+1, irythme) .eq.0.) then271 if ( MOD(zitau+1,diagfi_output_rate) .eq.0.) then 277 272 278 273 ! Compute/write/extend 'Time' coordinate (date given in days) … … 287 282 ntime=ntime+1 ! increment # of stored time steps 288 283 ! compute corresponding date (in days and fractions thereof) 289 date=(zitau +1.) /day_step284 date=(zitau +1.)*(dtphys/daysec) 290 285 ! Get NetCDF ID of 'Time' variable 291 286 ierr= NF_INQ_VARID(nid,"Time",varid) … … 619 614 endif ! of if (dim.eq.3) elseif(dim.eq.2)... 620 615 621 endif ! of if ( MOD(zitau+1, irythme) .eq.0.)616 endif ! of if ( MOD(zitau+1,diagfi_output_rate) .eq.0.) 622 617 623 618 if (is_master) then -
trunk/LMDZ.GENERIC/libf/phystd/writediagsoil.F90
r3522 r3725 20 20 use comsoil_h, only: nsoilmx, inertiedat 21 21 use geometry_mod, only: cell_area 22 use time_phylmdz_mod, only: ecritphy, day_step, iphysiq22 use time_phylmdz_mod, only: diagfi_output_rate, dtphys, daysec 23 23 use mod_phys_lmdz_para, only : is_mpi_root, is_master, gather 24 24 use mod_grid_phy_lmdz, only : klon_glo, Grid1Dto2D_glo, & … … 92 92 stop 93 93 endif 94 95 ! Set output sample rate96 isample=int(ecritphy) ! same as for diagfi outputs97 ! Note ecritphy is known from control.h98 94 99 95 ! Create output NetCDF file … … 158 154 if (name.eq.firstname) then 159 155 ! if we run across 'firstname', then it is a new time step 160 zitau=zitau+iphysiq 161 ! Note iphysiq is known from control.h 156 zitau=zitau+1 162 157 endif 163 158 164 159 ! 3. Write data, if the time index matches the sample rate 165 if (mod(zitau+1, isample).eq.0) then160 if (mod(zitau+1,diagfi_output_rate).eq.0) then 166 161 167 162 ! 3.1 If first call at this date, update 'time' variable 168 163 if (name.eq.firstname) then 169 164 ntime=ntime+1 170 date= float(zitau+1)/float(day_step)165 date=(zitau+1)*(dtphys/daysec) 171 166 ! Note: day_step is known from control.h 172 167 … … 394 389 395 390 endif ! of if (dimpx.eq.3) elseif (dimpx.eq.2) ... 396 endif ! of if (mod(zitau+1, isample).eq.0)391 endif ! of if (mod(zitau+1,diagfi_output_rate).eq.0) 397 392 398 393 ! 4. Close the NetCDF file -
trunk/LMDZ.GENERIC/libf/phystd/writediagspecIR.F
r2732 r3725 48 48 & nbp_lon, nbp_lat, grid_type, 49 49 & unstructured 50 use time_phylmdz_mod, only: ecritphy, iphysiq, day_step, day_ini 50 use time_phylmdz_mod, only: diagfi_output_rate, dtphys, daysec 51 use time_phylmdz_mod, only: day_ini 51 52 use callkeys_mod, only: iradia 52 53 … … 67 68 ! real dx0 68 69 69 integer i rythme70 integer isample 70 71 integer ierr 71 72 integer iq … … 106 107 !Sortie des variables au rythme voulu 107 108 108 i rythme = ecritphy*iradia ! sortie au rythme de ecritphy*iradia109 !EM+JL if the spe tra need to be output more frequently, need to define a ecritSpec...110 ! irythme = iphysiq ! sortie a tous les pas physique109 isample = diagfi_output_rate*iradia ! output rate 110 !EM+JL if the spectra need to be output more frequently, need to define a 111 ! diagspec_output_rate... 111 112 112 113 … … 193 194 !------------------------------------------------------------------------ 194 195 if (nom.eq.firstnom) then 195 zitau = zitau + iphysiq196 zitau = zitau + 1 196 197 end if 197 198 … … 200 201 !-------------------------------------------------------- 201 202 202 if ( MOD(zitau+1,i rythme) .eq.0.) then203 if ( MOD(zitau+1,isample) .eq.0.) then 203 204 204 205 ! Compute/write/extend 'Time' coordinate (date given in days) … … 213 214 ntime=ntime+1 ! increment # of stored time steps 214 215 ! compute corresponding date (in days and fractions thereof) 215 date= float (zitau +1)/float (day_step)216 date=(zitau +1)*(dtphys/daysec) 216 217 217 218 if (is_master) then … … 333 334 endif ! of if (dimpx.eq.3) 334 335 335 endif ! of if ( MOD(zitau+1,i rythme) .eq.0.)336 endif ! of if ( MOD(zitau+1,isample) .eq.0.) 336 337 337 338 ! Close the NetCDF file -
trunk/LMDZ.GENERIC/libf/phystd/writediagspecVI.F
r2732 r3725 48 48 & nbp_lon, nbp_lat, grid_type, 49 49 & unstructured 50 use time_phylmdz_mod, only: ecritphy, iphysiq, day_step, day_ini 50 use time_phylmdz_mod, only: diagfi_output_rate, dtphys, daysec 51 use time_phylmdz_mod, only: day_ini 51 52 use callkeys_mod, only: iradia 52 53 … … 67 68 ! real dx0 68 69 69 integer i rythme70 integer isample 70 71 integer ierr 71 72 integer iq … … 106 107 !Sortie des variables au rythme voulu 107 108 108 i rythme = ecritphy*iradia ! sortie au rythme de ecritphy109 !EM+JL if the spe tra need to be output more frequently, need to define a ecritSpec...110 ! irythme = iphysiq ! sortie a tous les pas physique109 isample = diagfi_output_rate*iradia ! output rate 110 !EM+JL if the spectra need to be output more frequently, need to define a 111 ! diagspec_output_rate... 111 112 112 113 !*************************************************************** … … 192 193 !------------------------------------------------------------------------ 193 194 if (nom.eq.firstnom) then 194 zitau = zitau + iphysiq195 zitau = zitau + 1 195 196 end if 196 197 … … 199 200 !-------------------------------------------------------- 200 201 201 if ( MOD(zitau+1,i rythme) .eq.0.) then202 if ( MOD(zitau+1,isample) .eq.0.) then 202 203 203 204 ! Compute/write/extend 'Time' coordinate (date given in days) … … 212 213 ntime=ntime+1 ! increment # of stored time steps 213 214 ! compute corresponding date (in days and fractions thereof) 214 date= float (zitau +1)/float (day_step)215 date=(zitau +1)*(dtphys/daysec) 215 216 216 217 if (is_master) then … … 332 333 endif ! of if (dimpx.eq.3) 333 334 334 endif ! of if ( MOD(zitau+1,i rythme) .eq.0.)335 endif ! of if ( MOD(zitau+1,isample) .eq.0.) 335 336 336 337 ! Close the NetCDF file
Note: See TracChangeset
for help on using the changeset viewer.