44 |
|
|
45 |
#ifdef ALLOW_SHAP_FILT |
#ifdef ALLOW_SHAP_FILT |
46 |
|
|
47 |
|
LOGICAL DIFFERENT_MULTIPLE |
48 |
|
EXTERNAL DIFFERENT_MULTIPLE |
49 |
|
|
50 |
C !LOCAL VARIABLES: |
C !LOCAL VARIABLES: |
|
#ifdef USE_OLD_SHAPIRO_FILTERS |
|
51 |
C == Local variables == |
C == Local variables == |
52 |
|
#ifdef USE_OLD_SHAPIRO_FILTERS |
53 |
C bi,bj,k :: loop index |
C bi,bj,k :: loop index |
54 |
INTEGER bi, bj, k |
INTEGER bi, bj, k |
55 |
#endif /* USE_OLD_SHAPIRO_FILTERS */ |
#endif /* USE_OLD_SHAPIRO_FILTERS */ |
56 |
|
CHARACTER*(MAX_LEN_MBUF) suff |
57 |
|
|
58 |
CEOP |
CEOP |
59 |
|
|
99 |
ELSE |
ELSE |
100 |
STOP 'SHAP_FILT_APPLY_UV: Ooops! Bad Shap_funct in UV block' |
STOP 'SHAP_FILT_APPLY_UV: Ooops! Bad Shap_funct in UV block' |
101 |
ENDIF |
ENDIF |
102 |
|
|
103 |
|
C----- Diagnostic of Shapiro Filter effect on Momentum : |
104 |
|
C Note: Shap_tmpFld1,2 from shap_filt_tracer_s2 (and not s1, s4) |
105 |
|
C are directly proportional to Delta-U,V due to the Filter |
106 |
|
IF ( Shap_funct.NE.1 .AND. Shap_funct.NE.4 .AND. |
107 |
|
& DIFFERENT_MULTIPLE(Shap_diagFreq,myTime,myTime-deltaTClock) |
108 |
|
& ) THEN |
109 |
|
_BARRIER |
110 |
|
_BEGIN_MASTER( myThid ) |
111 |
|
WRITE(suff,'(I10.10)') myIter |
112 |
|
CALL WRITE_FLD_XYZ_RL( 'shap_dU.', suff, Shap_tmpFld1, |
113 |
|
& myIter,myThid) |
114 |
|
CALL WRITE_FLD_XYZ_RL( 'shap_dV.', suff, Shap_tmpFld2, |
115 |
|
& myIter,myThid) |
116 |
|
_END_MASTER( myThid ) |
117 |
|
_BARRIER |
118 |
|
ENDIF |
119 |
|
|
120 |
ENDIF |
ENDIF |
121 |
|
|
122 |
#endif /* USE_OLD_SHAPIRO_FILTERS */ |
#endif /* USE_OLD_SHAPIRO_FILTERS */ |