Changeset 131 for trunk/mars/libf/phymars/meso_physiq.F
- Timestamp:
- May 25, 2011, 4:19:36 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/mars/libf/phymars/meso_physiq.F
r86 r131 1060 1060 DO ig=1,ngrid 1061 1061 zdtsurf(ig) = zdtsurf(ig) + zdtsurfc(ig) 1062 !!!**WRF: newphys: ici la pression n'est plus mise a jour ds le GCM1063 !!!**WRF: mais il faut retablir ca dans le cas du mesoscale ?1064 !!!**WRF: ...non probablement OK1065 ! ps(ig)=pplev(ig,1) + pdpsrf(ig)*ptimestep1066 1062 ENDDO 1067 1063 … … 1091 1087 c ---------------------------------------- 1092 1088 IF (water) THEN 1093 1094 c **WRF: new arguments here rnuclei,rice,nuice1095 c plus no more iqmin +igcm_h2o_vap replaces iq, what are the consequences?1096 c checks needed when tracers simulations1097 1089 1098 1090 call watercloud(ngrid,nlayer,ptimestep, … … 1216 1208 zdqssed(1:ngrid,1:nq)=0 1217 1209 1218 c1219 c **WRF: new arguments rnuclei, rice, need checks1220 c1221 1210 call callsedim(ngrid,nlayer, ptimestep, 1222 1211 & pplev,zzlev, pt, rdust, rice, … … 1293 1282 ENDDO 1294 1283 1295 ccc1296 ccc **WRF very specific to GCM1297 ccc1298 1284 c Prescribe a cold trap at south pole (except at high obliquity !!) 1299 1285 c Temperature at the surface is set there to be the temperature … … 1325 1311 c ALWAYS PLACE these lines after newcondens !!! 1326 1312 c ------------------------------------------------------------- 1327 c1328 1313 c **WRF : OK avec le mesoscale, pas d'indices bizarres au pole 1329 c1330 1314 do ig=1,ngrid 1331 1315 if ((co2ice(ig).eq.0).and. … … 1478 1462 if (water) then 1479 1463 1480 !!1481 !!***WRF: ok, des nouveaux trucs cools de la nouvelle physique1482 !!1483 1464 call zerophys(ngrid,mtot) 1484 1465 call zerophys(ngrid,icetot)
Note: See TracChangeset
for help on using the changeset viewer.