Changeset 632 for trunk/NEMO/OPA_SRC/DIA/diawri.F90

Timestamp:

2007-03-02T17:45:17+01:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_bugfix_026 : CT : -add key_flx_core to save in restart files the nfbulk parameter and the gsst(:,:) field when using CORE forcing

• save surface fluxes qla, qsb and qlw

File:

• trunk/NEMO/OPA_SRC/DIA/diawri.F90 (modified) (9 diffs)

trunk/NEMO/OPA_SRC/DIA/diawri.F90

@@ -120 +120 @@
       CHARACTER (len=80) :: clname
       !!----------------------------------------------------------------------
-      
+
       ! 0. Initialisation
       ! -----------------
-      
+
       ! local variable for debugging
       ll_print = .FALSE.
@@ -160 +160 @@
 
          ! Define the NETCDF files (one per grid)
-         
+
          ! Compute julian date from starting date of the run
          CALL ymds2ju( nyear, nmonth, nday, 0.e0, zjulian )
@@ -175 +175 @@
          WRITE(inum,*) clhstnam
          CLOSE(inum)
-         
+
          ! Define the T grid FILE ( nid_T )
-         
+
          CALL dia_nam( clhstnam, nwrite, 'grid_T' )
          IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam    ! filename
@@ -224 +224 @@
          CALL histvert( nid_W, "depthw", "Vertical W levels",      &  ! Vertical grid: gdepw
             &          "m", ipk, gdepw_0, nz_W )
- 
+
 
          ! Declare all the output fields as NETCDF variables
@@ -280 +280 @@
             &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout )
 #endif
+
+#if defined key_flx_core
+         CALL histdef( nid_T, "solhflup", "Latent Heat Flux Upward"         , "W/m2"   ,   &  ! qla
+            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout )
+         CALL histdef( nid_T, "solwfldo", "Longwave Radiation downward"     , "W/m2"   ,   &  ! qlw
+            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout )
+         CALL histdef( nid_T, "sosbhfup", "Sensible Heat Flux upward"       , "W/m2"   ,   &  ! qsb
+            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout )
+#endif
+
 
 #if ( defined key_coupled && ! defined key_ice_lim ) 
@@ -432 +442 @@
       zw2d(:,:) = erp(:,:) * sn(:,:,1) * tmask(:,:,1)
       CALL histwrite( nid_T, "sosafldp", it, zw2d          , ndim_hT, ndex_hT )   ! salt flux damping
+#endif
+#if defined key_flx_core
+      CALL histwrite( nid_T, "solhflup", it, qla           , ndim_hT, ndex_hT )   ! latent heat flux
+      CALL histwrite( nid_T, "solwfldo", it, qlw           , ndim_hT, ndex_hT )   ! longwave heat flux
+      CALL histwrite( nid_T, "sosbhfup", it, qsb           , ndim_hT, ndex_hT )   ! sensible heat flux
 #endif
 #if ( defined key_coupled && ! defined key_ice_lim ) 
@@ -556 +571 @@
       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~   and forcing fields file created '
       IF(lwp) WRITE(numout,*) '                and named :', cdfile_name, '.nc'
-      
+
       ! 0. Initialisation
       ! -----------------
-      
+
       ! Define frequency of output and means
       zdt  = rdt
@@ -617 +632 @@
       ! donne le nombre d'elements, et idex la liste des indices a sortir
       idex(1) = 1   ! init to avoid compil warning
-      
+
       ! Write all fields on T grid
       CALL histwrite( id_i, "votemper", 1, tn    , jpi*jpj*jpk, idex )    ! now temperature
@@ -641 +656 @@
 
    END SUBROUTINE dia_wri_state
- 
+
 #endif
    !!======================================================================