Changeset 1482 for trunk/NEMO/LIM_SRC_2/limwri_2.F90

Timestamp:

2009-07-03T17:28:06+02:00 (15 years ago)

Author:

smasson

Message:

distribution of iom_put + cleaning of LIM2 outputs, see ticket:437

File:

• trunk/NEMO/LIM_SRC_2/limwri_2.F90 (modified) (7 diffs)

trunk/NEMO/LIM_SRC_2/limwri_2.F90

@@ -4 +4 @@
    !!         Ice diagnostics :  write ice output files
    !!======================================================================
-   !! history :  2.0  ! 03-08  (C. Ethe) original code
-   !!            2.0  ! 04-10  (C. Ethe )  1D configuration
+   !! history :  2.0  ! 2003-08  (C. Ethe)      original code
+   !!            2.0  ! 2004-10  (C. Ethe )     1D configuration
+   !!             -   ! 2009-06  (B. Lemaire )  iom_put + lim_wri_state_2
    !!-------------------------------------------------------------------
 #if defined key_lim2
@@ -14 +15 @@
    !!   lim_wri_2      : write of the diagnostics variables in ouput file 
    !!   lim_wri_init_2 : initialization and namelist read
+   !!   lim_wri_state_2 : write for initial state or/and abandon:
+   !!                     > output.init.nc (if ninist = 1 in namelist)
+   !!                     > output.abort.nc
    !!----------------------------------------------------------------------
    USE phycst
@@ -32 +36 @@
    PRIVATE
 
-   PUBLIC   lim_wri_2      ! routine called by sbc_ice_lim_2
+#if ! defined key_iomput
+   PUBLIC   lim_wri_2         ! called by sbc_ice_lim_2
+#endif
+   PUBLIC   lim_wri_state_2   ! called by dia_wri_state 
 
    INTEGER, PARAMETER                       ::   jpnoumax = 40   ! maximum number of variable for ice output
@@ -61 +68 @@
 CONTAINS
 
-#if defined key_dimgout
+#if ! defined key_iomput
+# if defined key_dimgout
    !!----------------------------------------------------------------------
    !!   'key_dimgout'                                    Direct Access file
    !!----------------------------------------------------------------------
 # include "limwri_dimg_2.h90"
-#else
-   !!----------------------------------------------------------------------
-   !!   Default option                                          NetCDF file
-   !!----------------------------------------------------------------------
-
+# else
    SUBROUTINE lim_wri_2( kt )
       !!-------------------------------------------------------------------
@@ -92 +96 @@
       REAL(wp), DIMENSION(jpi,jpj,jpnoumax) ::   zcmo
       !!-------------------------------------------------------------------
-
-      CALL iom_setkt( kt + nn_fsbc - 1 )
-      !                                          !--------------------!
+                                                 !--------------------!
       IF( kt == nit000 ) THEN                    !   Initialisation   !
          !                                       !--------------------!
@@ -180 +182 @@
          
          IF( nc(jf) == 1 )   CALL histwrite( nice, nam(jf), niter, zfield, ndim, ndex51 )
-         CALL iom_put( nam(jf), zfield )
          
       END DO
       
       IF( ( nn_fsbc * niter ) >= nitend )   CALL histclo( nice ) 
-      !
-      CALL iom_setkt( kt )
 
    END SUBROUTINE lim_wri_2
-   
-#endif
-   
+     
+
    SUBROUTINE lim_wri_init_2
       !!-------------------------------------------------------------------
@@ -273 +271 @@
    END SUBROUTINE lim_wri_init_2
 
+# endif
+#endif
+
+   SUBROUTINE lim_wri_state_2( kt, kid, kh_i )
+      !!---------------------------------------------------------------------
+      !!                 ***  ROUTINE lim_wri_state_2  ***
+      !!        
+      !! ** Purpose :   create a NetCDF file named cdfile_name which contains 
+      !!      the instantaneous ice state and forcing fields for ice model
+      !!        Used to find errors in the initial state or save the last
+      !!      ocean state in case of abnormal end of a simulation
+      !!
+      !! History :
+      !!   2.0  !  2009-06  (B. Lemaire)
+      !!----------------------------------------------------------------------
+      INTEGER, INTENT( in ) ::   kt               ! ocean time-step index)
+      INTEGER, INTENT( in ) ::   kid , kh_i       
+      !!----------------------------------------------------------------------
+
+      CALL histdef( kid, "isnowthi", "Snow thickness"          , "m"      , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      CALL histdef( kid, "iicethic", "Ice thickness"           , "m"      , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      CALL histdef( kid, "iiceprod", "Ice produced"            , "m/kt"   , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      CALL histdef( kid, "ileadfra", "Ice concentration"       , "-"      , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      CALL histdef( kid, "iicetemp", "Ice temperature"         , "K"      , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      CALL histdef( kid, "ioceflxb", "flux at ice base"        , "w/m2"   , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
+      CALL histdef( kid, "iicevelu", "i-Ice speed (I-point)"   , "m/s"    , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      CALL histdef( kid, "iicevelv", "j-Ice speed (I-point)"   , "m/s"    , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
+      CALL histdef( kid, "isstempe", "Sea surface temperature" , "C"      , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
+      CALL histdef( kid, "isssalin", "Sea surface salinity"    , "PSU"    , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
+      CALL histdef( kid, "iicestru", "i-Wind stress over ice (I-pt)", "Pa", jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      CALL histdef( kid, "iicestrv", "j-Wind stress over ice (I-pt)", "Pa", jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
+      CALL histdef( kid, "iicesflx", "Solar flux over ice"     , "w/m2"   , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
+      CALL histdef( kid, "iicenflx", "Non-solar flux over ice" , "w/m2"   , jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt )
+      CALL histdef( kid, "isnowpre", "Snow precipitation"      , "kg/m2/s", jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
+
+      CALL histend( kid )   ! end of the file definition
+
+      CALL histwrite( kid, "isnowthi", kt, hsnif          , jpi*jpj, (/1/) )   
+      CALL histwrite( kid, "iicethic", kt, hicif          , jpi*jpj, (/1/) )    
+      CALL histwrite( kid, "iiceprod", kt, hicifp         , jpi*jpj, (/1/) )   
+      CALL histwrite( kid, "ileadfra", kt, 1. - frld(:,:) , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "iicetemp", kt, sist(:,:) - rt0, jpi*jpj, (/1/) )
+      CALL histwrite( kid, "ioceflxb", kt, fbif           , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "iicevelv", kt, u_ice          , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "iicevelu", kt, v_ice          , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "isstempe", kt, sst_m          , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "isssalin", kt, sss_m          , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "iicestru", kt, utau_ice       , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "iicestrv", kt, vtau_ice       , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "iicesflx", kt, qsr_ice(:,:,1) , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "iicenflx", kt, qns_ice(:,:,1) , jpi*jpj, (/1/) )
+      CALL histwrite( kid, "isnowpre", kt, sprecip        , jpi*jpj, (/1/) )
+
+    END SUBROUTINE lim_wri_state_2
+
 #else
    !!----------------------------------------------------------------------