Changeset 1345 for trunk/LMDZ.COMMON/libf/dyn3d_common/relaxetat0.F

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.

File:

• trunk/LMDZ.COMMON/libf/dyn3d_common/relaxetat0.F (modified) (2 diffs)

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)