Changeset 3471

Timestamp:

Mar 27, 2019, 2:35:38 PM (6 years ago)

Author:

yann meurdesoif

Message:

Unstructured case : adapt limit,nc to hadisst files with daily values instead of monthly for sst and sic.

YM

File:

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

LMDZ6/trunk/libf/phylmd/create_limit_unstruct.F90

@@ -6 +6 @@
 
 CONTAINS
+
 
   SUBROUTINE create_limit_unstruct
@@ -18 +19 @@
   IMPLICIT NONE
     INCLUDE "iniprint.h"
-    REAL,    DIMENSION(klon,lmdep)                 :: sic
-    REAL,    DIMENSION(klon,lmdep)                 :: sst
+    REAL,    DIMENSION(:,:),ALLOCATABLE            :: sic
+    REAL,    DIMENSION(:,:),ALLOCATABLE            :: sst
     REAL,    DIMENSION(klon,lmdep)                 :: rugos
     REAL,    DIMENSION(klon,lmdep)                 :: albedo
-    REAL,    DIMENSION(klon_mpi,lmdep)             :: sic_mpi
-    REAL,    DIMENSION(klon_mpi,lmdep)             :: sst_mpi
+    REAL,    DIMENSION(:,:),ALLOCATABLE            :: sic_mpi
+    REAL,    DIMENSION(:,:),ALLOCATABLE            :: sst_mpi
     REAL,    DIMENSION(klon_mpi,lmdep)             :: rugos_mpi
     REAL,    DIMENSION(klon_mpi,lmdep)             :: albedo_mpi
@@ -41 +42 @@
     INTEGER :: l,k
     INTEGER :: nbad
-    
+    INTEGER :: sic_time_axis_size 
+    INTEGER :: sst_time_axis_size
+    CHARACTER(LEN=99)                  :: mess            ! error message
+    
+     
     ndays=ioget_year_len(annee_ref)
+    
+    IF (is_omp_master) CALL xios_get_axis_attr("time_sic",n_glo=sic_time_axis_size)
+    CALL bcast_omp(sic_time_axis_size)
+    ALLOCATE(sic_mpi(klon_mpi,sic_time_axis_size))
+    ALLOCATE(sic(klon,sic_time_axis_size))
+    
+    
+    IF (is_omp_master) CALL xios_get_axis_attr("time_sst",n_glo=sst_time_axis_size)
+    CALL bcast_omp(sst_time_axis_size)
+    ALLOCATE(sst_mpi(klon_mpi,sst_time_axis_size))
+    ALLOCATE(sst(klon,sst_time_axis_size))
     
     IF (is_omp_master) THEN
@@ -62 +78 @@
     ALLOCATE(phy_bil(klon,ndays)); phy_bil=0.0
 
+
 ! sic
-    CALL time_interpolation(ndays,sic,'gregorian',sic_year)
+    IF (sic_time_axis_size==lmdep) THEN
+      CALL time_interpolation(ndays,sic,'gregorian',sic_year)
+    ELSE IF (sic_time_axis_size==ndays) THEN
+      sic_year=sic
+    ELSE
+      WRITE(mess,*) 'sic time axis is nor montly, nor daily. sic time interpolation ',&
+                    'is requiered but is not currently managed'
+      CALL abort_physic('create_limit_unstruct',TRIM(mess),1)
+    ENDIF
+    
     sic_year(:,:)=sic_year(:,:)/100.  ! convert percent to fraction
     WHERE(sic_year(:,:)>1.0) sic_year(:,:)=1.0    ! Some fractions have some time large negative values
@@ -69 +95 @@
     
 ! sst
-    CALL time_interpolation(ndays,sst,'gregorian',sst_year)
+    IF (sst_time_axis_size==lmdep) THEN
+      CALL time_interpolation(ndays,sst,'gregorian',sst_year)
+    ELSE IF (sst_time_axis_size==ndays) THEN
+      sst_year=sst
+    ELSE
+      WRITE(mess,*)'sic time axis is nor montly, nor daily. sic time interpolation ',&
+                   'is requiered but is not currently managed'
+      CALL abort_physic('create_limit_unstruct',TRIM(mess),1)
+    ENDIF
     WHERE(sst_year(:,:)<271.38) sst_year(:,:)=271.38
+
 
 ! rugos