Changeset 80
- Timestamp:
- Feb 28, 2011, 2:54:20 PM (14 years ago)
- Location:
- trunk
- Files:
-
- 3 added
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/000-USERS
r74 r80 17 17 **** Name **** **** Institute **** **** Admitted **** **** Use **** 18 18 Arnaud Colaitis LMD [CNRS] 26 - 01 - 2011 Mars GCM, MMM and LES 19 Franck Lefevre LATMOS [CNRS] 25 - 02 - 2011 Mars GCM 20 Francisco Gonzalez-Galindo IAA [CSIC, Spain] 28 - 02 - 2011 Mars GCM -
trunk/mars/libf/phymars/meso_physiq.F
r77 r80 432 432 ! PRINT*, tnom(iq), pq(:,:,iq) 433 433 !ENDDO 434 434 435 c 435 436 c ****WRF … … 482 483 CALL initracer(qsurf,co2ice) 483 484 ENDIF ! end tracer 485 486 !!!!!! WRF WRF WRF MARS MARS 487 !!!!!! TEST TEST TEST TEST AS+JBM 28/02/11 488 !!!!!! TEST TEST TEST TEST AS+JBM 28/02/11 489 !!!!!! TEST TEST TEST TEST AS+JBM 28/02/11 490 !!!! 491 !!!! principe: une option 'caps=T' specifique au mesoscale 492 !!!! ... en vue d'un meso_initracer ???? 493 !!!! 494 !!!! depots permanents => albedo TES du PDS 495 !!!! depots saisonniers => alb_surfice (~0.4, cf plus bas) 496 !!!! [!!!! y compris pour les depots saisonniers sur les depots permanents] 497 !!!! 498 !!!! --> todo: il faut garder les depots saisonniers qui viennent 499 !!!! du GCM lorsqu'ils sont consequents 500 !!!! 501 IF ( caps .and. (igcm_h2o_ice .ne. 0) ) THEN 502 PRINT *, 'OVERWRITING watercaptag DEFINITION in INITRACER' 503 PRINT *, 'lat>70 et alb>0.26 => watercaptag=T' 504 !! Perennial H20 north cap defined by watercaptag=true (allows surface to be 505 !! hollowed by sublimation in vdifc). 506 do ig=1,ngridmx 507 qsurf(ig,igcm_h2o_ice)=0. !! on jette les inputs GCM 508 if ( (lati(ig)*180./pi.gt.70.) .and. 509 . (albedodat(ig).ge.0.26) ) then 510 watercaptag(ig)=.true. 511 dryness(ig) = 1. 512 else 513 watercaptag(ig)=.false. 514 dryness(ig) = 1. 515 endif ! (lati, albedodat) 516 end do ! (ngridmx) 517 ELSE ! (caps) 518 print *,'Blork !!!' 519 print *,'caps=T avec water=F ????' 520 ENDIF ! (caps) 521 !!!!!! TEST TEST TEST TEST AS+JBM 28/02/11 522 !!!!!! TEST TEST TEST TEST AS+JBM 28/02/11 523 !!!!!! TEST TEST TEST TEST AS+JBM 28/02/11 484 524 485 525 -
trunk/mesoscale/LMD_MM_MARS/SIMU/runmeso
r78 r80 21 21 ########################################### 22 22 ########################################### 23 24 datafile = '/u/forget/WWW/datagcm/datafile'25 23 26 24 changeregis=0 ## 1 if changed registry … … 517 515 cd ${simu_folder} 518 516 rm -rf ${simu_folder}/wrfout* 519 ln -sf ${datafile} datafile520 517 if [[ ${compilo} -eq 3 ]] 521 518 then -
trunk/mesoscale/TESTS/newphys_tracers/namelist.input
r77 r80 20 20 21 21 &domains 22 time_step = 50 !! Dynamical timestep !!! BAISSER NE CHANGE RIEN22 time_step = 20 !! Dynamical timestep 23 23 dx = 30000, !! (p2) Horizontal resolution 24 24 dy = 30000, !! (p2) Horizontal resolution (should be equal to dx) 25 e_we = 61,!! (r)(p2) Number of longitude grid points26 e_sn = 61,!! (r)(p2) Number of latitude grid points25 e_we = 121, !! (r)(p2) Number of longitude grid points 26 e_sn = 121, !! (r)(p2) Number of latitude grid points 27 27 e_vert = 61, !! (r)(p2) Number of vertical levels 28 28 p_top_requested = 5 !! (p3) Chosen value of pressure at the top of the model … … 34 34 35 35 &physics 36 mars = 11 ! 0 1 2 3 1136 mars = 11 37 37 num_soil_layers = 18 !!! NOUVELLE PHYSIQUE 38 38 !! … … 44 44 tiso = 100. 45 45 pd_scalar = T !!! tres IMPORTANT, sinon valeurs negatives 46 !! 47 !! OPTIONAL48 !! 49 non_hydrostatic = F, !! Integrate in non-hydrostatic/hydrostatic mode50 !! --- necessaire pour garder un pas de tps long51 !! --- et parce qu'un bout d'Arsia est au bord46 !!! 47 !!! OPTIONAL 48 !!! 49 !non_hydrostatic = F, !! Integrate in non-hydrostatic/hydrostatic mode 50 ! !! --- necessaire pour garder un pas de tps long 51 ! !! --- et parce qu'un bout d'Arsia est au bord 52 52 / 53 53 -
trunk/mesoscale/TESTS/newphys_tracers/namelist.wps
r76 r80 16 16 i_parent_start = 0, !! (n)(p2) x-position of the bottom-left nest corner 17 17 j_parent_start = 0, !! (n)(p2) y-position of the bottom-left nest corner 18 e_we = 61, !! (r)(p2) Number of longitude grid points19 e_sn = 61, !! (r)(p2) Number of latitude grid points18 e_we = 121, !! (r)(p2) Number of longitude grid points 19 e_sn = 121, !! (r)(p2) Number of latitude grid points 20 20 dx = 30000 21 21 dy = 30000 … … 24 24 !! NB: please glance at geogrid/GEOGRID.TBL 25 25 map_proj = 'mercator', !! map projection: 'mercator', 'lambert' or 'polar' 26 ref_lat = 7., !! north latitude of the center of the domain27 ref_lon = -1 09., !! east longitude of the center of the domain28 truelat1 = 7., !! (lambert or polar) lat position of projection cone26 ref_lat = 10., !! north latitude of the center of the domain 27 ref_lon = -125., !! east longitude of the center of the domain 28 truelat1 = 10., !! (lambert or polar) lat position of projection cone 29 29 truelat2 = 0.0, !! [do not modify] 30 stand_lon = -1 09., !! (lambert or polar) lon position of projection cone30 stand_lon = -125., !! (lambert or polar) lon position of projection cone 31 31 geog_data_path = './WPS_GEOG', !! [do not modify: symbolic link in the WPS folder] 32 32 /
Note: See TracChangeset
for help on using the changeset viewer.