Changeset 3867 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/physiq.F90

Timestamp:

Jan 8, 2016, 9:52:14 AM (9 years ago)

Author:

ymipsl

Message:

• Switch to XIOS 2
• Append main part of forcing configuration (LMDZ stand alone)
• Create etat0 and limit from usual LMDZ input files, using XIOS 2 interpolation functionnalities
  • missing parametrization of gravity waves
  • missing aerosol

File:

• dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/physiq.F90 (modified) (10 diffs)

dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/physiq.F90

@@ -50 +50 @@
                               set_timestep, pdtphys
   USE tracinca_mod, ONLY : config_inca                   
-  USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev, grid1dTo2d_glo, klon_glo
+  USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, nbp_lev, grid1dTo2d_glo, klon_glo, grid_type, unstructured
 #ifdef REPROBUS
   USE CHEM_REP, ONLY : Init_chem_rep_xjour
@@ -64 +64 @@
   USE CFMIP_point_locations
   use FLOTT_GWD_rando_m, only: FLOTT_GWD_rando
-
+  USE etat0_limit_unstruct_mod
+  USE xios
+  USE regr_lat_time_climoz_m
+  USE climoz_mod
+  USE limit_read_mod, ONLY : init_limit_read
   IMPLICIT none
   !>======================================================================
@@ -247 +251 @@
   real clw(klon,klev),elij(klon,klev,klev)
   real epmlmMm(klon,klev,klev),eplaMm(klon,klev)
+  real delta_epsilon
   ! RomP <<<
   !IM definition dynamique o_trac dans phys_output_open
@@ -798 +803 @@
   !     ascending order
 
-  integer, save:: co3i = 0
+  integer, save:: co3i = -1
   !     time index in NetCDF file of current ozone fields
   !$OMP THREADPRIVATE(co3i) 
@@ -876 +881 @@
   mydebug=.FALSE.
   CALL set_timestep(pdtphys_)
-
+  IF (.NOT. debut)  CALL xios_update_calendar(itap+1)
   !======================================================================
   ! Ecriture eventuelle d'un profil verticale en entree de la physique.
@@ -904 +909 @@
 
   if (first) then 
-
+     CALL init_etat0_limit_unstruct
+     CALL init_limit_read(days_elapsed)
+     
      !CR:nvelles variables convection/poches froides
 
@@ -1055 +1062 @@
      radpas = NINT( 86400./pdtphys/nbapp_rad)
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
+     CALL init_iophy_new(lat_degrees,lon_degrees)
+
+     !=============================================================
+     !   Initialisation des sorties
+     !=============================================================
+
+#ifdef CPP_IOIPSL
+
+     !$OMP MASTER
+! FH : if ok_sync=.true. , the time axis is written at each time step
+! in the output files. Only at the end in the opposite case
+     ok_sync_omp=.false.
+     CALL getin('ok_sync',ok_sync_omp)
+     call phys_output_open(lon_degrees,lat_degrees,nCFMIP,tabijGCM, &
+          iGCM,jGCM,lonGCM,latGCM, &
+          jjmp1,nlevSTD,clevSTD,rlevSTD, pdtphys,ok_veget, &
+          type_ocean,iflag_pbl,iflag_pbl_split,ok_mensuel,ok_journe, &
+          ok_hf,ok_instan,ok_LES,ok_ade,ok_aie,  &
+          read_climoz, phys_out_filestations, &
+          new_aod, aerosol_couple, &
+          flag_aerosol_strat, pdtphys, paprs, pphis,  &
+          pplay, lmax_th, ptconv, ptconvth, ivap,  &
+          d_t, qx, d_qx, zmasse, ok_sync_omp)
+     !$OMP END MASTER
+     !$OMP BARRIER
+     ok_sync=ok_sync_omp
+
+     freq_outNMC(1) = ecrit_files(7)
+     freq_outNMC(2) = ecrit_files(8)
+     freq_outNMC(3) = ecrit_files(9)
+     WRITE(lunout,*)'OK freq_outNMC(1)=',freq_outNMC(1)
+     WRITE(lunout,*)'OK freq_outNMC(2)=',freq_outNMC(2)
+     WRITE(lunout,*)'OK freq_outNMC(3)=',freq_outNMC(3)
+
+     include "ini_histday_seri.h"
+
+     include "ini_paramLMDZ_phy.h"
+
+#endif
+     ecrit_reg = ecrit_reg * un_jour
+     ecrit_tra = ecrit_tra * un_jour
+
+     !XXXPB Positionner date0 pour initialisation de ORCHIDEE
+     date0 = jD_ref 
+     WRITE(*,*) 'physiq date0 : ',date0
+     
+!--- OZONE CLIMATOLOGY
+!    create ozone file for unstuctured grid
+     CALL create_climoz(read_climoz)
+
+     CALL phys_output_write(itap, pdtphys, paprs, pphis,                    &
+                             pplay, lmax_th, aerosol_couple,                 &
+                             ok_ade, ok_aie, ivap, new_aod, ok_sync,         &
+                             ptconv, read_climoz, clevSTD,                   &
+                             ptconvth, d_t, qx, d_qx, zmasse,                &
+                             flag_aerosol, flag_aerosol_strat, ok_cdnc)
+
+     CALL xios_update_calendar(1)    
+
+     CALL create_etat0_limit_unstruct
+     
+    
      CALL phyetat0 ("startphy.nc",clesphy0,tabcntr0)
      IF (klon_glo==1) THEN
