source: trunk/LMDZ.MARS/libf/phymars/surfdat.h @ 224

Last change on this file since 224 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: 1.2 KB
Line 
1      COMMON/surfdat/albedodat(ngridmx),                                &
2     &   phisfi(ngridmx),albedice(2),emisice(2),emissiv,                &
3     &   TESice_Ncoef,TESice_Scoef,                                     &
4     &   iceradius(2) , dtemisice(2),                                   &
5     &   zmea(ngridmx),zstd(ngridmx),                                   &
6     &   zsig(ngridmx),zgam(ngridmx),zthe(ngridmx),                     &
7     &   z0(ngridmx),z0_default
8
9      COMMON/surfdatl/TESicealbedo
10
11      real albedodat ! albedo of bare ground
12      real phisfi ! geopotential at ground level
13      real albedice ! default albedo for ice (1: North H. 2: South H.)
14      real emisice ! ice emissivity; 1:Northern hemisphere 2:Southern hemisphere
15      real emissiv ! emissivity of bare ground
16      logical TESicealbedo ! use TES ice cap albedoes (if set to .true.)
17      real TESice_Ncoef ! coefficient for TES ice albedo in Northern hemisphere
18      real TESice_Scoef ! coefficient for TES ice albedo in Southern hemisphere
19      real iceradius , dtemisice
20      real zmea,zstd,zsig,zgam,zthe
21      real z0 ! surface roughness lenght (m)
22      real z0_default ! default (constant over planet) surface roughness (m)
Note: See TracBrowser for help on using the repository browser.