Changeset 14072 for NEMO/trunk/src/OCE/SBC/sbcapr.F90
- Timestamp:
- 2020-12-04T08:48:38+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/SBC/sbcapr.F90
r14053 r14072 6 6 !! History : 3.3 ! 2010-09 (J. Chanut, C. Bricaud, G. Madec) Original code 7 7 !!---------------------------------------------------------------------- 8 8 9 9 !!---------------------------------------------------------------------- 10 !! sbc_apr : read atmospheric pressure in netcdf files 10 !! sbc_apr : read atmospheric pressure in netcdf files 11 11 !!---------------------------------------------------------------------- 12 12 USE dom_oce ! ocean space and time domain … … 25 25 PUBLIC sbc_apr ! routine called in sbcmod 26 26 PUBLIC sbc_apr_init ! routine called in sbcmod 27 27 28 28 ! !!* 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 30 30 LOGICAL, PUBLIC :: ln_ref_apr !: ref. pressure: global mean Patm (F) or a constant (F) 31 31 REAL(wp) :: rn_pref ! reference atmospheric pressure [N/m2] … … 34 34 REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) :: ssh_ibb ! Inverse barometer before sea surface height [m] 35 35 REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) :: apr ! atmospheric pressure at kt [N/m2] 36 36 37 37 REAL(wp) :: tarea ! whole domain mean masked ocean surface 38 38 REAL(wp) :: r1_grau ! = 1.e0 / (grav * rho0) 39 39 40 40 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_apr ! structure of input fields (file informations, fields read) 41 41 … … 54 54 !! 55 55 !! ** Method : - Read namelist namsbc_apr 56 !! - Read Patm fields in netcdf files 56 !! - Read Patm fields in netcdf files 57 57 !! - Compute reference atmospheric pressure 58 58 !! - Compute inverse barometer ssh … … 60 60 !! ssh_ib : inverse barometer ssh at kt 61 61 !!--------------------------------------------------------------------- 62 INTEGER :: ierror ! local integer 62 INTEGER :: ierror ! local integer 63 63 INTEGER :: ios ! Local integer output status for namelist read 64 64 !! … … 103 103 IF(lwp) WRITE(numout,*) ' Inverse barometer added to OBC ssh data' 104 104 ENDIF 105 !jc: stop below should rather be a warning 105 !jc: stop below should rather be a warning 106 106 IF( ln_apr_obc .AND. .NOT.ln_apr_dyn ) & 107 107 CALL ctl_warn( 'sbc_apr: use inverse barometer ssh at open boundary ONLY requires ln_apr_dyn=T' ) … … 116 116 !! 117 117 !! ** Method : - Read namelist namsbc_apr 118 !! - Read Patm fields in netcdf files 118 !! - Read Patm fields in netcdf files 119 119 !! - Compute reference atmospheric pressure 120 120 !! - Compute inverse barometer ssh … … 148 148 ! ! ---------------------------------------- ! 149 149 ! !* 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 151 151 IF(lwp) WRITE(numout,*) 'sbc_apr: ssh_ibb read in the restart file' 152 152 CALL iom_get( numror, jpdom_auto, 'ssh_ibb', ssh_ibb ) ! before inv. barometer ssh … … 167 167 ! 168 168 END SUBROUTINE sbc_apr 169 169 170 170 !!====================================================================== 171 171 END MODULE sbcapr
Note: See TracChangeset
for help on using the changeset viewer.