New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 14072 for NEMO/trunk/src/OCE/SBC/sbcapr.F90 – NEMO

Ignore:
Timestamp:
2020-12-04T08:48:38+01:00 (3 years ago)
Author:
laurent
Message:

Merging branch "2020/dev_r13648_ASINTER-04_laurent_bulk_ice", ticket #2369

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/SBC/sbcapr.F90

    r14053 r14072  
    66   !! History :  3.3  !   2010-09  (J. Chanut, C. Bricaud, G. Madec)  Original code 
    77   !!---------------------------------------------------------------------- 
    8     
     8 
    99   !!---------------------------------------------------------------------- 
    10    !!   sbc_apr        : read atmospheric pressure in netcdf files  
     10   !!   sbc_apr        : read atmospheric pressure in netcdf files 
    1111   !!---------------------------------------------------------------------- 
    1212   USE dom_oce         ! ocean space and time domain 
     
    2525   PUBLIC   sbc_apr       ! routine called in sbcmod 
    2626   PUBLIC   sbc_apr_init  ! routine called in sbcmod 
    27     
     27 
    2828   !                                          !!* namsbc_apr namelist (Atmospheric PRessure) * 
    29    LOGICAL, PUBLIC ::   ln_apr_obc = .false.   !: inverse barometer added to OBC ssh data  
     29   LOGICAL, PUBLIC ::   ln_apr_obc = .false.   !: inverse barometer added to OBC ssh data 
    3030   LOGICAL, PUBLIC ::   ln_ref_apr             !: ref. pressure: global mean Patm (F) or a constant (F) 
    3131   REAL(wp)        ::   rn_pref                !  reference atmospheric pressure   [N/m2] 
     
    3434   REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) ::   ssh_ibb   ! Inverse barometer before sea surface height   [m] 
    3535   REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) ::   apr       ! atmospheric pressure at kt                 [N/m2] 
    36     
     36 
    3737   REAL(wp) ::   tarea                ! whole domain mean masked ocean surface 
    3838   REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rho0) 
    39     
     39 
    4040   TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf_apr   ! structure of input fields (file informations, fields read) 
    4141 
     
    5454      !! 
    5555      !! ** Method  : - Read namelist namsbc_apr 
    56       !!              - Read Patm fields in netcdf files  
     56      !!              - Read Patm fields in netcdf files 
    5757      !!              - Compute reference atmospheric pressure 
    5858      !!              - Compute inverse barometer ssh 
     
    6060      !!                ssh_ib   : inverse barometer ssh at kt 
    6161      !!--------------------------------------------------------------------- 
    62       INTEGER            ::   ierror  ! local integer  
     62      INTEGER            ::   ierror  ! local integer 
    6363      INTEGER            ::   ios     ! Local integer output status for namelist read 
    6464      !! 
     
    103103         IF(lwp) WRITE(numout,*) '         Inverse barometer added to OBC ssh data' 
    104104      ENDIF 
    105 !jc: stop below should rather be a warning  
     105!jc: stop below should rather be a warning 
    106106      IF( ln_apr_obc .AND. .NOT.ln_apr_dyn   )   & 
    107107            CALL ctl_warn( 'sbc_apr: use inverse barometer ssh at open boundary ONLY requires ln_apr_dyn=T' ) 
     
    116116      !! 
    117117      !! ** Method  : - Read namelist namsbc_apr 
    118       !!              - Read Patm fields in netcdf files  
     118      !!              - Read Patm fields in netcdf files 
    119119      !!              - Compute reference atmospheric pressure 
    120120      !!              - Compute inverse barometer ssh 
     
    148148         !                                      ! ---------------------------------------- ! 
    149149         !                                            !* Restart: read in restart file 
    150          IF( ln_rstart .AND. .NOT.l_1st_euler ) THEN  
     150         IF( ln_rstart .AND. .NOT.l_1st_euler ) THEN 
    151151            IF(lwp) WRITE(numout,*) 'sbc_apr:   ssh_ibb read in the restart file' 
    152152            CALL iom_get( numror, jpdom_auto, 'ssh_ibb', ssh_ibb )   ! before inv. barometer ssh 
     
    167167      ! 
    168168   END SUBROUTINE sbc_apr 
    169        
     169 
    170170   !!====================================================================== 
    171171END MODULE sbcapr 
Note: See TracChangeset for help on using the changeset viewer.