Changeset 12029


Ignore:
Timestamp:
2019-12-02T17:24:47+01:00 (7 months ago)
Author:
laurent
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk.F90

    r12028 r12029  
    506506      !!    (see Josey, Gulev & Yu, 2013) / doi=10.1016/B978-0-12-391851-2.00005-2 
    507507      !!    (since reanalysis products provide absolute temperature "T" at z, not theta !) 
    508       ztpot = sf(jp_tair)%fnow(:,:,1) + gamma_moist( sf(jp_tair)%fnow(:,:,1), zqair(:,:) ) * rn_zqt 
    509  
     508      !!#LB: because AGRIF hates functions that return something else than a scalar, need to 
     509      !!     use vectorial version of gamma_moist() ... 
     510      DO jj = 2, jpjm1 
     511         DO ji = fs_2, fs_jpim1   ! vect. opt. 
     512            ztpot(ji,jj) = sf(jp_tair)%fnow(ji,jj,1) + gamma_moist( sf(jp_tair)%fnow(ji,jj,1), zqair(ji,jj) ) * rn_zqt 
     513         END DO 
     514      END DO 
     515       
    510516 
    511517      !! Time to call the user-selected bulk parameterization for 
Note: See TracChangeset for help on using the changeset viewer.