Changeset 3223 for LMDZ6/trunk/libf/phylmd/dyn1d/lmdz1d.F90

Timestamp:

Feb 17, 2018, 6:02:23 PM (6 years ago)

Author:

fhourdin

Message:

Corrections forcages cas 1D sandu

File:

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

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

@@ -286 +286 @@
        hthturb_gcssold = 0.
        hqturb_gcssold = 0.
+
+
+
 
 !---------------------------------------------------------------------
@@ -985 +988 @@
 !
 !=====================================================================
+       CALL iophys_ini
 ! START OF THE TEMPORAL LOOP :
 !=====================================================================
@@ -1109 +1113 @@
       endif     
 !RC
+      if (forcing_sandu) then
+         ug(1:llm)=u_mod(1:llm)
+         vg(1:llm)=v_mod(1:llm)
+      endif
 
       IF (prt_level >= 5) print*, 'fcoriolis, xlat,mxcalc ', &
                                    fcoriolis, xlat,mxcalc
 
-       du_age(1:mxcalc)=fcoriolis*(v(1:mxcalc)-vg(1:mxcalc))
-       dv_age(1:mxcalc)=-fcoriolis*(u(1:mxcalc)-ug(1:mxcalc))
 !       print *,'u-ug=',u-ug
 
 !!!!!!!!!!!!!!!!!!!!!!!!
 ! Geostrophic wind
+! Le calcul ci dessous est insuffisamment precis
+!      du_age(1:mxcalc)=fcoriolis*(v(1:mxcalc)-vg(1:mxcalc))
+!      dv_age(1:mxcalc)=-fcoriolis*(u(1:mxcalc)-ug(1:mxcalc))
 !!!!!!!!!!!!!!!!!!!!!!!!
        sfdt = sin(0.5*fcoriolis*timestep)