Changeset 1164

Timestamp:

May 27, 2009, 12:35:49 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Reglage d'une configuration a 39 niveaux pour la discretisation verticale FH
LF

Location:

LMDZ4/branches/LMDZ4-dev/libf

Files:

• dyn3d/disvert.F (modified) (3 diffs)
• dyn3dpar/disvert.F (modified) (3 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/disvert.F

@@ -36 +36 @@
 c
       INTEGER l
-      REAL snorm
+      REAL snorm,dsigmin
       REAL alpha,beta,gama,delta,deltaz,h
       INTEGER np,ierr
@@ -96 +96 @@
       WRITE(LUNOUT,*)'WARNING!!! Ancienne discretisation verticale'
 
+      if (ok_strato) then
+         if (llm==39) then
+            dsigmin=0.3
+         else if (llm==50) then
+            dsigmin=1.
+         else
+            WRITE(LUNOUT,*) 'ATTENTION discretisation z a ajuster'
+            dsigmin=1.
+         endif
+         WRITE(LUNOUT,*) 'Discretisation verticale DSIGMIN=',dsigmin
+      endif
+
       h=7.
       snorm  = 0.
@@ -149 +161 @@
 c
       ENDDO
+
+      bp(1)=1.
+      ap(1)=0.
+
       ap(llmp1) = pa * ( sig(llmp1) - bp(llmp1) )
 

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/disvert.F

@@ -36 +36 @@
 c
       INTEGER l
-      REAL snorm
+      REAL snorm,dsigmin
       REAL alpha,beta,gama,delta,deltaz,h
       INTEGER np,ierr
@@ -96 +96 @@
       WRITE(LUNOUT,*)'WARNING!!! Ancienne discretisation verticale'
 
+      if (ok_strato) then
+         if (llm==39) then
+            dsigmin=0.3
+         else if (llm==50) then
+            dsigmin=1.
+         else
+            WRITE(LUNOUT,*) 'ATTENTION discretisation z a ajuster'
+            dsigmin=1.
+         endif
+         WRITE(LUNOUT,*) 'Discretisation verticale DSIGMIN=',dsigmin
+      endif
+
       h=7.
       snorm  = 0.
@@ -149 +161 @@
 c
       ENDDO
+
+      bp(1)=1.
+      ap(1)=0.
+
       ap(llmp1) = pa * ( sig(llmp1) - bp(llmp1) )