Changes between Version 1 and Version 2 of Documentation/11layer

Timestamp:

2013-11-29T08:46:03+01:00 (11 years ago)

Author:

mmcgrath

Comment:

--

Documentation/11layer

@@ -1 @@
-You can find instructions on how to install the coupled LMDZOR model here.  This page assumes you already have the model checked out, and simply want to active the 11 hydrology, stomate, and/or the new physics.  This all comes from tests done by Fuxing Wang.
+You can find instructions on how to install the coupled LMDZOR model here.  This page assumes you already have the model checked out, and simply want to active the 11 hydrology, stomate, and/or the new physics.  This all comes from tests done by Fuxing Wang without stomate.  If you are not using stomate, you need to use an LAI map.
 
 Most of the changes are done in the configuration files after the model is compiled, with the exception of the following.
 
+== Source code changes ==
 
 '''New Physics'''
@@ -34 +35 @@
 }}}
 
+'''11 layer hydrology (CWRR)'''
+
+In modeles/LMDZ/libf/phylmd/surf_land_orchidee_mod.F90:
+
+Original
+
+{{{
+    zlev(1:knon) = (100.*plev(1:knon))/((ps(1:knon)/RD*temp_air(1:knon))*RG)
+}}}
+
+Modified
+{{{
+    zlev(1:knon) = plev(1:knon)*RD*temp_air(1:knon)/((ps(1:knon)*100.0)*RG)
+}}}
+
+In modeles/LMDZ/libf/phylmd/coefcdrag.F90:
+
+Original
+{{{
+       ztsolv(i) = ts(i)
+       ztvd(i) = t(i) * (psol(i)/pref(i))**RKAPPA
+       trm0(i) = 1. + RETV * max(qsurf(i),0.0)
+       trm1(i) = 1. + RETV * max(q(i),0.0)
+       ztsolv(i) = ztsolv(i) * trm0(i)
+       ztvd(i) = ztvd(i) * trm1(i)
+       zri1(i) = zdphi(i)*(ztvd(i)-ztsolv(i))/(zdu2(i)*ztvd(i))
+ 
+}}}
+
+Modified
+{{{
+       ztsolv(i) = ts(i)
+!       ztvd(i) = t(i) * (psol(i)/pref(i))**RKAPPA
+       trm0(i) = 1. + RETV * max(qsurf(i),0.0)
+       trm1(i) = 1. + RETV * max(q(i),0.0)
+       ztsolv(i) = ztsolv(i) * trm0(i)
+!       ztvd(i) = ztvd(i) * trm1(i)
+       ztvd(i) = (t(i)+zdphi(i)/RCPD/(1.+RVTMP2*q(i))) &
+          *(1.+RETV*q(i))
+       zri1(i) = zdphi(i)*(ztvd(i)-ztsolv(i))/(zdu2(i)*ztvd(i))
+}}}
+
+Once you have done this, you can create the makefile and compile as normal (found elsewhere on the Wiki).
+
+== Configuration Files ==
+
+These changes are done in the directory where you will be submitting the job (config/LMDZOR_v5.2/JOBNAME, for example).
+
+''' 11 layer hydrology (CWRR)'''
+
+Confirm the following lines are present.  If a line exists but with something different after the = sign, then change it.
+
+In PARAM/orchidee.def_CWRR:
+{{{
+   RIVER_ROUTING = n
+}}}
+
+In COMP/orchidee.card:
+{{{
+   DefSuffix=CWRR
+}}}
+
+''' New physics'''
+
+Confirm the following lines are present.
+
+In PARAM/orchidee.def_CWRR:
+{{{
+   evapot_corr=y
+}}}
+
+In COMP/lmdz.card:
+{{{
+   LMDZ_Physics=NPv3.1
+}}}
+
+In PARAM/physiq.def_L39_NPv3.1 (assuming you are using 39 levels, i.e. you compile with gmake LMD????-L39):
+{{{
+ iflag_thermals=17
+ iflag_pbl=11 
+}}}
+
+''' STOMATE '''
+
+In PARAM/orchidee.def_CWRR:
+{{{
+STOMATE_OK_STOMATE = y
+}}}
+
+In config.card:
+{{{
+SBG= (stomate, ORCHIDEE_1_9_5)
+}}}
+
+In stomate.card (useful mostly for debugging):
+{{{
+[OutputText]
+List=   (out_orchidee, run.def)
+}}}
+
+''' No STOMATE '''
+
+In PARAM/orchidee.def_CWRR:
+{{{
+LAI_MAP = y
+STOMATE_OK_STOMATE = n
+}}}
+
+In COMP/orchidee.card:
+{{{
+ListNonDel= (${R_BC}/SRF/${config_UserChoices_TagName}/lai2D.nc, .)
+}}}
+
+In config.card, remove the following line:
+{{{
+SBG= (stomate, ORCHIDEE_1_9_5)
+}}}