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