Changeset 1959 for LMDZ5/trunk/libf/dyn3d_common

Timestamp:

Feb 7, 2014, 10:37:38 AM (11 years ago)

Author:

fhourdin

Message:

Nouvelles options dans disvert.F90 :
=------------------------------------

1. avec vert_sampling="sigma" -> niveaux sigma plutôt qu'hybride.

P_i,j,k= Bhyb_k * Ps , Ahyb_k=0.

2. disgmin, l'épaisseur approximative du premier niveau en coordonnée

sigma peut être défini dans les fichiers .def.
Cette option est valable pour vert_sampling=sigma, tropo ou strato.
Valeur par défaut : 1, sauf pour (llm==39 .and. vert_sampling=="strato")

ce qui correspond aux configurations CMIP5

=======================================================================
New options in disvert.F90 :
=---------------------------

1. with vert_samplig="sigma" -> sigma levels instead of hybrid

P_i,j,k= Bhyb_k * Ps , Ahyb_k=0.

2. disgmin, the approximate thickness of the first layer in sigma

coordinate can be specified from the .def files.
Active for vert_sampling=sigma, tropo or strato.
Defaut value : 1, except for (llm==39 .and. vert_sampling=="strato")

which corresponds to the CMIP5 configurations.

--Cette ligne, et les suivantes ci-dessous, seront ignorées-

M dyn3d_common/disvert.F90

File:

• LMDZ5/trunk/libf/dyn3d_common/disvert.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/dyn3d_common/disvert.F90

@@ -50 +50 @@
   call getin('vert_sampling', vert_sampling)
   print *, 'vert_sampling = ' // vert_sampling
+  if (llm==39 .and. vert_sampling=="strato") then
+     dsigmin=0.3 ! Vieille option par défaut pour CMIP5
+  else
+     dsigmin=1.
+  endif
+  call getin('dsigmin', dsigmin)
+  WRITE(LUNOUT,*) trim(modname), 'Discretisation verticale DSIGMIN=',dsigmin
+
 
   select case (vert_sampling)
@@ -86 +94 @@
 
      ap = pa * (sig - bp)
+  case("sigma")
+     DO l = 1, llm
+        x = 2*asin(1.) * (l - 0.5) / (llm + 1)
+        dsig(l) = dsigmin + 7.0 * SIN(x)**2
+     ENDDO
+     dsig = dsig / sum(dsig)
+     sig(llm+1) = 0.
+     DO l = llm, 1, -1
+        sig(l) = sig(l+1) + dsig(l)
+     ENDDO
+
+     bp(1)=1.
+     bp(2: llm) = sig(2:llm)
+     bp(llmp1) = 0.
+     ap(:)=0.
   case("tropo")
      DO l = 1, llm
         x = 2*asin(1.) * (l - 0.5) / (llm + 1)
-        dsig(l) = 1.0 + 7.0 * SIN(x)**2
+        dsig(l) = dsigmin + 7.0 * SIN(x)**2
      ENDDO
      dsig = dsig / sum(dsig)
@@ -104 +127 @@
      ap(2: llm + 1) = pa * (sig(2: llm + 1) - bp(2: llm + 1))
   case("strato")
-     if (llm==39) then
-        dsigmin=0.3
-     else if (llm==50) then
-        dsigmin=1.
-     else
-        write(lunout,*) trim(modname), ' ATTENTION discretisation z a ajuster'
-        dsigmin=1.
-     endif
-     WRITE(LUNOUT,*) trim(modname), 'Discretisation verticale DSIGMIN=',dsigmin
-
      DO l = 1, llm
         x = 2*asin(1.) * (l - 0.5) / (llm + 1)