Changes between Version 9 and Version 10 of Branches/Driver_Atm_Lev

Timestamp:

2012-06-05T10:20:03+02:00 (12 years ago)

Author:

jpolcher

Comment:

--

Branches/Driver_Atm_Lev

@@ -70 +70 @@
 If none of the above cases is found, the variables  '''HEIGHT_LEV1''' and '''HEIGHT_LEVW''' found in run.def will be used to fill the global fields.
 
+== Proposed convention ==
+
+We will test for the following variables in the forcing file :
+
+=== Sigma and Sigma_uv===
+
+Expected is one real value which allows to compute the height of the first level with the following formula :
+
+zlev_vec(i,j) = rau(i,j) * cte_grav * (psurf(i,j) - Sigma * psurf(i,j))
+zlevuv_vec() = rau(i,j) * cte_grav * (psurf(i,j) - Sigma_uv * psurf(i,j))
+
+Should only Sigma exist, Sigma=Sigma_uv will be assumed.
+
+=== HybSigA, HybSigB, HybSigA_uv and HybSigB_uv ===
+
+Expected is one real value for each variable. They allow to compute the height of the first level with the following formula :
+
+zlev_vec(i,j) = rau(i,j) * cte_grav * (psurf(i,j) - (HybSigA + HybSigB * psurf(i,j)))
+zlevuv_vec(i,j) = rau(i,j) * cte_grav * (psurf(i,j) - (HybSigA_uv + HybSigB_uv * psurf(i,j)))
+
+If the UV values do not exist we assume that they are the same as for T and Q.
+
+=== Levels and levels_uv === 
+
+This is the case when the levels are provied as 2D time evolving fields. The height is supposed to be in meters.
+
+The following assignation will be done for the values read at each forcing time step :
+zlev_vec(i,j) = Levles(i,j)
+
+If levels_uv is present :
+zlevuv_vec(i,j) = levels_uv(i,j)
+Else :
+zlevuv_vec(i,j) = levels(i,j)
+
+=== Height_Lev1 and Height_Levuv ===
+
+This comes back to the same case as when these variables are provided through the run.def file.
+
+The following assignations will be done :
+
+lev_vec(i,j) = Levles(i,j)
+
+If levels_uv is present :
+zlevuv_vec(i,j) = levels_uv(i,j)
+Else :
+zlevuv_vec(i,j) = levels(i,j)
+
+
 == Proposed modifications in the code ==
 
@@ -115 +163 @@
 === readdim2.f90, forcing_just_read ===
 
-Once surface pressure is read we can compute or fill zlev and zlevuv
+Once surface pressure is read we compute or fill zlev and zlevuv
 
 This will be done on the zoomed variables.