Changeset 1786 for branches/dev_005_AWL/NEMO/NST_SRC/agrif_user.F90

Timestamp:

2009-12-02T09:50:10+01:00 (14 years ago)

Author:

sga

Message:

add new AGRIF without LIM code to NEMO branch dev_005_AWL taken from NOCS NEMO branch noc_dev_024_AWL revision 1051 (patch file differences from NOCS NEMO trunk revision 1043)

File:

• branches/dev_005_AWL/NEMO/NST_SRC/agrif_user.F90 (modified) (7 diffs)

branches/dev_005_AWL/NEMO/NST_SRC/agrif_user.F90

@@ -88 +88 @@
       REAL(wp) :: tabtrtemp(jpi,jpj,jpk,jptra)
 #endif 
+#if defined key_agrif_nolim
+      REAL(wp) :: tabflxtemp(jpi,jpj)
+#endif
       LOGICAL check_namelist
       !!----------------------------------------------------------------------
@@ -135 +138 @@
 #endif
       
+#if defined key_lim2 && defined key_agrif_nolim
+      Call Agrif_Set_type(utau_nst,(/1,2/),(/2,3/))
+      Call Agrif_Set_type(vtau_nst,(/2,1/),(/3,2/))
+      Call Agrif_Set_type(taum_nst,(/2,2/),(/3,3/))
+      Call Agrif_Set_type(qsr_nst,(/2,2/),(/3,3/))
+      Call Agrif_Set_type(qns_nst,(/2,2/),(/3,3/))
+      Call Agrif_Set_type(emp_nst,(/2,2/),(/3,3/))
+      Call Agrif_Set_type(emps_nst,(/2,2/),(/3,3/))
+      Call Agrif_Set_type(wndm_nst,(/2,2/),(/3,3/))
+      Call Agrif_Set_type(fri_nst,(/2,2/),(/3,3/))
+      Call Agrif_Set_type(tag_nst,(/2,2/),(/3,3/))
+#endif
+
       ! 2. Space directions for each variables
       !---------------------------------------
@@ -164 +180 @@
 #endif
 
+#if defined key_lim2 && defined key_agrif_nolim
+      Call Agrif_Set_raf(utau_nst,(/'x','y'/))
+      Call Agrif_Set_raf(vtau_nst,(/'x','y'/))
+      Call Agrif_Set_raf(taum_nst,(/'x','y'/))
+      Call Agrif_Set_raf(qsr_nst,(/'x','y'/))
+      Call Agrif_Set_raf(qns_nst,(/'x','y'/))
+      Call Agrif_Set_raf(emp_nst,(/'x','y'/))
+      Call Agrif_Set_raf(emps_nst,(/'x','y'/))
+      Call Agrif_Set_raf(wndm_nst,(/'x','y'/))
+      Call Agrif_Set_raf(fri_nst,(/'x','y'/))
+      Call Agrif_Set_raf(tag_nst,(/'x','y'/))
+#endif
+
       ! 3. Type of interpolation
       !------------------------- 
@@ -186 +215 @@
 #endif
 
+#if defined key_lim2 && defined key_agrif_nolim
+      Call Agrif_Set_interp(utau_nst,interp1=Agrif_linear,interp2=AGRIF_ppm)
+      Call Agrif_Set_interp(vtau_nst,interp1=AGRIF_ppm,interp2=Agrif_linear)
+      Call Agrif_Set_interp(taum_nst,interp=AGRIF_linear)
+      Call Agrif_Set_interp(qsr_nst,interp=AGRIF_linear)
+      Call Agrif_Set_interp(qns_nst,interp=AGRIF_linear)
+      Call Agrif_Set_interp(emp_nst,interp=AGRIF_linear)
+      Call Agrif_Set_interp(emps_nst,interp=AGRIF_linear)
+      Call Agrif_Set_interp(wndm_nst,interp=AGRIF_linear)
+      Call Agrif_Set_interp(fri_nst,interp=AGRIF_linear)
+      Call Agrif_Set_interp(tag_nst,interp=AGRIF_linear)
+
+      Call Agrif_Set_bcinterp(utau_nst,interp1=Agrif_linear,interp2=AGRIF_ppm)
+      Call Agrif_Set_bcinterp(vtau_nst,interp1=AGRIF_ppm,interp2=Agrif_linear)
+      Call Agrif_Set_bcinterp(taum_nst,interp=AGRIF_linear)
+      Call Agrif_Set_bcinterp(qsr_nst,interp=AGRIF_linear)
+      Call Agrif_Set_bcinterp(qns_nst,interp=AGRIF_linear)
+      Call Agrif_Set_bcinterp(emp_nst,interp=AGRIF_linear)
+      Call Agrif_Set_bcinterp(emps_nst,interp=AGRIF_linear)
+      Call Agrif_Set_bcinterp(wndm_nst,interp=AGRIF_linear)
+      Call Agrif_Set_bcinterp(fri_nst,interp=AGRIF_linear)
+      Call Agrif_Set_bcinterp(tag_nst,interp=AGRIF_linear)
+#endif
+
       ! 4. Location of interpolation
       !-----------------------------
