source: trunk/LMDZ.MARS/libf/phymars/planete.h @ 317

Last change on this file since 317 was 224, checked in by emillour, 13 years ago

Mars GCM:

Implemented using 'z0' roughness length map (important: 'z0' reference

field is in datafile surface.nc, which has also been updated).

  • made z0 a z0(ngridmx) array and moved 'z0' from 'planete.h' to 'surfdat.h'; added a 'z0_default' (common in surfdat.h) corresponding to the 'control' array value (contole(19) in startfi.nc).
  • adapted 'tabfi.F' to use 'z0_default'.
  • adapted 'phyetat0.F' to look for a 'z0' field in startfi.nc. If 'z0' is not found in the startfi.nc file, then the uniform default value (z0_default) is used.
  • modified 'physdem1.F' to write 'z0' field to restart.nc
  • adapted use of z0() in 'physiq.F' (diagnostic computation of surface stress), 'vdifc.F' and 'vdif_cd.F'.
  • adapted 'dustdevil.F' to use 'z0_default'.
  • 'testphys1d.F' now uses 'z0_default', and the value to use can be set in run.def (with "z0=TheValueYouWant?").
  • modified 'datareadnc.F' to load reference map of 'z0' from surface.nc, and added a 'z0' option in 'newstart.F' to force a uniform value of z0. Note that the use of the z0 map is automatic when using newstart, but only when it loads a start_archive.nc file.
File size: 706 bytes
Line 
1c-----------------------------------------------------------------------
2c INCLUDE planet.h
3
4      COMMON/planet/aphelie,periheli,year_day,peri_day,                 &
5     &       obliquit,                                                  &
6     &       lmixmin,emin_turb,coefvis,coefir,                          &
7     &       timeperi,e_elips,p_elips,unitastr
8
9      REAL aphelie,periheli,year_day,peri_day,                          &
10     &     obliquit,                                                    &
11     &     lmixmin,emin_turb,coefvis,coefir,                            &
12     &       timeperi,e_elips,p_elips,unitastr
13
14c-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.