Changeset 1345 for trunk/LMDZ.COMMON

Timestamp:

Sep 8, 2014, 10:09:19 AM (10 years ago)

Author:

aslmd

Message:

LMDZ.COMMON. introduced an attenuation profile for nudging towards wind. new parameter: phalf: pressure at which attenuation factor is divided by two.

Location:

trunk/LMDZ.COMMON/libf

Files:

• dyn3d_common/comuforc_h.F90 (modified) (1 diff)
• dyn3d_common/relaxetat0.F (modified) (2 diffs)
• dyn3dpar/leapfrog_p.F (modified) (3 diffs)

trunk/LMDZ.COMMON/libf/dyn3d_common/comuforc_h.F90

@@ -2 +2 @@
        IMPLICIT NONE
        REAL,ALLOCATABLE,DIMENSION(:,:) :: uforc
+       REAL,ALLOCATABLE,DIMENSION(:) :: attenua
        REAL :: facwind
+       REAL :: phalf
        end module comuforc_h
 

trunk/LMDZ.COMMON/libf/dyn3d_common/relaxetat0.F

@@ -22 +22 @@
 ! VARIABLES
       REAL ucov(iip1, jjp1,llm)
+      REAL press(llm)
       INTEGER ierr, nid, nvarid
       INTEGER i,j,l,ij
@@ -63 +64 @@
       CALL getin('facwind',facwind)
 
+      !! CALCULATE ATTENUATION FACTOR WITH ALTITUDE
+      ierr = NF_INQ_VARID (nid, "presnivs", nvarid)
+      IF (ierr .NE. NF_NOERR) THEN
+         write(lunout,*)"dynetat0: Le champ <presnivs> est absent"
+         CALL abort
+      ENDIF
+      ierr = nf90_get_var(nid, nvarid, press)
+      IF (ierr .NE. NF_NOERR) THEN
+         write(lunout,*)"dynetat0: Lecture echouee pour <press>"
+         CALL abort
+      ENDIF
+      phalf = 1.e5 !! at this pressure, attenua is 0.5
+      CALL getin('phalf',phalf)
+      !! attenua goes from 1 at lowermost levels to 0 at uppermost levels
+      IF (.NOT.ALLOCATED(attenua)) ALLOCATE(attenua(llm))
+      attenua = 0.5*(1.+tanh(log10(press/phalf)))
+
       !! END
       ierr = NF_CLOSE(nid)

trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F

@@ -678 +678 @@
 
 #ifdef NODYN
-      WRITE(lunout,*)"NO DYN !!!!!"
+      !WRITE(lunout,*)"NO DYN !!!!!"
       dv(:,:) = 0.D+0
       du(:,:) = 0.D+0
@@ -684 +684 @@
       dq(:,:,:) = 0.D+0
       dp(:) = 0.D+0
+      if (planet_type.eq."generic") then
+      if (ok_guide) then
+       DO l=1,llm
+       ucov(:,l) = ucov(:,l) + dt*attenua(l)*
+     &   ((uforc(:,l)-ucov(:,l))/facwind)
+       ENDDO
+      endif
+      endif
 #else
 ! ADAPTATION GCM POUR CP(T)
@@ -719 +727 @@
       if (planet_type.eq."generic") then
        if (ok_guide) then
-         du(:,:) = du(:,:) + ((uforc(:,:)-ucov(:,:)) / facwind)
+         DO l=1,llm
+          du(:,l)=du(:,l)+attenua(l)*((uforc(:,l)-ucov(:,l))/facwind) 
+         ENDDO
        endif
       endif