Changeset 3579 for LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90

Timestamp:

Oct 9, 2019, 3:11:07 PM (5 years ago)

Author:

Laurent Fairhead

Message:

Make aquaplanets run again (on jean-zay)
EM & MP

File:

• LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90

@@ -8 +8 @@
 CONTAINS
 
-  SUBROUTINE iniaqua(nlon, iflag_phys)
+  SUBROUTINE iniaqua(nlon,year_len,iflag_phys)
 
     ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -37 +37 @@
     USE mod_grid_phy_lmdz
     USE ioipsl_getin_p_mod, ONLY : getin_p
-    USE phys_cal_mod , ONLY: year_len
+    USE phys_cal_mod , ONLY: calend, year_len_phy => year_len 
     IMPLICIT NONE
 
@@ -44 +44 @@
     include "dimsoil.h"
 
-    INTEGER, INTENT (IN) :: nlon, iflag_phys
+    INTEGER, INTENT (IN) :: nlon, year_len, iflag_phys
     ! IM ajout latfi, lonfi
 !    REAL, INTENT (IN) :: lonfi(nlon), latfi(nlon)
@@ -129 +129 @@
     ! -------------------------------
 
-
+    !IF (calend .EQ. "earth_360d") Then
+      year_len_phy = year_len
+    !END IF
+    
     if (year_len.ne.360) then
+      write (*,*) year_len
       write (*,*) 'iniaqua: 360 day calendar is required !'
       stop
@@ -351 +355 @@
     alp_bl =0.
     treedrg(:,:,:)=0.
+
+    u10m = 0.
+    v10m = 0.
+
+    ql_ancien   = 0.
+    qs_ancien   = 0.
+    u_ancien    = 0.
+    v_ancien    = 0.
+    prw_ancien  = 0.
+    prlw_ancien = 0.
+    prsw_ancien = 0.  
+
+    ale_wake    = 0.
+    ale_bl_stat = 0.  
+
 
 !ym error : the sub surface dimension is the third not second : forgotten for iniaqua