Changeset 1345 for trunk/LMDZ.COMMON/libf/dyn3dpar

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/dyn3dpar/leapfrog_p.F (modified) (3 diffs)

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