Changeset 4990 for trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90
- Timestamp:
- 2014-12-15T17:42:49+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90
r4624 r4990 19 19 USE phycst ! physical constants 20 20 USE sbc_oce ! surface boundary condition variables 21 USE sbcisf ! PM we could remove it I think 21 22 USE closea ! closed seas 22 23 USE fldread ! read input field at current time step … … 24 25 USE iom ! I/O module 25 26 USE lib_mpp ! MPP library 27 USE eosbn2 28 USE wrk_nemo ! Memory allocation 26 29 27 30 IMPLICIT NONE … … 98 101 INTEGER :: z_err = 0 ! dummy integer for error handling 99 102 !!---------------------------------------------------------------------- 103 REAL(wp), DIMENSION(:,:), POINTER :: ztfrz ! freezing point used for temperature correction 104 ! 105 CALL wrk_alloc( jpi,jpj, ztfrz) 106 100 107 ! 101 108 IF( kt == nit000 ) CALL sbc_rnf_init ! Read namelist and allocate structures … … 134 141 WHERE( sf_t_rnf(1)%fnow(:,:,1) == -999._wp ) ! if missing data value use SST as runoffs temperature 135 142 rnf_tsc(:,:,jp_tem) = sst_m(:,:) * rnf(:,:) * r1_rau0 143 END WHERE 144 WHERE( sf_t_rnf(1)%fnow(:,:,1) == -222._wp ) ! where fwf comes from melting of ice shelves or iceberg 145 ztfrz(:,:) = -1.9 !tfreez( sss_m(:,:) ) !PM to be discuss (trouble if sensitivity study) 146 rnf_tsc(:,:,jp_tem) = ztfrz(:,:) * rnf(:,:) * r1_rau0 - rnf(:,:) * lfusisf * r1_rau0_rcp 136 147 END WHERE 137 148 ELSE ! use SST as runoffs temperature … … 175 186 CALL iom_rstput( kt, nitrst, numrow, 'rnf_sc_b', rnf_tsc(:,:,jp_sal) ) 176 187 ENDIF 188 CALL wrk_dealloc( jpi,jpj, ztfrz) 177 189 ! 178 190 END SUBROUTINE sbc_rnf
Note: See TracChangeset
for help on using the changeset viewer.