Changeset 2839 for LMDZ5/branches/testing/libf/phylmd/open_climoz_m.F90

Timestamp:

Mar 30, 2017, 4:16:38 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r2785:2838 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/open_climoz_m.F90 (modified) (1 diff)

LMDZ5/branches/testing/libf/phylmd/open_climoz_m.F90

@@ -1 +1 @@
 ! $Id$
-module open_climoz_m
+MODULE open_climoz_m
 
-  implicit none
+  USE print_control_mod, ONLY: lunout
+  IMPLICIT NONE
 
-contains
+CONTAINS
 
-  subroutine open_climoz(ncid, press_in_edg)
+!-------------------------------------------------------------------------------
+!
+SUBROUTINE open_climoz(ncid, press_in_cen, press_in_edg, time_in, daily, adjust)
+!
+!-------------------------------------------------------------------------------
+  USE netcdf95, ONLY: nf95_open, nf95_close, nf95_gw_var, nf95_inq_varid
+  USE netcdf,   ONLY: nf90_nowrite
+  USE mod_phys_lmdz_mpi_data,      ONLY: is_mpi_root
+  USE mod_phys_lmdz_mpi_transfert, ONLY: bcast_mpi
+  USE phys_cal_mod,                ONLY: calend, year_len, year_cur
+!-------------------------------------------------------------------------------
+! Purpose: This procedure should be called once per "gcm" run, by a single
+!          thread of each MPI process.
+!          The root MPI process opens "climoz_LMDZ.nc", reads the pressure
+! levels and the times and broadcasts them to the other processes.
+!          We assume that, in "climoz_LMDZ.nc", the pressure levels are in hPa
+! and in strictly ascending order.
+!-------------------------------------------------------------------------------
+! Arguments (OpenMP shared):
+  INTEGER, INTENT(OUT):: ncid      !--- "climoz_LMDZ.nc" identifier
+  REAL, POINTER :: press_in_cen(:) !--- at cells centers
+  REAL, POINTER :: press_in_edg(:) !--- at the interfaces (pressure intervals)
+  REAL, POINTER :: time_in(:)      !--- records times, in days since Jan. 1st
+  LOGICAL, INTENT(IN) :: daily     !--- daily files (calendar dependent days nb)
+  LOGICAL, INTENT(IN) :: adjust    !--- tropopause adjustement required
+!   pressure levels press_in_cen/edg are in Pa a,d strictly ascending order.
+!   time_in is only used for monthly files (14 records)
+!-------------------------------------------------------------------------------
+! Local variables:
+  INTEGER :: varid                 !--- NetCDF variables identifier
+  INTEGER :: nlev, ntim            !--- pressure levels and time records numbers
+  CHARACTER(LEN=80)  :: sub
+  CHARACTER(LEN=320) :: msg
+!-------------------------------------------------------------------------------
+  sub="open_climoz"
+  WRITE(lunout,*)"Entering routine "//TRIM(sub)
 
-    ! This procedure should be called once per "gcm" run, by a single
-    ! thread of each MPI process.
-    ! The root MPI process opens "climoz_LMDZ.nc", reads the pressure
-    ! levels and broadcasts them to the other processes.
+  IF(is_mpi_root) THEN
 
-    ! We assume that, in "climoz_LMDZ.nc", the pressure levels are in hPa
-    ! and in strictly ascending order.
+    !--- OPEN FILE, READ PRESSURE LEVELS AND TIME VECTOR
+    CALL nf95_open("climoz_LMDZ.nc", nf90_nowrite, ncid)
+    CALL nf95_inq_varid(ncid, "plev", varid)
+    CALL nf95_gw_var(ncid, varid, press_in_cen)
+    ! Convert from hPa to Pa because "paprs" and "pplay" are in Pa:
+    press_in_cen = press_in_cen * 100.
+    nlev = SIZE(press_in_cen)
+    CALL NF95_INQ_VARID(ncID, "time", varID)
+    CALL NF95_GW_VAR(ncid, varid, time_in)
+    ntim = SIZE(time_in)
 
-    use netcdf95, only: nf95_open, nf95_close, nf95_gw_var, nf95_inq_varid
-    use netcdf, only: nf90_nowrite
+    !--- BUILD EDGES OF PRESSURE INTERVALS: HALFWAY IN LOGARITHMS
+    ALLOCATE(press_in_edg(nlev+1))
+    press_in_edg=[0.,SQRT(press_in_cen(1:nlev-1)*press_in_cen(2:nlev)),HUGE(0.)]
 
