Changeset 12384 for NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/TRA/tramle.F90

Timestamp:

2020-02-13T17:21:21+01:00 (4 years ago)

Author:

agn

Message:

corrections to tramle & zdfosm to avoid division by 0

File:

• NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/TRA/tramle.F90 (modified) (1 diff)

NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/TRA/tramle.F90

@@ -254 +254 @@
     DO jj = 1, jpjm1
        DO ji = 1, fs_jpim1   ! vector opt.
-          zhu(ji,jj) = 1._wp / zhu(ji,jj)                   ! hu --> 1/hu
-          zhv(ji,jj) = 1._wp / zhv(ji,jj)
+          zhu(ji,jj) = 1._wp / MAX(zhu(ji,jj), rsmall)                   ! hu --> 1/hu
+          zhv(ji,jj) = 1._wp / MAX(zhv(ji,jj), rsmall) 
        END DO
     END DO