Changeset 11182 for NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_phy.F90

Timestamp:

2019-06-26T12:08:17+02:00 (5 years ago)

Author:

laurent

Message:

LB: adding possibility to use dew-point temperature rather than specific humidity as "air humidity" => ln_humi_dpt@namsbc_blk

File:

• NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_phy.F90 (moved) (moved from NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/src/OCE/SBC/sbcblk_phymbl.F90) (5 diffs)

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

@@ -1 @@
-MODULE sbcblk_phymbl
+MODULE sbcblk_phy
    !!======================================================================
-   !!                       ***  MODULE  sbcblk_phymbl  ***
+   !!                       ***  MODULE  sbcblk_phy  ***
    !! A set of functions to compute air themodynamics parameters
    !!                     needed by Aerodynamic Bulk Formulas
@@ -15 +15 @@
    !!   visc_air      : kinematic viscosity (aka Nu_air) of air from temperature
    
-   USE oce            ! ocean dynamics and tracers
    USE dom_oce        ! ocean space and time domain
    USE phycst         ! physical constants
-   !USE sbc_oce        ! Surface boundary condition: ocean fields
-   !USE sbcdcy         ! surface boundary condition: diurnal cycle
-   !USE sbcwave , ONLY :   cdn_wave ! wave module
-   !USE sbc_ice        ! Surface boundary condition: ice fields
-   USE lib_fortran    ! to use key_nosignedzero
-   !
-   !USE iom            ! I/O manager library
-   !USE in_out_manager ! I/O manager
-   USE lib_mpp        ! distribued memory computing library
-   !USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
-   !USE prtctl         ! Print control
 
    IMPLICIT NONE
@@ -98 +86 @@
       !
       INTEGER  ::   ji, jj         ! dummy loop indices
-      REAL(wp) ::   ze_sat, ztmp   ! local scalar
+      REAL(wp) ::   zta, ze_sat, ztmp   ! local scalar
       !!----------------------------------------------------------------------------------
       !
@@ -104 +92 @@
          DO ji = 1, jpi
             !
-            ztmp = rt0 / ptak(ji,jj)
+            zta = MAX( ptak(ji,jj) , 180._wp )   ! air temp., prevents fpe0 errors dute to unrealistically low values over masked regions...
+            ztmp = rt0 / zta
             !
             ! Vapour pressure at saturation [hPa] : WMO, (Goff, 1957)
-            ze_sat = 10.**( 10.79574*(1. - ztmp) - 5.028*LOG10(ptak(ji,jj)/rt0)        &
-               &    + 1.50475*10.**(-4)*(1. - 10.**(-8.2969*(ptak(ji,jj)/rt0 - 1.)) )  &
+            ze_sat = 10.**( 10.79574*(1. - ztmp) - 5.028*LOG10(zta/rt0)        &
+               &    + 1.50475*10.**(-4)*(1. - 10.**(-8.2969*(zta/rt0 - 1.)) )  &
                &    + 0.42873*10.**(-3)*(10.**(4.76955*(1. - ztmp)) - 1.) + 0.78614  )
             !
@@ -191 +180 @@
 
    !!======================================================================
-END MODULE sbcblk_phymbl
+END MODULE sbcblk_phy