- Timestamp:
- 2015-06-04T17:01:01+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5218_CNRS17_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90
r5343 r5352 17 17 USE sbcapr ! surface boundary condition: atmospheric pressure 18 18 USE eosbn2 ! equation of state and related derivatives 19 USE traqsr, ONLY: ln_qsr_ice,fraqsr_1lev 19 20 ! 20 21 USE in_out_manager ! I/O manager … … 59 60 REAL(wp), DIMENSION(jpi,jpj,jpts) :: zts 60 61 !!--------------------------------------------------------------------- 61 62 62 63 ! !* surface T-, U-, V- ocean level variables (T, S, depth, velocity) 63 64 DO jj = 1, jpj … … 81 82 ENDIF 82 83 ! 83 IF( lk_vvl ) fse3t_m(:,:) = fse3t_n(:,:,1) 84 IF( lk_vvl ) e3t_m(:,:) = fse3t_n(:,:,1) 85 ! 86 IF( ln_qsr_ice ) frq_m(:,:) = fraqsr_1lev(:,:) 84 87 ! 85 88 ELSE … … 101 104 ENDIF 102 105 ! 103 IF( lk_vvl ) fse3t_m(:,:) = zcoef * fse3t_n(:,:,1) 106 IF( lk_vvl ) e3t_m(:,:) = zcoef * fse3t_n(:,:,1) 107 ! 108 IF( ln_qsr_ice ) frq_m(:,:) = zcoef * fraqsr_1lev(:,:) 104 109 ! ! ---------------------------------------- ! 105 110 ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN ! Initialisation: New mean computation ! … … 110 115 sss_m(:,:) = 0.e0 111 116 ssh_m(:,:) = 0.e0 112 IF( lk_vvl ) fse3t_m(:,:) = 0.e0 117 IF( lk_vvl ) e3t_m(:,:) = 0.e0 118 IF( ln_qsr_ice ) frq_m(:,:) = 0.e0 113 119 ENDIF 114 120 ! ! ---------------------------------------- ! … … 126 132 ENDIF 127 133 ! 128 IF( lk_vvl ) fse3t_m(:,:) = fse3t_m(:,:) + fse3t_n(:,:,1) 134 IF( lk_vvl ) e3t_m(:,:) = fse3t_m(:,:) + fse3t_n(:,:,1) 135 ! 136 IF( ln_qsr_ice ) frq_m(:,:) = frq_m(:,:) + fraqsr_1lev(:,:) 129 137 130 138 ! ! ---------------------------------------- ! … … 137 145 ssv_m(:,:) = ssv_m(:,:) * zcoef ! 138 146 ssh_m(:,:) = ssh_m(:,:) * zcoef ! mean SSH [m] 139 IF( lk_vvl ) fse3t_m(:,:) = fse3t_m(:,:) * zcoef ! mean vertical scale factor [m] 147 IF( lk_vvl ) e3t_m(:,:) = fse3t_m(:,:) * zcoef ! mean vertical scale factor [m] 148 IF( ln_qsr_ice ) frq_m(:,:) = frq_m(:,:) * zcoef ! mean fraction of solar net radiation absorbed in the 1st T level [-] 140 149 ! 141 150 ENDIF … … 154 163 CALL iom_rstput( kt, nitrst, numrow, 'sss_m' , sss_m ) 155 164 CALL iom_rstput( kt, nitrst, numrow, 'ssh_m' , ssh_m ) 156 IF( lk_vvl ) THEN 157 CALL iom_rstput( kt, nitrst, numrow, 'fse3t_m' , fse3t_m(:,:) ) 158 END IF 165 IF( lk_vvl ) CALL iom_rstput( kt, nitrst, numrow, 'e3t_m' , e3t_m ) 166 IF( ln_qsr_ice ) CALL iom_rstput( kt, nitrst, numrow, 'frq_m' , frq_m ) 159 167 ! 160 168 ENDIF … … 168 176 CALL iom_put( 'sss_m', sss_m ) 169 177 CALL iom_put( 'ssh_m', ssh_m ) 170 IF( lk_vvl ) CALL iom_put( 'e3t_m', fse3t_m(:,:) ) 178 IF( lk_vvl ) CALL iom_put( 'e3t_m', e3t_m ) 179 IF( ln_qsr_ice ) CALL iom_put( 'frq_m', frq_m ) 171 180 ENDIF 172 181 ! … … 204 213 CALL iom_get( numror, jpdom_autoglo, 'sss_m' , sss_m ) ! " " salinity (T-point) 205 214 CALL iom_get( numror, jpdom_autoglo, 'ssh_m' , ssh_m ) ! " " height (T-point) 206 IF( lk_vvl ) CALL iom_get( numror, jpdom_autoglo, 'fse3t_m', fse3t_m(:,:) ) 215 IF( lk_vvl ) CALL iom_get( numror, jpdom_autoglo, 'e3t_m', e3t_m ) 216 ! fraction of solar net radiation absorbed in 1st T level 217 IF( iom_varid( numror, 'frq_m', ldstop = .FALSE. ) > 0 ) THEN 218 CALL iom_get( numror, jpdom_autoglo, 'frq_m' , frq_m ) 219 ELSE 220 frq_m(:,:) = 1._wp ! default definition 221 ENDIF 207 222 ! 208 223 IF( zf_sbc /= REAL( nn_fsbc, wp ) ) THEN ! nn_fsbc has changed between 2 runs … … 215 230 sss_m(:,:) = zcoef * sss_m(:,:) 216 231 ssh_m(:,:) = zcoef * ssh_m(:,:) 217 IF( lk_vvl ) fse3t_m(:,:) = zcoef * fse3t_m(:,:) 232 IF( lk_vvl ) e3t_m(:,:) = zcoef * fse3t_m(:,:) 233 frq_m(:,:) = zcoef * frq_m(:,:) ! bug: must not be done if ln_qsr_ice = .false. 218 234 ELSE 219 235 IF(lwp) WRITE(numout,*) '~~~~~~~ mean fields read in the ocean restart file' … … 232 248 sss_m(:,:) = tsn(:,:,1,jp_sal) 233 249 ssh_m(:,:) = sshn(:,:) 234 IF( lk_vvl ) fse3t_m(:,:) = fse3t_n(:,:,1) 250 IF( lk_vvl ) e3t_m(:,:) = fse3t_n(:,:,1) 251 frq_m(:,:) = 1._wp 235 252 ! 236 253 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