@@ -1230 +1298 @@
      !IM cf. AM 081204 END
      !
-     !=============================================================
-     !   Initialisation des sorties
-     !=============================================================
-
-#ifdef CPP_IOIPSL
-
-     !$OMP MASTER
-! FH : if ok_sync=.true. , the time axis is written at each time step
-! in the output files. Only at the end in the opposite case
-     ok_sync_omp=.false.
-     CALL getin('ok_sync',ok_sync_omp)
-     call phys_output_open(lon_degrees,lat_degrees,nCFMIP,tabijGCM, &
-          iGCM,jGCM,lonGCM,latGCM, &
-          jjmp1,nlevSTD,clevSTD,rlevSTD, pdtphys,ok_veget, &
-          type_ocean,iflag_pbl,iflag_pbl_split,ok_mensuel,ok_journe, &
-          ok_hf,ok_instan,ok_LES,ok_ade,ok_aie,  &
-          read_climoz, phys_out_filestations, &
-          new_aod, aerosol_couple, &
-          flag_aerosol_strat, pdtphys, paprs, pphis,  &
-          pplay, lmax_th, ptconv, ptconvth, ivap,  &
-          d_t, qx, d_qx, zmasse, ok_sync_omp)
-     !$OMP END MASTER
-     !$OMP BARRIER
-     ok_sync=ok_sync_omp
-
-     freq_outNMC(1) = ecrit_files(7)
-     freq_outNMC(2) = ecrit_files(8)
-     freq_outNMC(3) = ecrit_files(9)
-     WRITE(lunout,*)'OK freq_outNMC(1)=',freq_outNMC(1)
-     WRITE(lunout,*)'OK freq_outNMC(2)=',freq_outNMC(2)
-     WRITE(lunout,*)'OK freq_outNMC(3)=',freq_outNMC(3)
-
-     include "ini_histday_seri.h"
-
-     include "ini_paramLMDZ_phy.h"
-
-#endif
-     ecrit_reg = ecrit_reg * un_jour
-     ecrit_tra = ecrit_tra * un_jour
-
-     !XXXPB Positionner date0 pour initialisation de ORCHIDEE
-     date0 = jD_ref 
-     WRITE(*,*) 'physiq date0 : ',date0
+
      !
      !
@@ -1611 +1637 @@
      ! Ozone from a file
      ! Update required ozone index:
-     ro3i = int((days_elapsed + jh_cur - jh_1jan) / year_len * 360.) + 1
+!ym     ro3i = int((days_elapsed + jh_cur - jh_1jan) / year_len * 360.) + 1
+!     ro3i = int((days_elapsed + jh_cur - jh_1jan )) + 1
+     
+     delta_epsilon=jh_cur-pdtphys/un_jour
+     if (ABS(delta_epsilon)<1e-6) delta_epsilon=0
+     ro3i = int((days_elapsed+delta_epsilon - jh_1jan ))+1
+     
+             
      if (ro3i == 361) ro3i = 360
      ! (This should never occur, except perhaps because of roundup
@@ -1618 +1651 @@
         ! Update ozone field:
         if (read_climoz == 1) then
-           call regr_pr_av(ncid_climoz, (/"tro3"/), julien=ro3i, &
+
+          IF (grid_type==unstructured) THEN
+            CALL get_ozone_var((/"tro3_read"/),press_climoz,paprs,wo)
+          ELSE
+             CALL regr_pr_av(ncid_climoz, (/"tro3"/), julien=ro3i, &
                 press_in_edg=press_climoz, paprs=paprs, v3=wo)
+          ENDIF
+      
         else
            ! read_climoz == 2
-           call regr_pr_av(ncid_climoz, (/"tro3         ", "tro3_daylight"/), &
-                julien=ro3i, press_in_edg=press_climoz, paprs=paprs, v3=wo)
+
+          IF (grid_type==unstructured) THEN
+            CALL get_ozone_var((/"tro3_read","tro3_daylight_read"/),press_climoz,paprs,wo)
+          ELSE
+            CALL regr_pr_av(ncid_climoz, (/"tro3         ", "tro3_daylight"/), &
+                 julien=ro3i, press_in_edg=press_climoz, paprs=paprs, v3=wo)
+          ENDIF
+
         end if
         ! Convert from mole fraction of ozone to column density of ozone in a