Changeset 3593 for LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90

Timestamp:

Oct 28, 2019, 5:35:54 PM (5 years ago)

Author:

fhourdin

Message:

Significant progress for the SCM standard format

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90

@@ -818 +818 @@
        endif
 !
-!---------------------------------------------------------------------
-!    Initialize target profile for RHT nudging if needed
-!---------------------------------------------------------------------
-      if (nudge(inudge_RHT)) then
-        call nudge_RHT_init(plev,play,temp,q(:,1),t_targ,rh_targ)
-      endif
-      if (nudge(inudge_UV)) then
-        call nudge_UV_init(plev,play,u,v,u_targ,v_targ)
-      endif
-!
 !=====================================================================
        CALL iophys_ini
@@ -863 +853 @@
    d_q_adv(:,1)=d_q_adv(:,1)+dq_dyn(:,1)
 
+   print*,'OMEGA ',omega_mod_cas(10),omega(10)
 !---------------------------------------------------------------------
 !  Geopotential :
@@ -1017 +1008 @@
 !!!!!!!!!!!!!!!!!!!!!!!!
       d_t_nudge(:) = 0.
-      d_q_nudge(:,:) = 0.
       d_u_nudge(:) = 0.
       d_v_nudge(:) = 0.
-      if (nudge(inudge_RHT)) then
-        call nudge_RHT(timestep,plev,play,t_targ,rh_targ,temp,q(:,1),     &
-    &                  d_t_nudge,d_q_nudge(:,1))
-      endif
-      if (nudge(inudge_UV)) then
-        call nudge_UV(timestep,plev,play,u_targ,v_targ,u,v,     &
-    &                  d_u_nudge,d_v_nudge)
-      endif
+      d_q_nudge(:,:) = 0.
+      print*,'GUIDAGE ABSURDE POUR SANDU '
+      nudging_t=0.
+      nudging_qv=0.
+      IF ( nint(nudging_t) /= 0 ) d_t_nudge(:)=(temp_nudg_mod_cas(:)-temp(:))/nudging_t
+      IF ( nint(nudging_qv) /= 0 ) d_q_nudge(:,1)  =(qv_nudg_mod_cas(:)  -q(:,1) )/nudging_qv
+      IF ( nint(nudging_u*nudging_v) /= 0 ) THEN
+              d_u_nudge(:)=(u_nudg_mod_cas(:)-u(:))/nudging_u
+              d_v_nudge(:)=(v_nudg_mod_cas(:)-v(:))/nudging_u
+      ENDIF
+      !CALL iophys_ecrit('qv',klev,'qv','g/kg',1000*q(:,1))
+      !CALL iophys_ecrit('qvnud',klev,'qvnud','g/kg',1000*u_nudg_mod_cas)
+      !CALL iophys_ecrit('u',klev,'u','m/s',u)
+      !CALL iophys_ecrit('unud',klev,'unud','m/s',u_nudg_mod_cas)
+      !CALL iophys_ecrit('v',klev,'v','m/s',v)
+      !CALL iophys_ecrit('vnud',klev,'vnud','m/s',v_nudg_mod_cas)
+
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!