Changeset 2804 for branches/2011/dev_r2802_LOCEAN10_agrif_lim/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

Timestamp:

2011-07-18T17:05:52+02:00 (13 years ago)

Author:

rblod

Message:

dev_r2802_LOCEAN10_agrif_lim: first implementation see ticket #848

File:

• branches/2011/dev_r2802_LOCEAN10_agrif_lim/NEMOGCM/NEMO/NST_SRC/agrif_user.F90 (modified) (4 diffs)

branches/2011/dev_r2802_LOCEAN10_agrif_lim/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

@@ -84 +84 @@
 # if ! defined key_offline
       CALL Agrif_InitValues_cont
-# endif       
+#  if defined key_lim2
+!      CALL Agrif_InitValues_cont_lim2
+#  endif       
+# endif
 # if defined key_top
       CALL Agrif_InitValues_cont_top
@@ -279 +282 @@
 
    END SUBROUTINE agrif_declare_var
+
+
+#  if defined key_lim2
+   SUBROUTINE Agrif_InitValues_cont_lim2
+      !!----------------------------------------------------------------------
+      !!                 *** ROUTINE Agrif_InitValues_cont_lim2 ***
+      !!
+      !! ** Purpose :: Initialisation of variables to be interpolated for LIM2
+      !!----------------------------------------------------------------------
+      USE Agrif_Util
+      USE ice_2
+      USE in_out_manager
+      USE agrif_lim2_update
+      USE agrif_lim2_interp
+      USE lib_mpp 
+      !
+      IMPLICIT NONE
+      !
+      REAL(wp), DIMENSION(:,:)  , ALLOCATABLE :: zvel
+      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: zadv, zsadv
+      !!----------------------------------------------------------------------
+
+      ALLOCATE( zvel(jpi,jpj), zadv(jpi,jpj,7), zsadv(jpi,jpj,42) )
+      
+      ! 1. Declaration of the type of variable which have to be interpolated
+      !---------------------------------------------------------------------
+      CALL agrif_declare_var_lim2
+
+      ! 2. First interpolations of potentially non zero fields
+      !-------------------------------------------------------
+      Agrif_SpecialValue=-9999.
+      Agrif_UseSpecialValue = .TRUE.
+      Call Agrif_Bc_variable(zadv ,adv_ice_id ,calledweight=1.,procname=interp_adv_ice )    
+      Call Agrif_Bc_variable(zsadv,sadv_ice_id,calledweight=1.,procname=interp_sadv_ice)    
+      Call Agrif_Bc_variable(zvel ,u_ice_id   ,calledweight=1.,procname=interp_u_ice   )
+      Call Agrif_Bc_variable(zvel ,v_ice_id   ,calledweight=1.,procname=interp_v_ice   )
+      Agrif_SpecialValue=0.
+      Agrif_UseSpecialValue = .FALSE.
+
+      ! 3. Some controls
+      !-----------------
+       
+#   if ! defined key_lim2_vp
+      childfreq = 1.
+      CALL agrif_dyn_lim(0 ,1 , 'V')
+      CALL agrif_dyn_lim(0 ,1 , 'U')
+      childfreq = 0.
+#   endif
+!RB mandatory but why ???
+      IF( nbclineupdate /= nn_fsbc .AND. nn_ice == 2 )THEN
+         CALL ctl_warn ('With ice model on child grid, nbclineupdate is set to nn_fsbc')
+         nbclineupdate = nn_fsbc
+       ENDIF
+      CALL Agrif_Update_lim2(0)
+      !
+      DEALLOCATE( zvel, zadv, zsadv )
+      !
+   END SUBROUTINE Agrif_InitValues_cont_lim2
+
+
+   SUBROUTINE agrif_declare_var_lim2
+      !!----------------------------------------------------------------------
+      !!                 *** ROUTINE agrif_declare_var_lim2 ***
+      !!
+      !! ** Purpose :: Declaration of variables to be interpolated for LIM2
+      !!----------------------------------------------------------------------
+      USE agrif_util
+      USE ice_2
+
+      IMPLICIT NONE
+      !!----------------------------------------------------------------------
+   
+      ! 1. Declaration of the type of variable which have to be interpolated
+      !---------------------------------------------------------------------
+      CALL agrif_declare_variable((/2,2,0/),(/3,3,0/),(/'x','y','N'/),(/1,1,1/),(/jpi,jpj, 7/),adv_ice_id )         
+      CALL agrif_declare_variable((/2,2,0/),(/3,3,0/),(/'x','y','N'/),(/1,1,1/),(/jpi,jpj,42/),sadv_ice_id)         
+#   if defined key_lim2_vp
+      CALL agrif_declare_variable((/1,1/),(/3,3/),(/'x','y'/),(/1,1/),(/jpi,jpj/),u_ice_id)
+      CALL agrif_declare_variable((/1,1/),(/3,3/),(/'x','y'/),(/1,1/),(/jpi,jpj/),v_ice_id)
+#   else
+      CALL agrif_declare_variable((/1,2/),(/2,3/),(/'x','y'/),(/1,1/),(/jpi,jpj/),u_ice_id)
+      CALL agrif_declare_variable((/2,1/),(/3,2/),(/'x','y'/),(/1,1/),(/jpi,jpj/),v_ice_id)
+#   endif
+          
+      ! 2. Type of interpolation
+      !-------------------------
+      CALL Agrif_Set_bcinterp(adv_ice_id ,interp=AGRIF_linear)   
+      CALL Agrif_Set_bcinterp(sadv_ice_id,interp=AGRIF_linear)   
+      Call Agrif_Set_bcinterp(u_ice_id,interp1=Agrif_linear,interp2=AGRIF_ppm)
+      Call Agrif_Set_bcinterp(v_ice_id,interp1=AGRIF_ppm,interp2=Agrif_linear)
+
+      ! 3. Location of interpolation
+      !-----------------------------
+      Call Agrif_Set_bc(adv_ice_id ,(/0,1/))
+      Call Agrif_Set_bc(sadv_ice_id,(/0,1/))
+      Call Agrif_Set_bc(u_ice_id,(/0,1/))
+      Call Agrif_Set_bc(v_ice_id,(/0,1/))
+
+      ! 5. Update type
+      !--------------- 
+      Call Agrif_Set_Updatetype(adv_ice_id , update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(sadv_ice_id, update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(u_ice_id,update1 = Agrif_Update_Copy, update2 = Agrif_Update_Average)
+      Call Agrif_Set_Updatetype(v_ice_id,update1 = Agrif_Update_Average, update2 = Agrif_Update_Copy)
+
+   END SUBROUTINE agrif_declare_var_lim2
+#  endif
 # endif
    
@@ -459 +569 @@
       !!----------------------------------------------------------------------
       USE agrif_oce 
+      USE agrif_ice 
       USE in_out_manager
       USE lib_mpp
@@ -486 +597 @@
       visc_dyn      = rn_sponge_dyn
       !
-      IF( agrif_oce_alloc()  > 0 )   CALL ctl_warn('agrif sol_oce_alloc: allocation of arrays failed')
+      IF( agrif_oce_alloc()  > 0 )   CALL ctl_stop('agrif agrif_oce_alloc: allocation of arrays failed')
+# if defined key_lim2
+      IF( agrif_ice_alloc()  > 0 )   CALL ctl_stop('agrif agrif_ice_alloc: allocation of arrays failed')
+# endif
       !
     END SUBROUTINE agrif_nemo_init