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

Timestamp:

Jun 25, 2019, 4:50:13 PM (5 years ago)

Author:

Laurent Fairhead

Message:

Modifications needed for "real" calendar in 1D model
MPL/EM

File:

• LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90 (modified) (10 diffs, 1 prop)

LMDZ6/trunk/libf/phylmd/phyaqua_mod.F90

@@ +1 @@
+!
+! $Id$
+!
 MODULE phyaqua_mod
   ! Routines complementaires pour la physique planetaire.
@@ -34 +37 @@
     USE mod_grid_phy_lmdz
     USE ioipsl_getin_p_mod, ONLY : getin_p
+    USE phys_cal_mod , ONLY: year_len
     IMPLICIT NONE
 
@@ -72 +76 @@
     CHARACTER *2 cnbl
 
-    REAL phy_nat(nlon, 360)
-    REAL phy_alb(nlon, 360)
-    REAL phy_sst(nlon, 360)
-    REAL phy_bil(nlon, 360)
-    REAL phy_rug(nlon, 360)
-    REAL phy_ice(nlon, 360)
-    REAL phy_fter(nlon, 360)
-    REAL phy_foce(nlon, 360)
-    REAL phy_fsic(nlon, 360)
-    REAL phy_flic(nlon, 360)
+    REAL phy_nat(nlon, year_len)
+    REAL phy_alb(nlon, year_len)
+    REAL phy_sst(nlon, year_len)
+    REAL phy_bil(nlon, year_len)
+    REAL phy_rug(nlon, year_len)
+    REAL phy_ice(nlon, year_len)
+    REAL phy_fter(nlon, year_len)
+    REAL phy_foce(nlon, year_len)
+    REAL phy_fsic(nlon, year_len)
+    REAL phy_flic(nlon, year_len)
 
     INTEGER, SAVE :: read_climoz = 0 ! read ozone climatology
@@ -125 +129 @@
     ! -------------------------------
 
+
+    if (year_len.ne.360) then
+      write (*,*) 'iniaqua: 360 day calendar is required !'
+      stop
+    endif
 
     type_aqua = iflag_phys/100
@@ -223 +232 @@
     ! endif !alb_ocean
 
-    DO i = 1, 360
+    DO i = 1, year_len
       ! IM Terraplanete   phy_sst(:,i) = 260.+50.*cos(rlatd(:))**2
       ! IM ajout calcul profil sst selon le cas considere (cf. FBr)
@@ -544 +553 @@
     USE mod_grid_phy_lmdz, ONLY: klon_glo
     USE mod_phys_lmdz_transfert_para, ONLY: gather
+    USE phys_cal_mod, ONLY: year_len
     IMPLICIT NONE
     include "netcdf.inc"
 
     INTEGER, INTENT (IN) :: klon
-    REAL, INTENT (IN) :: phy_nat(klon, 360)
-    REAL, INTENT (IN) :: phy_alb(klon, 360)
-    REAL, INTENT (IN) :: phy_sst(klon, 360)
-    REAL, INTENT (IN) :: phy_bil(klon, 360)
-    REAL, INTENT (IN) :: phy_rug(klon, 360)
-    REAL, INTENT (IN) :: phy_ice(klon, 360)
-    REAL, INTENT (IN) :: phy_fter(klon, 360)
-    REAL, INTENT (IN) :: phy_foce(klon, 360)
-    REAL, INTENT (IN) :: phy_flic(klon, 360)
-    REAL, INTENT (IN) :: phy_fsic(klon, 360)
-
-    REAL :: phy_glo(klon_glo, 360) ! temporary variable, to store phy_***(:)
+    REAL, INTENT (IN) :: phy_nat(klon, year_len)
+    REAL, INTENT (IN) :: phy_alb(klon, year_len)
+    REAL, INTENT (IN) :: phy_sst(klon, year_len)
+    REAL, INTENT (IN) :: phy_bil(klon, year_len)
+    REAL, INTENT (IN) :: phy_rug(klon, year_len)
+    REAL, INTENT (IN) :: phy_ice(klon, year_len)
+    REAL, INTENT (IN) :: phy_fter(klon, year_len)
+    REAL, INTENT (IN) :: phy_foce(klon, year_len)
+    REAL, INTENT (IN) :: phy_flic(klon, year_len)
+    REAL, INTENT (IN) :: phy_fsic(klon, year_len)
+
+    REAL :: phy_glo(klon_glo, year_len) ! temporary variable, to store phy_***(:)
       ! on the whole physics grid
     INTEGER :: k
@@ -665 +675 @@
 
       ! write the 'times'
-      DO k = 1, 360
+      DO k = 1, year_len
 #ifdef NC_DOUBLE
         ierr = nf_put_var1_double(nid, id_tim, k, dble(k))
@@ -809 +819 @@
   SUBROUTINE profil_sst(nlon, rlatd, type_profil, phy_sst)
     USE dimphy
+    USE phys_cal_mod , ONLY: year_len
     IMPLICIT NONE
 
     INTEGER nlon, type_profil, i, k, j
-    REAL :: rlatd(nlon), phy_sst(nlon, 360)
+    REAL :: rlatd(nlon), phy_sst(nlon, year_len)
     INTEGER imn, imx, amn, amx, kmn, kmx
     INTEGER p, pplus, nlat_max
@@ -825 +836 @@
     ENDIF
     WRITE (*, *) ' profil_sst: type_profil=', type_profil
-    DO i = 1, 360
+    DO i = 1, year_len
       ! phy_sst(:,i) = 260.+50.*cos(rlatd(:))**2
 
@@ -1018 +1029 @@
     imx = 1
     kmx = 1
-    DO k = 1, 360
+    DO k = 1, year_len
       DO i = 2, nlon
         IF (phy_sst(i,k)<amn) THEN