Changeset 1202 for trunk/MESOSCALE/LMD_MM_MARS/SRC
- Timestamp:
- Mar 7, 2014, 11:51:22 AM (11 years ago)
- Location:
- trunk/MESOSCALE/LMD_MM_MARS/SRC
- Files:
-
- 5 edited
-
WPS/wps_mars/geogrid/GEOGRID.TBL.ARW_Mars (modified) (6 diffs)
-
WRFV2/mars_lmd/libf/phymars/inifis.F (modified) (1 diff)
-
WRFV2/mars_lmd/libf/phymars/physiq.F (modified) (4 diffs)
-
WRFV2/mars_lmd/libf/phymars/surfini.F (modified) (1 diff)
-
WRFV2/meso_callkeys.h (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/MESOSCALE/LMD_MM_MARS/SRC/WPS/wps_mars/geogrid/GEOGRID.TBL.ARW_Mars
r1198 r1202 29 29 interp_option = tyler:four_pt 30 30 rel_path= tyler:mola_topo32 31 interp_option = tylerhires:four_pt 32 rel_path= tylerhires:mola_topo64 31 33 =============================== 32 34 name = HGT_U … … 55 57 interp_option = tyler:four_pt 56 58 rel_path= tyler:mola_topo32 59 interp_option = tylerhires:four_pt 60 rel_path= tylerhires:mola_topo64 57 61 =============================== 58 62 name = HGT_V … … 81 85 interp_option = tyler:four_pt 82 86 rel_path= tyler:mola_topo32 87 interp_option = tylerhires:four_pt 88 rel_path= tylerhires:mola_topo64 83 89 =============================== 84 90 name = THERMAL_INERTIA … … 106 112 interp_option = tyler:sixteen_pt+four_pt 107 113 rel_path = tyler:tylerall/ti/ 114 interp_option = tylerhires:sixteen_pt+four_pt 115 rel_path = tylerhires:tylercap 108 116 =============================== 109 117 name = ALBEDO_GCM … … 131 139 interp_option = tyler:sixteen_pt+four_pt 132 140 rel_path = tyler:tylerall/albedo/ 141 interp_option = tylerhires:sixteen_pt+four_pt 142 rel_path = tylerhires:tylercap 133 143 =============================== 134 144 name = SURF_ROUGHNESS … … 154 164 fill_missing=0.12 ## north polar 155 165 interp_option = tyler:nearest_neighbor 156 rel_path= tyler:rough_GCM/ 166 rel_path= tyler:rough_GCM/ 167 interp_option = tylerhires:nearest_neighbor 168 rel_path= tylerhires:rough_GCM/ 157 169 #=============================== 158 170 #name=LANDUSEF -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/inifis.F
r331 r1202 504 504 endif 505 505 506 READ(99,*) 507 READ(99,*) 508 509 !! prescribed crocus line for caps 510 READ(99,fmt='(a)') ch1 511 READ(99,*,iostat=ierr) tituscap 512 if (ierr.ne.0) then 513 write(*,*) 'no titus cap requested. set it to false' 514 tituscap = .false. 515 endif 516 WRITE(*,8000) ch1,tituscap 517 506 518 PRINT*,'--------------------------------------------' 507 519 PRINT* -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/physiq.F
r1180 r1202 330 330 integer igmin, lmin 331 331 logical tdiag 332 logical tituscap333 332 334 333 REAL zplev(ngrid,nlayermx+1),zplay(ngrid,nlayermx) … … 345 344 346 345 c======================================================================= 346 347 347 348 348 349 … … 1064 1065 c ------------------------------------------- 1065 1066 1066 tituscap = .true. 1067 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1068 !!! get the actual co2 seasonal cap from Titus observations 1069 !!! (inherited from GCM only at first step, but not a big deal) 1067 1070 IF (tituscap) THEN 1068 PRINT*,'TITUS CAP !!!'1069 !!! get the actual co2 seasonal cap from Titus observations1070 1071 CALL geticecover( ngrid, 180.*zls/pi, 1071 1072 . 180.*long/pi, 180.*lati/pi, co2ice ) 1072 co2ice = co2ice * 10000. 1073 co2ice(:) = co2ice(:)*10000. 1074 emis(:) = 0.95 ! so that points outside the cap are indeed at 0.95 1075 ! avoid unwanted patchiness from GCM initial state 1073 1076 ENDIF 1077 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1074 1078 1075 1079 IF(callcond) THEN … … 1081 1085 $ fluxsurf_sw) 1082 1086 1087 !! Titus cap security. 1088 !! -- neglect pressure changes caused by sublimation/condensation 1083 1089 IF (tituscap) THEN 1084 1090 DO ig=1,ngrid -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/surfini.F
r1041 r1202 37 37 !! si elles sont consequentes 38 38 enddo 39 if ( ( lati(ig)*180./pi .gt. 70. ) .and. 40 . ( albedodat(ig) .ge. 0.26 ) ) then 41 write(*,*)"outlier ",ig,albedodat(ig),inertiedat(ig) 39 40 !! AS: my previous stuff 41 ! if ( ( lati(ig)*180./pi .gt. 70. ) .and. 42 ! . ( albedodat(ig) .ge. 0.26 ) ) then 43 ! write(*,*)"outlier ",ig,albedodat(ig),inertiedat(ig) 44 ! watercaptag(ig) = .true. 45 ! dryness(ig) = 1. 46 ! albedodat(ig) = 0.45 !albedo_h2o_ice !! pour output 47 ! inertiedat(ig) = 800. 48 ! write(*,*)"new values ",ig,albedodat(ig),inertiedat(ig) 49 ! else 50 ! watercaptag(ig) = .false. 51 ! dryness(ig) = 1. 52 ! endif 53 ! if ( inertiedat(ig) .ge. 800. ) then 54 ! write(*,*)"change inertie from ",inertiedat(ig)," to 800." 55 ! inertiedat(ig) = 800. 56 ! endif 57 58 !! using Tyler and Barnes maps. 59 dryness(ig) = 1. 60 if ( inertiedat(ig) .ge. 600. ) then 61 write(*,*)"ice ",ig,albedodat(ig),inertiedat(ig) 42 62 watercaptag(ig) = .true. 43 dryness(ig) = 1.44 albedodat(ig) = 0.45 !albedo_h2o_ice !! pour output45 inertiedat(ig) = 800.46 write(*,*)"new values ",ig,albedodat(ig),inertiedat(ig)47 63 else 48 64 watercaptag(ig) = .false. 49 dryness(ig) = 1.50 endif51 if ( inertiedat(ig) .ge. 800. ) then52 write(*,*)"change inertie from ",inertiedat(ig)," to 800."53 inertiedat(ig) = 800.54 65 endif 55 66 -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/meso_callkeys.h
r11 r1202 19 19 . , sedimentation,iceparty,activice,water,caps 20 20 . , photochem,nqchem_min 21 . , tituscap 21 22 22 23 LOGICAL callrad,calldifv,calladj,callcond,callsoil, … … 25 26 s ,callnirco2,callnlte,callthermos,callconduct, 26 27 s calleuv,callmolvis,callmoldiff,thermochem,thermoswater 27 28 s ,tituscap 28 29 29 30 logical callemis
Note: See TracChangeset
for help on using the changeset viewer.
