| 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,ngrid |
|---|
| 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 ! (ngrid) |
|---|
| 37 | ELSE ! (caps) |
|---|
| 38 | print *,'Blork !!!' |
|---|
| 39 | print *,'caps=T avec water=F ????' |
|---|
| 40 | ENDIF ! (caps) |
|---|
