Changeset 781

Timestamp:

11/21/18 10:41:53 (5 years ago)

Author:

sfromang

Message:

New vertical profil for dissipation coefficient (via vert_prof_dissip variable) for compatibility with LMDZ

File:

• codes/icosagcm/trunk/src/dissip/dissip_gcm.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/dissip/dissip_gcm.f90

@@ -87 +87 @@
   INTEGER,ALLOCATABLE   :: seed(:)
   
-            
+  ! New variables added for dissipation vertical profile (SF, 19/09/18)
+  INTEGER    :: vert_prof_dissip
+  REAL(rstd) :: dissip_factz,dissip_deltaz,dissip_zref,pseudoz
+                        
   INTEGER :: i,j,ij,ind,it,iter,M
 
@@ -429 +432 @@
     IF (is_master) PRINT *,"cdivgrad : ",cdivgrad
 
-     
-    fact=2
-    DO l=1,llm
-       IF(ap_bp_present) THEN ! height-dependent dissipation
-          zz= 1. - preff/presnivs(l)
-       ELSE
-          zz = 0.
-       END IF
-       zvert=fact-(fact-1)/(1+zz*zz)
-       tau_graddiv(l) = tau_graddiv(l)/zvert
-       tau_gradrot(l) = tau_gradrot(l)/zvert
-       tau_divgrad(l) = tau_divgrad(l)/zvert
-    ENDDO
+    ! Added CMIP5 dissipation vertical profile (SF, 19/09/18)
+    vert_prof_dissip = 1
+    CALL getin('vert_prof_dissip',vert_prof_dissip)
+    IF (vert_prof_dissip == 1) then
+       dissip_zref = 30.
+       dissip_deltaz = 10.
+       dissip_factz = 4.
+       CALL getin('dissip_factz',dissip_factz )
+       CALL getin('dissip_deltaz',dissip_deltaz )
+       CALL getin('dissip_zref',dissip_zref )
+       DO l=1,llm
+          pseudoz=8.*log(preff/presnivs(l))
+          zvert=1+ &
+               (tanh((pseudoz-dissip_zref)/dissip_deltaz)+1.)/2. &
+               *(dissip_factz-1.)
+          tau_graddiv(l) = tau_graddiv(l)/zvert
+          tau_gradrot(l) = tau_gradrot(l)/zvert
+          tau_divgrad(l) = tau_divgrad(l)/zvert
+       ENDDO
+    ELSE     
+       fact=2
+       DO l=1,llm
+          IF(ap_bp_present) THEN ! height-dependent dissipation
+             zz= 1. - preff/presnivs(l)
+          ELSE
+             zz = 0.
+          END IF
+          zvert=fact-(fact-1)/(1+zz*zz)
+          tau_graddiv(l) = tau_graddiv(l)/zvert
+          tau_gradrot(l) = tau_gradrot(l)/zvert
+          tau_divgrad(l) = tau_divgrad(l)/zvert
+       ENDDO
+    ENDIF
 
     mintau=tau_graddiv(1)