source: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_caps.F @ 308

Last change on this file since 308 was 285, checked in by aslmd, 14 years ago

08/09/11 == AS

LMDZ.MARS + MESOSCALE.
---> Setting up a more realistic water ice source at the poles (notably outliers)

surfini.F ?
Main changes and bug fixes

  • reference to comcstfi.h was wrong. big problem because e.g. pi was not known.
  • commented about a problem to be fixed, due to surfini being called before initracer.
  • MESOSCALE: put the mesoscale north cap definition into a precompiling flag #MESOSCALE

for the moment: if [alb_mean_TES > 0.26 and lat > 70] then outliers
(previously done in meso_inc_caps.F)

inifis.F ?
Just changed a comment with wrong formatting

--> below, only MESOSCALE

soil.F ?
if somewhere IT > IT_outliers, then makes it = IT_outliers

physiq.F ?
meso_inc/meso_inc_caps.F ?
meso_inc/meso_inc_ini.F ?
meso_inc_caps no longer called. keep for reference for the moment.

meso_inc/meso_inc_var.F ?
deleted lines with *_lim variables, now useless
File size: 1.7 KB
Line 
1      !!! This has to go after initracer which change dryness and watercaptag
2
3      !!!!!! MARS MESOSCALE MODELING
4      !!!!!! TEST TEST TEST TEST  AS+JBM 28/02/11
5      !!!!!! TEST TEST TEST TEST  AS+JBM 28/02/11
6      !!!!!! TEST TEST TEST TEST  AS+JBM 28/02/11
7      !!!!
8      !!!! principe: une option 'caps=T' specifique au mesoscale
9      !!!! ... en vue d'un meso_initracer ????
10      !!!!
11      !!!! depots permanents => albedo TES du PDS
12      !!!! depots saisonniers => alb_surfice (~0.4, cf plus bas)
13      !!!!     [!!!! y compris pour les depots saisonniers sur les depots permanents]
14      !!!!
15      !!!! --> todo: il faut garder les depots saisonniers qui viennent
16      !!!!           du GCM lorsqu'ils sont consequents
17      !!!!
18      !IF ( caps .and. (igcm_h2o_ice .ne. 0) ) THEN
19
20      IF ( caps .and. water ) THEN
21          PRINT *, 'OVERWRITING watercaptag DEFINITION in INITRACER'
22          PRINT *, 'lat > lat_lim et alb > alb_lim => watercaptag=T'
23          PRINT *, 'ind for water ice: ', igcm_h2o_ice
24          !! Perennial H20 north cap defined by watercaptag=true (allows surface to be
25          !! hollowed by sublimation in vdifc).
26          do ig=1,ngridmx
27            qsurf(ig,igcm_h2o_ice)=0.  !! on jette les inputs GCM
28            if ( ( lati(ig)*180./pi .gt. 70. ) .and.
29     .           ( albedodat(ig) .ge. 0.26   ) )  then
30                    watercaptag(ig)  = .true.
31                    dryness(ig)      = 1.
32            else
33                    watercaptag(ig)  = .false.
34                    dryness(ig)      = 1.
35            endif  ! (lati, albedodat)
36          end do ! (ngridmx)
37      ELSE  ! (caps)
38          print *,'Blork !!!'
39          print *,'caps=T avec water=F ????'
40      ENDIF ! (caps)
Note: See TracBrowser for help on using the repository browser.