Changeset 12029

Timestamp:

2019-12-02T17:24:47+01:00 (5 years ago)

Author:

laurent

Message:

Making "sbcblk.F90" compliant for AGRIF (function returns a scalar)

File:

• NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90 (modified) (1 diff)

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90

@@ -506 +506 @@
       !!    (see Josey, Gulev & Yu, 2013) / doi=10.1016/B978-0-12-391851-2.00005-2
       !!    (since reanalysis products provide absolute temperature "T" at z, not theta !)
-      ztpot = sf(jp_tair)%fnow(:,:,1) + gamma_moist( sf(jp_tair)%fnow(:,:,1), zqair(:,:) ) * rn_zqt
-
+      !!#LB: because AGRIF hates functions that return something else than a scalar, need to
+      !!     use vectorial version of gamma_moist() ...
+      DO jj = 2, jpjm1
+         DO ji = fs_2, fs_jpim1   ! vect. opt.
+            ztpot(ji,jj) = sf(jp_tair)%fnow(ji,jj,1) + gamma_moist( sf(jp_tair)%fnow(ji,jj,1), zqair(ji,jj) ) * rn_zqt
+         END DO
+      END DO
+      
 
       !! Time to call the user-selected bulk parameterization for