Changeset 15749

Timestamp:

2022-03-10T17:40:56+01:00 (2 years ago)

Author:

jwhile

Message:

Fix for pcbias boundary.

File:

• branches/UKMO/dev_r5518_GO6_package_FOAMv14_ASMFix/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_GO6_package_FOAMv14_ASMFix/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

@@ -524 +524 @@
          ! Use the inertial ramp.
          lenscl_bias = ( rn_maxlat_bias - rn_minlat_bias )*2._wp
+         fbcoef(:,:) = 1._wp - exp( -( abs( gphit(:,:) ) - rn_minlat_bias ) &
+                           * ( abs( gphit(:,:) ) - rn_minlat_bias ) / lenscl_bias )  
          WHERE ( abs( gphit(:,:) ) <= rn_minlat_bias )         
-            fbcoef(:,:) = 0._wp          
-         ELSEWHERE ( abs( gphit(:,:) ) >= rn_maxlat_bias )         
+            fbcoef(:,:) = 0._wp
+         ENDWHERE         
+         WHERE ( abs( gphit(:,:) ) >= rn_maxlat_bias )         
             fbcoef(:,:) = 1._wp                    
-         ELSEWHERE        
-            fbcoef(:,:) = 1._wp - exp( -( abs( gphit(:,:) ) - rn_minlat_bias ) &
-                           * ( abs( gphit(:,:) ) - rn_minlat_bias ) / lenscl_bias )                         
          ENDWHERE 
+         
       ELSEIF ( nn_lat_ramp == 2 ) THEN   
          ! Use a linear ramp consist with the geostrophic velocity balance ramp in NEMOVAR
-     
+         
+         fbcoef(:,:) = 1._wp - ((rn_maxlat_bias - abs( gphit(:,:)))/(rn_maxlat_bias - rn_minlat_bias))
          WHERE ( abs( gphit(:,:) ) <= rn_minlat_bias )
             fbcoef(:,:) = 0._wp
-         ELSEWHERE ( abs( gphit(:,:) ) >= rn_maxlat_bias ) 
+         ENDWHERE
+         WHERE ( abs( gphit(:,:) ) >= rn_maxlat_bias ) 
             fbcoef(:,:) = 1._wp
-         ELSEWHERE
-            fbcoef(:,:) = 1._wp - ((rn_maxlat_bias - abs( gphit(:,:)))/(rn_maxlat_bias - rn_minlat_bias))
          ENDWHERE
+            
+         
       ELSE
          fbcoef(:,:) = 0.0_wp