Changeset 81 for trunk/mars/libf/phymars
- Timestamp:
- Feb 28, 2011, 6:53:54 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/mars/libf/phymars/meso_physiq.F
r80 r81 1613 1613 ccc**************** WRF OUTPUT ************************** 1614 1614 ccc**************** WRF OUTPUT ************************** 1615 do ig=1,ngrid 1616 wtsurf(ig) = tsurf(ig) !! surface temperature 1617 wco2ice(ig) = co2ice(ig) !! co2 ice 1618 1619 !!! specific to WRF WRF WRF 1620 !!! just to output water ice on surface 1621 !!! uncomment the Registry entry 1622 IF (igcm_h2o_ice .ne. 0) qsurfice(ig) = qsurf(ig,igcm_h2o_ice) 1623 1624 enddo 1615 !do ig=1,ngrid 1616 ! wtsurf(ig) = tsurf(ig) !! surface temperature 1617 ! wco2ice(ig) = co2ice(ig) !! co2 ice 1618 ! 1619 ! !!! specific to WRF WRF WRF 1620 ! !!! just to output water ice on surface 1621 ! !!! uncomment the Registry entry 1622 ! IF (igcm_h2o_ice .ne. 0) qsurfice(ig) = qsurf(ig,igcm_h2o_ice) 1623 ! 1624 ! !!! "VMR_ICE" "VOL. MIXING RATIO ICE" "ppm" 1625 ! IF (igcm_h2o_ice .ne. 0) THEN 1626 ! vmr=zq(1:ngridmx,1:nlayermx,igcm_h2o_ice)*mugaz/mmol(igcm_h2o_ice) 1627 ! ENDIF 1628 ! 1629 !enddo 1630 wtsurf(1:ngrid) = tsurf(1:ngrid) !! surface temperature 1631 wco2ice(1:ngrid) = co2ice(1:ngrid) !! co2 ice 1632 mtot(1:ngrid) = mtot(1:ngrid) * 1.e6 / 917. 1633 IF (igcm_h2o_ice .ne. 0) THEN 1634 qsurfice(1:ngrid) = qsurf(1:ngrid,igcm_h2o_ice) 1635 vmr=1.e6 * zq(1:ngrid,1:nlayer,igcm_h2o_ice) 1636 . * mugaz / mmol(igcm_h2o_ice) 1637 ENDIF 1638 1625 1639 c 1626 1640 c THIS INCLUDE IS AUTOMATICALLY GENERATED FROM REGISTRY
Note: See TracChangeset
for help on using the changeset viewer.