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

Timestamp:

Oct 29, 2019, 7:32:03 AM (5 years ago)

Author:

fhourdin

Message:

Improvement for SCM standard format

File:

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

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

@@ -1011 +1011 @@
       d_v_nudge(:) = 0.
       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)
+      DO l=1,llm
+         IF ( play(l) < p_nudging_u .AND. nint(nudging_u) /= 0 ) &
+             & d_u_nudge(l)=(u_nudg_mod_cas(l)-u(l))/nudging_u
+         IF ( play(l) < p_nudging_v .AND. nint(nudging_v) /= 0 ) &
+             & d_v_nudge(l)=(v_nudg_mod_cas(l)-v(l))/nudging_v
+         IF ( play(l) < p_nudging_t .AND. nint(nudging_t) /= 0 ) &
+             & d_t_nudge(l)=(temp_nudg_mod_cas(l)-temp(l))/nudging_t
+         IF ( play(l) < p_nudging_qv .AND. nint(nudging_qv) /= 0 ) &
+             & d_q_nudge(l,1)=(qv_nudg_mod_cas(l)-q(l,1))/nudging_qv
+      ENDDO
+      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)
+      CALL iophys_ecrit('dtnud',klev,'dtnud','K/day',86400*d_t_nudge)
+      CALL iophys_ecrit('dqnud',klev,'dqnud','K/day',1000*86400*d_q_nudge(:,1))
 
 !
@@ -1109 +1112 @@
      &             +dt_cooling(1:mxcalc))  ! Taux de chauffage ou refroid.
 
-       print*,'MXCALC d_t_adv ',mxcalc,d_t_adv(1:20)*86400
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!