Changeset 2517

Timestamp:

May 20, 2016, 11:57:25 AM (8 years ago)

Author:

oboucher

Message:

Adding runofflic diagnostic in outputs
This is the freshwater flux going into the coupler

Location:

LMDZ5/trunk/libf/phylmd

Files:

• fonte_neige_mod.F90 (modified) (4 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/fonte_neige_mod.F90

@@ -313 +313 @@
 !
   SUBROUTINE fonte_neige_get_vars(pctsrf, fqcalving_out, &
-       fqfonte_out, ffonte_out)
-
+       fqfonte_out, ffonte_out, run_off_lic_out)
 
 
@@ -330 +329 @@
     REAL, DIMENSION(klon), INTENT(OUT)      :: fqfonte_out
     REAL, DIMENSION(klon), INTENT(OUT)      :: ffonte_out
+    REAL, DIMENSION(klon), INTENT(OUT)      :: run_off_lic_out
 
     INTEGER   :: nisurf
@@ -343 +343 @@
        fqcalving_out(:) = fqcalving_out(:) + fqcalving_global(:,nisurf)*pctsrf(:,nisurf)
     ENDDO
+    run_off_lic_out(:)=run_off_lic(:)
 
   END SUBROUTINE fonte_neige_get_vars
@@ -349 +350 @@
 !
 END MODULE fonte_neige_mod
-
-
-

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -273 +273 @@
   TYPE(ctrl_out), SAVE :: o_fqfonte = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
     'fqfonte', 'Land ice melt', 'kg/m2/s', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_runofflic = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
+    'runofflic', 'Land ice melt to ocean', 'kg/m2/s', (/ ('', i=1, 9) /))
   TYPE(ctrl_out), SAVE :: o_taux = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11 /), &
     'taux', 'Zonal wind stress', 'Pa', (/ ('', i=1, 9) /))

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -50 +50 @@
          o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &
          o_bils_latent, o_bils_enthalp, o_sens, &
-         o_fder, o_ffonte, o_fqcalving, o_fqfonte, &
+         o_fder, o_ffonte, o_fqcalving, o_fqfonte, o_runofflic, &
          o_taux, o_tauy, o_snowsrf, o_qsnow, &
          o_snowhgt, o_toice, o_sissnow, o_runoff, &
@@ -196 +196 @@
          zu10m, zv10m, zq2m, zustar, zxqsurf, &
          rain_lsc, rain_num, snow_lsc, bils, sens, fder, &
-         zxffonte, zxfqcalving, zxfqfonte, fluxu, &
+         zxffonte, zxfqcalving, zxfqfonte, zxrunofflic, fluxu, &
          fluxv, zxsnow, qsnow, snowhgt, to_ice, &
          sissnow, runoff, albsol3_lic, evap_pot, &
@@ -547 +547 @@
        CALL histwrite_phy(o_fqcalving, zxfqcalving)
        CALL histwrite_phy(o_fqfonte, zxfqfonte)
+       CALL histwrite_phy(o_runofflic, zxrunofflic)
        IF (vars_defined) THEN
           zx_tmp_fi2d=0.

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -68 +68 @@
 
 #ifdef CPP_RRTM
-    USE YOERAD   , ONLY : NRADLP
+    USE YOERAD, ONLY : NRADLP
 #endif
     USE ioipsl_getin_p_mod, ONLY : getin_p
@@ -3430 +3430 @@
        !albedo SB <<<
 
-
-       if (mydebug) then
-          call writefield_phy('u_seri',u_seri,nbp_lev)
-          call writefield_phy('v_seri',v_seri,nbp_lev)
-          call writefield_phy('t_seri',t_seri,nbp_lev)
-          call writefield_phy('q_seri',q_seri,nbp_lev)
-       endif
-
        !
        !sonia : If Iflag_radia >=2, pertubation of some variables
@@ -3483 +3475 @@
              print *,' ->radlwsw, number 1 '
           ENDIF
+
           !
           CALL radlwsw &
@@ -4254 +4247 @@
     ! Get some variables from module fonte_neige_mod
     CALL fonte_neige_get_vars(pctsrf,  &
-         zxfqcalving, zxfqfonte, zxffonte)
-
-
+         zxfqcalving, zxfqfonte, zxffonte, zxrunofflic)