Changeset 3465 for LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.F90

Timestamp:

Mar 14, 2019, 10:34:31 AM (5 years ago)

Author:

Laurent Fairhead

Message:

Further modifications for DYNAMICO/LMDZ convergence. These are based
on Yann's LMDZ6_V2 sources. Compiles on irene and converges with revision 3459
in a bucket configuration
YM/LF

File:

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

LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.F90

@@ -218 +218 @@
       END DO
 
-    !--- Ensure first input longitudes interval contains first output point boundslon_reg(1,west)
-    dx1=locate(v1,boundslon_reg(1,west))-1
-    v1=CSHIFT(v1,SHIFT=dx1,DIM=1); v1(nlon_in-dx1+2:)=v1(nlon_in-dx1+2:)+2.*pi
-    !--- Extend input longitudes vector until last interval contains boundslon_reg(nlon_ou,east)
-    dx2=0; DO WHILE(v1(1+dx2)+2.*pi<=boundslon_reg(nlon_ou,east)); dx2=dx2+1; END DO
-
+      !--- Prepare quantities for time interpolation
+      tmidmonth=mid_month(annee_ref, cal_in)
+      IF(interpt) THEN
+        ntim_ou=ioget_year_len(annee_ref)
+        ALLOCATE(tmidday(ntim_ou))
+        tmidday=[(REAL(k)-0.5,k=1,ntim_ou)]
+        CALL ioget_calendar(cal_ou)
+      ELSE
+        ntim_ou=14
+        cal_ou=cal_in
+      END IF
+    ENDIF
+
+    IF (grid_type==unstructured) THEN
+      CALL bcast_mpi(nlon_in)
+      CALL bcast_mpi(nlat_in)
+      CALL bcast_mpi(nlev_in)
+      CALL bcast_mpi(l3d)
+      CALL bcast_mpi(tmidmonth)
+      CALL bcast_mpi(tmidday)
+      CALL bcast_mpi(ntim_ou)
+
+#ifdef CPP_XIOS    
+      IF (is_mpi_root) THEN
+        CALL xios_set_domain_attr("domain_climoz",nj_glo=nlat_in, nj=nlat_in, jbegin=0, latvalue_1d=lat_in/deg2rad)
+        IF (l3D) THEN
+          CALL xios_set_domain_attr("domain_climoz",ni_glo=nlon_in, ni=nlon_in, ibegin=0, lonvalue_1d=lon_in/deg2rad)
+        ELSE
+          CALL xios_set_domain_attr("domain_climoz",ni_glo=8, ni=8, ibegin=0, lonvalue_1d = (/ 0.,45.,90.,135.,180.,225.,270., 315. /))
+        ENDIF
+      ELSE
+        CALL xios_set_domain_attr("domain_climoz",nj_glo=nlat_in, nj=0, jbegin=0, latvalue_1d=null_array )
+        IF (l3D) THEN
+          CALL xios_set_domain_attr("domain_climoz",ni_glo=nlon_in, ni=0, ibegin=0, lonvalue_1d=null_array)
+        ELSE
+          CALL xios_set_domain_attr("domain_climoz",ni_glo=8, ni=0, ibegin=0, lonvalue_1d=null_array)
+        ENDIF
+      ENDIF
+      CALL  xios_set_axis_attr("axis_climoz", n_glo=nlev_in)
+      CALL  xios_set_axis_attr("time_axis_climoz", n_glo=ntim_ou)
+      CALL  xios_set_axis_attr("time_axis_climoz", n_glo=ntim_ou)
+      CALL  xios_set_axis_attr("tr_climoz", n_glo=read_climoz)
+      CALL  xios_set_field_attr("tro3_out", enabled=.TRUE.)
+      CALL  xios_set_field_attr("tro3_out", enabled=.TRUE.)
+#endif
+      
       IF (first) THEN
         first=.FALSE.
@@ -254 +294 @@
           dx1=locate(v1,boundslon_reg(1,west))-1
           v1=CSHIFT(v1,SHIFT=dx1,DIM=1) 
-          v1(nlon_in-dx1+1:)=v1(nlon_in-dx1+1:)+2.*pi
+          v1(nlon_in-dx1+2:)=v1(nlon_in-dx1+2:)+2.*pi
     
           !--- Extend input longitudes vector until last interval contains boundslon_reg(nlat_ou,east)
-          dx2=0; DO WHILE(v1(1+dx2)+2.*pi<boundslon_reg(nlon_ou,east)); dx2=dx2+1; END DO
+          dx2=0; DO WHILE(v1(1+dx2)+2.*pi<=boundslon_reg(nlon_ou,east)); dx2=dx2+1; END DO
 
           !--- Final edges longitudes vector (with margin and end point)