@@ -208 +261 @@
 #endif
 
+#if defined key_lim2 && defined key_agrif_nolim
+      Call Agrif_Set_bc(utau_nst,(/0,1/))
+      Call Agrif_Set_bc(vtau_nst,(/0,1/))
+      Call Agrif_Set_bc(taum_nst,(/0,1/))
+      Call Agrif_Set_bc(qsr_nst,(/0,1/))
+      Call Agrif_Set_bc(qns_nst,(/0,1/))
+      Call Agrif_Set_bc(emp_nst,(/0,1/))
+      Call Agrif_Set_bc(emps_nst,(/0,1/))
+      Call Agrif_Set_bc(wndm_nst,(/0,1/))
+      Call Agrif_Set_bc(fri_nst,(/0,1/))
+      Call Agrif_Set_bc(tag_nst,(/0,1/))
+#endif
+
       ! 5. Update type
       !--------------- 
@@ -229 +295 @@
       Call Agrif_Set_Updatetype(e1u,update1 = Agrif_Update_Copy, update2=Agrif_Update_Average)
       Call Agrif_Set_Updatetype(e2v,update1 = Agrif_Update_Average, update2=Agrif_Update_Copy)
+
+#if defined key_lim2 && defined key_agrif_nolim
+      Call Agrif_Set_Updatetype(utau_nst,update1 = Agrif_Update_Copy, update2 = Agrif_Update_Average)
+      Call Agrif_Set_Updatetype(vtau_nst,update1 = Agrif_Update_Average, update2 = Agrif_Update_Copy)
+      Call Agrif_Set_Updatetype(taum_nst , update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(qsr_nst , update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(qns_nst , update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(emp_nst , update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(emps_nst, update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(wndm_nst, update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(fri_nst, update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(tag_nst , update = AGRIF_Update_Average)
+#endif
 
       ! 6. First interpolations of potentially non zero fields
@@ -248 +327 @@
       Call Agrif_Bc_variable(tabtrtemp,trn,calledweight=1.)
       Call Agrif_Bc_variable(tabtrtemp,tra,calledweight=1.,procname=interptrn)
+#endif
+
+#if defined key_lim2 && defined key_agrif_nolim
+
+!     Call Agrif_Interp_variable(tabflxtemp,utau_nst,procname=interputau)
+!     Call Agrif_Interp_variable(tabflxtemp,vtau_nst,procname=interpvtau)
+      Call Agrif_Interp_variable(tabflxtemp,utau_nst)
+      Call Agrif_Interp_variable(tabflxtemp,vtau_nst)
+      Call Agrif_Interp_variable(tabflxtemp,taum_nst)
+      Call Agrif_Interp_variable(tabflxtemp,qsr_nst)
+      Call Agrif_Interp_variable(tabflxtemp,qns_nst)
+      Call Agrif_Interp_variable(tabflxtemp,emp_nst)
+      Call Agrif_Interp_variable(tabflxtemp,emps_nst)
+      Call Agrif_Interp_variable(tabflxtemp,wndm_nst)
+      Call Agrif_Interp_variable(tabflxtemp,fri_nst)
+      Call Agrif_Interp_variable(tabflxtemp,tag_nst)
+!     Call Agrif_Bc_variable(tabflxtemp,utau_nst,calledweight=1.,procname=interputau)
+!     Call Agrif_Bc_variable(tabflxtemp,vtau_nst,calledweight=1.,procname=interpvtau)
+      Call Agrif_Bc_variable(tabflxtemp,utau_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,vtau_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,taum_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,qsr_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,qns_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,emp_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,emps_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,wndm_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,fri_nst,calledweight=1.)
+      Call Agrif_Bc_variable(tabflxtemp,tag_nst,calledweight=1.)
 #endif
       Agrif_UseSpecialValue = .FALSE.