Changeset 3110 for LMDZ6/trunk/libf

Timestamp:

Dec 4, 2017, 5:19:33 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Added sza, solar zenithal angle, to outputs for CMIP6 dataresquest

Location:

LMDZ6/trunk/libf/phylmd

Files:

• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_var_mod.F90 (modified) (3 diffs)
• phys_output_write_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -51 +51 @@
 
 !!! 2D
+  TYPE(ctrl_out), SAVE :: o_sza = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11, 11/), &
+    'sza', 'Solar zenithal angle', 'degrees', (/ ('', i=1, 10) /))
 
 ! Marine

LMDZ6/trunk/libf/phylmd/phys_output_var_mod.F90

@@ -14 +14 @@
   REAL, SAVE, ALLOCATABLE :: snow_o(:), zfra_o(:)
   !$OMP THREADPRIVATE(snow_o, zfra_o)
+  REAL, SAVE, ALLOCATABLE :: sza_o(:) ! solar zenithal angle
+  !$OMP THREADPRIVATE(sza_o)
   INTEGER, SAVE, ALLOCATABLE ::  itau_con(:)       ! Nombre de pas ou rflag <= 1
   !$OMP THREADPRIVATE(itau_con)
@@ -144 +146 @@
 
     allocate(snow_o(klon), zfra_o(klon))
+    allocate(sza_o(klon) )
     allocate(itau_con(klon))
     allocate(sens_prec_liq_o(klon,2))
@@ -198 +201 @@
 
     deallocate(snow_o,zfra_o,itau_con)
+    deallocate(sza_o)
     deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
     deallocate (d_qw_col, d_ql_col, d_qs_col, d_qt_col, d_ek_col, d_h_dair_col &

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -178 +178 @@
          o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
          o_lat_prec_sol_oce, o_lat_prec_sol_sic, &
+         o_sza, &
 ! Marine
          o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &
@@ -319 +320 @@
          itau_con, nfiles, clef_files, nid_files, dryaod_diag, &
          zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
-         zustr_gwd_front, zvstr_gwd_front,     &
+         zustr_gwd_front, zvstr_gwd_front, sza_o,    &
          sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
          cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
@@ -614 +615 @@
      ENDIF
 
+       CALL histwrite_phy(o_sza, sza_o)
        CALL histwrite_phy(o_flat, zxfluxlat)
        CALL histwrite_phy(o_ptstar, ptstar)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1122 +1122 @@
     logical, parameter :: mass_fixer=.false.
     real qql1(klon),qql2(klon),corrqql
+
+    REAL pi
+
+    pi = 4. * ATAN(1.)
 
     ! Ehouarn: set value of jjmp1 since it is no longer a "fixed parameter"
@@ -2190 +2194 @@
        END SELECT
     ENDIF
+    sza_o = ACOS (rmu0) *180./pi
 
     IF (mydebug) THEN