Changeset 4042 for branches/2013/dev_r4028_CNRS_LIM3

Timestamp:

2013-09-25T15:40:39+02:00 (11 years ago)

Author:

clem

Message:

new LIM3

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90 (modified) (12 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

@@ -11 +11 @@
    !!            3.3  ! 2010-11  (G. Madec) ice-ocean stress always computed at each ocean time-step
    !!            3.4  ! 2011-01  (A Porter)  dynamical allocation
+   !!             -   ! 2012-10  (C. Rousset)  add lim_diahsb
    !!----------------------------------------------------------------------
 #if defined key_lim3
@@ -42 +43 @@
    USE limsbc          ! sea surface boundary condition
    USE limdia          ! Ice diagnostics
+   USE limdiahsb       ! Ice budget diagnostics
    USE limwri          ! Ice outputs
    USE limrst          ! Ice restarts
-   USE limupdate       ! update of global variables
+   USE limupdate1       ! update of global variables
+   USE limupdate2       ! update of global variables
    USE limvar          ! Ice variables switch
 
@@ -51 +54 @@
    USE lib_mpp         ! MPP library
    USE wrk_nemo        ! work arrays
+   USE timing          ! Timing
    USE iom             ! I/O manager library
    USE in_out_manager  ! I/O manager
    USE prtctl          ! Print control
+
+#if defined key_bdy 
+   USE bdyice_lim       ! unstructured open boundary data  (bdy_ice_lim routine)
+#endif
 
    IMPLICIT NONE
@@ -98 +106 @@
       !!----------------------------------------------------------------------
 
+      IF( nn_timing == 1 )  CALL timing_start('sbc_ice_lim')
+
       CALL wrk_alloc( jpi,jpj,jpl, zalb_ice_os, zalb_ice_cs )
 
@@ -108 +118 @@
          !
          IF( ln_nicep ) THEN      ! control print at a given point
-            jiindx = 44   ;   jjindx = 140
+            jiindx = 15   ;   jjindx = 46
             WRITE(numout,*) ' The debugging point is : jiindx : ',jiindx, ' jjindx : ',jjindx
          ENDIF
@@ -178 +188 @@
          d_oa_i_thd (:,:,:)   = 0._wp   ;   d_oa_i_trp (:,:,:)   = 0._wp
          !
-         sfx    (:,:) = 0._wp
+         sfx    (:,:) = 0._wp   ;   sfx_thd  (:,:) = 0._wp
          sfx_bri(:,:) = 0._wp   ;   sfx_mec  (:,:) = 0._wp   ;   sfx_res  (:,:) = 0._wp
          fhbri  (:,:) = 0._wp   ;   fheat_mec(:,:) = 0._wp   ;   fheat_res(:,:) = 0._wp
@@ -185 +195 @@
          focea2D(:,:) = 0._wp
          fsup2D (:,:) = 0._wp
-         ! 
+
+         ! used in limthd.F90
+         rdvosif(:,:) = 0._wp   ! variation of ice volume at surface
+         rdvobif(:,:) = 0._wp   ! variation of ice volume at bottom
+         fdvolif(:,:) = 0._wp   ! total variation of ice volume
+         rdvonif(:,:) = 0._wp   ! lateral variation of ice volume
+         fstric (:,:) = 0._wp   ! part of solar radiation transmitted through the ice
+         ffltbif(:,:) = 0._wp   ! linked with fstric
+         qfvbq  (:,:) = 0._wp   ! linked with fstric
+         rdm_snw(:,:) = 0._wp   ! variation of snow mass per unit area
+         rdm_ice(:,:) = 0._wp   ! variation of ice mass per unit area
+         hicifp (:,:) = 0._wp   ! daily thermodynamic ice production. 
+         !
          diag_sni_gr(:,:) = 0._wp   ;   diag_lat_gr(:,:) = 0._wp
          diag_bot_gr(:,:) = 0._wp   ;   diag_dyn_gr(:,:) = 0._wp
          diag_bot_me(:,:) = 0._wp   ;   diag_sur_me(:,:) = 0._wp
+         diag_res_pr(:,:) = 0._wp   ;   diag_trp_vi(:,:) = 0._wp
          ! dynamical invariants
          delta_i(:,:) = 0._wp       ;   divu_i(:,:) = 0._wp       ;   shear_i(:,:) = 0._wp
@@ -199 +222 @@
                           CALL lim_dyn( kt )              ! Ice dynamics    ( rheology/dynamics )
                           CALL lim_trp( kt )              ! Ice transport   ( Advection/diffusion )
-                          CALL lim_var_agg(1)             ! aggregate categories, requested
                           CALL lim_var_glo2eqv            ! equivalent variables, requested for rafting
          IF( ln_nicep )   CALL lim_prt_state( jiindx, jjindx,-1, ' - ice dyn & trp - ' )   ! control print
                           CALL lim_itd_me                 ! Mechanical redistribution ! (ridging/rafting)
+                          CALL lim_var_agg( 1 ) 
+                          CALL lim_update1
          ENDIF
+!                         !- Change old values for new values
+                          old_u_ice(:,:)   = u_ice (:,:)
+                          old_v_ice(:,:)   = v_ice (:,:)
+                          old_a_i(:,:,:)   = a_i (:,:,:)
+                          old_v_s(:,:,:)   = v_s (:,:,:)
+                          old_v_i(:,:,:)   = v_i (:,:,:)
+                          old_e_s(:,:,:,:) = e_s (:,:,:,:)
+                          old_e_i(:,:,:,:) = e_i (:,:,:,:)
+                          old_oa_i(:,:,:)  = oa_i(:,:,:)
+                          old_smv_i(:,:,:) = smv_i (:,:,:)
          !                                           ! Ice thermodynamics 
                           CALL lim_var_glo2eqv            ! equivalent variables
@@ -217 +251 @@
          !                                           ! Global variables update
                           CALL lim_var_agg( 1 )           ! requested by limupdate
-                          CALL lim_update                 ! Global variables update
+                          CALL lim_update2                 ! Global variables update
+#if defined key_bdy
+                          CALL bdy_ice_lim( kt )          ! clem modif: bdy ice
+#endif
                           CALL lim_var_glo2eqv            ! equivalent variables (outputs)
                           CALL lim_var_agg(2)             ! aggregate ice thickness categories
@@ -229 +266 @@
          IF( ( MOD( kt+nn_fsbc-1, ninfo ) == 0 .OR. ntmoy == 1 ) .AND. .NOT. lk_mpp )   &
             &             CALL lim_dia 
+         IF (ln_limdiahsb) CALL lim_diahsb
                           CALL lim_wri( 1  )              ! Ice outputs 
+         IF( kt == nit000 )   CALL iom_close( numrir )  ! clem: close input ice restart file
          IF( lrst_ice )   CALL lim_rst_write( kt )        ! Ice restart file 
                           CALL lim_var_glo2eqv            ! ???
@@ -248 +287 @@
       !
       CALL wrk_dealloc( jpi,jpj,jpl, zalb_ice_os, zalb_ice_cs )
+      !
+      IF( nn_timing == 1 )  CALL timing_stop('sbc_ice_lim')
       !
    END SUBROUTINE sbc_ice_lim
@@ -607 +648 @@
 !       WRITE(numout,*) ' sfx_bri    : ', sfx_bri  (ki,kj)
 !       WRITE(numout,*) ' sfx        : ', sfx      (ki,kj)
-!       WRITE(numout,*) ' fsalt_res  : ', fsalt_res(ki,kj)
+!       WRITE(numout,*) ' sfx_res  : ', sfx_res(ki,kj)
         WRITE(numout,*) ' fmmec      : ', fmmec    (ki,kj)
         WRITE(numout,*) ' fhmec      : ', fhmec    (ki,kj)