- Timestamp:
- 2015-11-06T11:42:35+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90
r5845 r5866 78 78 ENDIF 79 79 ! 80 IF( lk_vvl) e3t_m(:,:) = e3t_n(:,:,1)80 IF( .NOT.ln_linssh ) e3t_m(:,:) = e3t_n(:,:,1) 81 81 ! 82 82 frq_m(:,:) = fraqsr_1lev(:,:) … … 100 100 ENDIF 101 101 ! 102 IF( lk_vvl) e3t_m(:,:) = zcoef * e3t_n(:,:,1)102 IF( .NOT.ln_linssh ) e3t_m(:,:) = zcoef * e3t_n(:,:,1) 103 103 ! 104 104 frq_m(:,:) = zcoef * fraqsr_1lev(:,:) … … 106 106 ELSEIF( MOD( kt - 2 , nn_fsbc ) == 0 ) THEN ! Initialisation: New mean computation ! 107 107 ! ! ---------------------------------------- ! 108 ssu_m(:,:) = 0. e0! reset to zero ocean mean sbc fields109 ssv_m(:,:) = 0. e0110 sst_m(:,:) = 0. e0111 sss_m(:,:) = 0. e0112 ssh_m(:,:) = 0. e0113 IF( lk_vvl ) e3t_m(:,:) = 0.e0114 frq_m(:,:) = 0. e0108 ssu_m(:,:) = 0._wp ! reset to zero ocean mean sbc fields 109 ssv_m(:,:) = 0._wp 110 sst_m(:,:) = 0._wp 111 sss_m(:,:) = 0._wp 112 ssh_m(:,:) = 0._wp 113 IF( .NOT.ln_linssh ) e3t_m(:,:) = 0._wp 114 frq_m(:,:) = 0._wp 115 115 ENDIF 116 116 ! ! ---------------------------------------- ! … … 128 128 ENDIF 129 129 ! 130 IF( lk_vvl) e3t_m(:,:) = e3t_m(:,:) + e3t_n(:,:,1)130 IF( .NOT.ln_linssh ) e3t_m(:,:) = e3t_m(:,:) + e3t_n(:,:,1) 131 131 ! 132 132 frq_m(:,:) = frq_m(:,:) + fraqsr_1lev(:,:) … … 141 141 ssv_m(:,:) = ssv_m(:,:) * zcoef ! 142 142 ssh_m(:,:) = ssh_m(:,:) * zcoef ! mean SSH [m] 143 IF( lk_vvl) e3t_m(:,:) = e3t_m(:,:) * zcoef ! mean vertical scale factor [m]143 IF( .NOT.ln_linssh ) e3t_m(:,:) = e3t_m(:,:) * zcoef ! mean vertical scale factor [m] 144 144 frq_m(:,:) = frq_m(:,:) * zcoef ! mean fraction of solar net radiation absorbed in the 1st T level [-] 145 145 ! … … 159 159 CALL iom_rstput( kt, nitrst, numrow, 'sss_m' , sss_m ) 160 160 CALL iom_rstput( kt, nitrst, numrow, 'ssh_m' , ssh_m ) 161 IF( lk_vvl) CALL iom_rstput( kt, nitrst, numrow, 'e3t_m' , e3t_m )161 IF( .NOT.ln_linssh ) CALL iom_rstput( kt, nitrst, numrow, 'e3t_m' , e3t_m ) 162 162 CALL iom_rstput( kt, nitrst, numrow, 'frq_m' , frq_m ) 163 163 ! … … 172 172 CALL iom_put( 'sss_m', sss_m ) 173 173 CALL iom_put( 'ssh_m', ssh_m ) 174 IF( lk_vvl) CALL iom_put( 'e3t_m', e3t_m )174 IF( .NOT.ln_linssh ) CALL iom_put( 'e3t_m', e3t_m ) 175 175 CALL iom_put( 'frq_m', frq_m ) 176 176 ENDIF … … 209 209 CALL iom_get( numror, jpdom_autoglo, 'sss_m' , sss_m ) ! " " salinity (T-point) 210 210 CALL iom_get( numror, jpdom_autoglo, 'ssh_m' , ssh_m ) ! " " height (T-point) 211 IF( lk_vvl )CALL iom_get( numror, jpdom_autoglo, 'e3t_m', e3t_m )211 IF( .NOT.ln_linssh ) CALL iom_get( numror, jpdom_autoglo, 'e3t_m', e3t_m ) 212 212 ! fraction of solar net radiation absorbed in 1st T level 213 213 IF( iom_varid( numror, 'frq_m', ldstop = .FALSE. ) > 0 ) THEN … … 226 226 sss_m(:,:) = zcoef * sss_m(:,:) 227 227 ssh_m(:,:) = zcoef * ssh_m(:,:) 228 IF( lk_vvl) e3t_m(:,:) = zcoef * e3t_m(:,:)228 IF( .NOT.ln_linssh ) e3t_m(:,:) = zcoef * e3t_m(:,:) 229 229 frq_m(:,:) = zcoef * frq_m(:,:) 230 230 ELSE … … 244 244 sss_m(:,:) = tsn(:,:,1,jp_sal) 245 245 ssh_m(:,:) = sshn(:,:) 246 IF( lk_vvl) e3t_m(:,:) = e3t_n(:,:,1)246 IF( .NOT.ln_linssh ) e3t_m(:,:) = e3t_n(:,:,1) 247 247 frq_m(:,:) = 1._wp 248 248 !
Note: See TracChangeset
for help on using the changeset viewer.