- Timestamp:
- 2019-12-05T10:57:46+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfcav.F90
r11987 r12062 25 25 USE par_oce , ONLY: jpi,jpj ! ocean space and time domain 26 26 USE phycst , ONLY: grav,rau0,r1_rau0_rcp ! physical constants 27 USE eosbn2 , ONLY: l _useCT ! l_useCT27 USE eosbn2 , ONLY: ln_teos10 ! use ln_teos10 or not 28 28 ! 29 29 USE in_out_manager ! I/O manager … … 84 84 ! 85 85 ! initialisation 86 IF (TRIM(cn_gammablk) == ' HJ99' ) zqoce_b (:,:) = ptsc(:,:,jp_tem) * rau0_rcp ! last time step total heat fluxes (to speed up convergence)86 IF (TRIM(cn_gammablk) == 'vel_stab' ) zqoce_b (:,:) = ptsc(:,:,jp_tem) * rau0_rcp ! last time step total heat fluxes (to speed up convergence) 87 87 ! 88 88 ! compute ice shelf melting … … 103 103 ! define if we need to iterate 104 104 SELECT CASE ( cn_gammablk ) 105 CASE ( 'spe',' ad15' )105 CASE ( 'spe','vel' ) 106 106 ! no convergence needed 107 107 lit = .FALSE. 108 CASE ( ' hj99' )108 CASE ( 'vel_stab' ) 109 109 ! compute error between 2 iterations 110 110 zerr = MAXVAL(ABS(zqoce(:,:) - zqoce_b(:,:))) … … 112 112 ! define if iteration needed 113 113 IF (nit >= 100) THEN ! too much iteration 114 CALL ctl_stop( 'STOP', 'isf_cav: HJ99gamma formulation had too many iterations ...' )114 CALL ctl_stop( 'STOP', 'isf_cav: vel_stab gamma formulation had too many iterations ...' ) 115 115 ELSE IF ( zerr <= 0.01_wp ) THEN ! convergence is achieve 116 116 lit = .FALSE. … … 163 163 INTEGER :: ierr 164 164 !!--------------------------------------------------------------------- 165 PRINT *, cn_isfcav_mlt166 165 ! 167 166 !============== … … 212 211 ! coeficient for linearisation of potential tfreez 213 212 ! Crude approximation for pressure (but commonly used) 214 IF ( l _useCT) THEN ! linearisation from Jourdain et al. (2017)213 IF ( ln_teos10 ) THEN ! linearisation from Jourdain et al. (2017) 215 214 risf_lamb1 =-0.0564_wp 216 215 risf_lamb2 = 0.0773_wp
Note: See TracChangeset
for help on using the changeset viewer.