Changeset 628 for trunk/NEMO/NST_SRC

Timestamp:

2007-02-28T12:22:51+01:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_update_006:RB: first integration of AGRIF for passive tracers

Location:

trunk/NEMO/NST_SRC

Files:

• agrif_top_interp.F90 (added)
• agrif_top_update.F90 (added)
• agrif_user.F90 (modified) (8 diffs)

trunk/NEMO/NST_SRC/agrif_user.F90

@@ -59 +59 @@
       USE ice_oce
 #endif
+#if defined key_passivetrc
+     USE agrif_top_update
+     USE agrif_top_interp
+     USE sms
+#endif
 #if defined key_agrif
      USE agrif_opa_update
@@ -72 +77 @@
 !
       REAL(wp) tabtemp(jpi,jpj,jpk)
+#if defined key_passivetrc
+      REAL(wp) tabtrtemp(jpi,jpj,jpk,jptra)
+#endif
 ! 
       LOGICAL check_namelist
@@ -115 +123 @@
       Call Agrif_Set_type(gcb,(/2,2/),(/3,3/))
 
+#if defined key_passivetrc
+      Call Agrif_Set_type(trb,(/2,2,0,0/),(/3,3,0,0/))
+      Call Agrif_Set_type(trn,(/2,2,0,0/),(/3,3,0,0/))
+      Call Agrif_Set_type(tra,(/2,2,0,0/),(/3,3,0,0/))
+#endif
+
+
+
 !
 !     Space directions for each variables
@@ -139 +155 @@
       Call Agrif_Set_raf(gcb,(/'x','y'/))
 
+#if defined key_passivetrc
+      Call Agrif_Set_raf(trb,(/'x','y','N','N'/))
+      Call Agrif_Set_raf(trn,(/'x','y','N','N'/))
+      Call Agrif_Set_raf(tra,(/'x','y','N','N'/))
+#endif
+
 !
 !     type of interpolation
@@ -157 +179 @@
       Call Agrif_Set_bcinterp(e2v,interp1=AGRIF_ppm,interp2=Agrif_linear)
 
+#if defined key_passivetrc
+      Call Agrif_Set_bcinterp(trn,interp=AGRIF_linear)
+      Call Agrif_Set_bcinterp(tra,interp=AGRIF_linear)
+#endif
+
 !
 !     Location of interpolation
@@ -175 +202 @@
       Call Agrif_Set_bc(va,(/-2*Agrif_irhox(),0/))
 
+#if defined key_passivetrc
+      Call Agrif_Set_bc(trn,(/0,1/))
+      Call Agrif_Set_bc(tra,(/-3*Agrif_irhox(),0/))
+#endif
+
 !     Update type
       
@@ -185 +217 @@
       Call Agrif_Set_Updatetype(sshn, update = AGRIF_Update_Average)
       Call Agrif_Set_Updatetype(gcb,update = AGRIF_Update_Average)
+
+#if defined key_passivetrc
+      Call Agrif_Set_Updatetype(trn, update = AGRIF_Update_Average)
+      Call Agrif_Set_Updatetype(trb, update = AGRIF_Update_Average)
+#endif
 
       Call Agrif_Set_Updatetype(un,update1 = Agrif_Update_Copy, update2 = Agrif_Update_Average)
@@ -206 +243 @@
        Call Agrif_Bc_variable(tabtemp,ua,calledweight=1.,procname=interpun)
        Call Agrif_Bc_variable(tabtemp,va,calledweight=1.,procname=interpvn)
+
+#if defined key_passivetrc
+       Call Agrif_Bc_variable(tabtrtemp,trn,calledweight=1.)
+!       Call Agrif_Bc_variable(tabtrtemp,tra,calledweight=1.,procname=interptrn)
+       Call Agrif_Bc_variable(tabtrtemp,tra,calledweight=1.)
+
+#endif
        Agrif_UseSpecialValue = .FALSE.
+
 !