-    use mod_phys_lmdz_mpi_data, only: is_mpi_root
-    use mod_phys_lmdz_mpi_transfert, only: bcast_mpi ! broadcast
+    !--- CHECK RECORDS NUMBER AND DISPLAY CORRESPONDING INFORMATION
+    IF(daily.AND.ntim/=year_len) THEN
+      WRITE(msg,'(a,3(i4,a))')TRIM(sub)//': Expecting a daily ozone file with',&
+     &year_len,' records (year ',year_cur,') ; found ',ntim,' instead'
+      CALL abort_physic(sub, msg, 1)
+    ELSE IF(ALL([360,14]/=ntim)) THEN
+      WRITE(msg,'(a,i4,a)')TRIM(sub)//': Expecting an ozone file with 14 (mont'&
+     &//'hly case) or 360 (old style files) records ; found ',ntim,' instead'
+      CALL abort_physic(sub, msg, 1)
+    ELSE
+      IF(daily) THEN
+         WRITE(msg,'(a,2(i4,a))')'daily file (',ntim,' days in ',year_cur,')'
+       ELSE IF(ntim==14) THEN
+         msg='14 records monthly file'
+       ELSE
+         msg='360 records files (old convention)'
+       END IF
+      WRITE(lunout,*)TRIM(sub)//': Using a '//TRIM(msg)
+    END IF
 
-    integer, intent(out):: ncid ! of "climoz_LMDZ.nc", OpenMP shared
+    !--- MESSAGE ABOUT OPTIONAL STRETCHING FOR TROPOPAUSES MATCHING
+    IF(adjust) THEN
+      WRITE(lunout,*)TRIM(sub)//': Adjusting O3 field to match gcm tropopause.'
+    ELSE
+      WRITE(lunout,*)TRIM(sub)//': Interpolating O3 field directly on gcm levels.'
+    END IF
 
-    real, pointer:: press_in_edg(:)
-    ! edges of pressure intervals for ozone climatology, in Pa, in strictly
-    ! ascending order, OpenMP shared
+  END IF
+  CALL bcast_mpi(nlev)
+  IF(.NOT.is_mpi_root) ALLOCATE(press_in_cen(nlev  )); CALL bcast_mpi(press_in_cen)
+  IF(.NOT.is_mpi_root) ALLOCATE(press_in_edg(nlev+1)); CALL bcast_mpi(press_in_edg)
+  CALL bcast_mpi(ntim)
+  IF(.NOT.is_mpi_root) ALLOCATE(time_in(ntim));        CALL bcast_mpi(time_in)
 
-    ! Variables local to the procedure:
+END SUBROUTINE open_climoz
+!
+!-------------------------------------------------------------------------------
 
-    real, pointer:: plev(:)
-    ! (pressure levels for ozone climatology, converted to Pa, in strictly
-    ! ascending order)
-
-    integer varid ! for NetCDF
-    integer n_plev ! number of pressure levels in the input data
-    integer k
-
-    !---------------------------------------
-
-    print *, "Call sequence information: open_climoz"
-
-    if (is_mpi_root) then
-       call nf95_open("climoz_LMDZ.nc", nf90_nowrite, ncid)
-
-       call nf95_inq_varid(ncid, "plev", varid)
-       call nf95_gw_var(ncid, varid, plev)
-       ! Convert from hPa to Pa because "paprs" and "pplay" are in Pa:
-       plev = plev * 100.
-       n_plev = size(plev)
-    end if
-
-    call bcast_mpi(n_plev)
-    if (.not. is_mpi_root) allocate(plev(n_plev))
-    call bcast_mpi(plev)
-    
-    ! Compute edges of pressure intervals:
-    allocate(press_in_edg(n_plev + 1))
-    if (is_mpi_root) then
-       press_in_edg(1) = 0.
-       ! We choose edges halfway in logarithm:
-       forall (k = 2:n_plev) press_in_edg(k) = sqrt(plev(k - 1) * plev(k))
-       press_in_edg(n_plev + 1) = huge(0.)
-       ! (infinity, but any value guaranteed to be greater than the
-       ! surface pressure would do)
-    end if
-    call bcast_mpi(press_in_edg)
-    deallocate(plev) ! pointer
-
-  end subroutine open_climoz
-
-end module open_climoz_m
+END MODULE open_climoz_m