Changeset 9780 for NEMO/trunk/src/NST


Ignore:
Timestamp:
2018-06-12T14:12:51+02:00 (3 years ago)
Author:
jchanut
Message:

Reorganize Agrif update in a single subroutine ; use adjoint stepping for initial state update

Location:
NEMO/trunk/src/NST
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/NST/agrif_ice_update.F90

    r9761 r9780  
    4242CONTAINS 
    4343 
    44    SUBROUTINE agrif_update_ice( kt ) 
     44   SUBROUTINE agrif_update_ice( ) 
    4545      !!---------------------------------------------------------------------- 
    4646      !!                     *** ROUTINE agrif_update_ice *** 
     
    4949      !! ** Action : - Update (u_ice,v_ice) and ice tracers 
    5050      !!---------------------------------------------------------------------- 
    51       INTEGER, INTENT(in) :: kt 
    52       !!---------------------------------------------------------------------- 
    5351      ! 
    5452      IF( Agrif_Root() .OR. nn_ice == 0 ) RETURN   ! do not update if inside Parent Grid or if child domain does not have ice 
    5553      ! 
    56       IF( ( MOD( (kt-nit000)/nn_fsbc + 1, Agrif_irhot() * Agrif_Parent(nn_fsbc)/nn_fsbc ) /=0 ) .AND. (kt /= 0) ) RETURN   ! update only at the parent ice time step 
     54      IF (lwp.AND.lk_agrif_debug) Write(*,*) 'Update sea ice from grid Number',Agrif_Fixed() 
     55      ! 
     56!      IF( ( MOD( (kt-nit000)/nn_fsbc + 1, Agrif_irhot() * Agrif_Parent(nn_fsbc)/nn_fsbc ) /=0 ) .AND. (kt /= 0) ) RETURN   ! update only at the parent ice time step 
    5757      ! 
    5858      Agrif_SpecialValueFineGrid    = -9999. 
  • NEMO/trunk/src/NST/agrif_oce_update.F90

    r9758 r9780  
    144144      Agrif_UseSpecialValueInUpdate = .FALSE. 
    145145      ! 
    146 #  if defined DECAL_FEEDBACK && defined VOL_REFLUX 
     146#  if defined VOL_REFLUX 
    147147      IF ( ln_dynspg_ts.AND.ln_bt_fw ) THEN 
    148148         ! Refluxing on ssh: 
     149#  if defined DECAL_FEEDBACK 
    149150         CALL Agrif_Update_Variable(ub2b_update_id,locupdate1=(/0, 0/),locupdate2=(/1, 1/),procname = reflux_sshu) 
    150151         CALL Agrif_Update_Variable(vb2b_update_id,locupdate1=(/1, 1/),locupdate2=(/0, 0/),procname = reflux_sshv) 
     152#  else 
     153         CALL Agrif_Update_Variable(ub2b_update_id,locupdate1=(/-1,-1/),locupdate2=(/ 0, 0/),procname = reflux_sshu) 
     154         CALL Agrif_Update_Variable(vb2b_update_id,locupdate1=(/ 0, 0/),locupdate2=(/-1,-1/),procname = reflux_sshv) 
     155#  endif 
    151156      END IF 
    152157#  endif 
     
    157162 
    158163 
    159    SUBROUTINE Agrif_Update_Tke( kt ) 
     164   SUBROUTINE Agrif_Update_Tke( ) 
    160165      !!--------------------------------------------- 
    161166      !!   *** ROUTINE Agrif_Update_Tke *** 
    162167      !!--------------------------------------------- 
    163168      !! 
    164       INTEGER, INTENT(in) :: kt  
    165169      !  
    166170      IF (Agrif_Root()) RETURN 
    167171      !        
    168       IF( ( Agrif_NbStepint() .NE. 0 ) .AND. (kt /= 0) ) RETURN 
    169172#  if defined TWO_WAY 
    170173 
  • NEMO/trunk/src/NST/agrif_user.F90

    r9761 r9780  
    6464   CALL Agrif_InitValues_cont_ice 
    6565# endif 
    66    ! 
    67    IF ( Agrif_Level().EQ.Agrif_MaxLevel() ) CALL agrif_Update_ini() 
    68  
    69    Agrif_UseSpecialValueInUpdate = .FALSE.      
    70  
     66   !     
    7167END SUBROUTINE Agrif_initvalues 
    7268 
     
    296292END SUBROUTINE Agrif_InitValues_cont 
    297293 
    298 RECURSIVE SUBROUTINE Agrif_Update_ini( ) 
    299    !!---------------------------------------------------------------------- 
    300    !!                 *** ROUTINE agrif_Update_ini *** 
    301    !! 
    302    !! ** Purpose :: Recursive update done at initialization 
    303    !!---------------------------------------------------------------------- 
    304    USE dom_oce 
    305    USE agrif_oce_update 
    306 #if defined key_top 
    307    USE agrif_top_update 
    308 #endif 
    309 #if defined key_si3 
    310    USE agrif_ice_update 
    311 #endif 
    312    ! 
    313    IMPLICIT NONE 
    314    !!---------------------------------------------------------------------- 
    315    ! 
    316    IF (Agrif_Root()) RETURN 
    317    ! 
    318    CALL Agrif_Update_ssh() 
    319    IF (.NOT.ln_linssh) CALL Agrif_Update_vvl() 
    320    CALL Agrif_Update_tra() 
    321 #if defined key_top 
    322    CALL Agrif_Update_Trc() 
    323 #endif 
    324    CALL Agrif_Update_dyn() 
    325 ! JC remove update because this precludes from perfect restartability 
    326 !!   CALL Agrif_Update_tke(0) 
    327  
    328 #if defined key_si3 
    329    CALL agrif_update_ice(0) 
    330 #endif 
    331     
    332    CALL Agrif_ChildGrid_To_ParentGrid() 
    333    CALL Agrif_Update_ini() 
    334    CALL Agrif_ParentGrid_To_ChildGrid() 
    335  
    336 END SUBROUTINE agrif_update_ini 
    337  
    338294SUBROUTINE agrif_declare_var 
    339295   !!---------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.