- Timestamp:
- 2016-11-16T16:04:30+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r7233_CMIP6_diags_trunk_version/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90
r6140 r7240 92 92 INTEGER :: ji, jj ! loop index 93 93 REAL(wp), DIMENSION (:,:), POINTER :: zt_frz, zdep ! freezing temperature (zt_frz) at depth (zdep) 94 REAL(wp), DIMENSION(:,:,:), POINTER :: zfwfisf3d, zqhcisf3d, zqlatisf3d 95 REAL(wp), DIMENSION(:,: ), POINTER :: zqhcisf2d 94 96 !!--------------------------------------------------------------------- 95 97 ! … … 161 163 CALL lbc_lnk(risf_tsc(:,:,jp_tem),'T',1.) 162 164 CALL lbc_lnk(risf_tsc(:,:,jp_sal),'T',1.) 163 CALL lbc_lnk(fwfisf(:,:) ,'T',1.) 164 CALL lbc_lnk(qisf(:,:) ,'T',1.) 165 CALL lbc_lnk(fwfisf(:,:) ,'T',1.) 166 CALL lbc_lnk(qisf(:,:) ,'T',1.) 167 168 !============================================================================================================================================= 169 IF ( iom_use('fwfisf3d') .OR. iom_use('qlatisf3d') .OR. iom_use('qhcisf3d') .OR. iom_use('qhcisf')) THEN 170 CALL wrk_alloc( jpi,jpj,jpk, zfwfisf3d, zqhcisf3d, zqlatisf3d ) 171 CALL wrk_alloc( jpi,jpj, zqhcisf2d ) 172 173 zfwfisf3d(:,:,:) = 0.0_wp ! 3d ice shelf melting (kg/m2/s) 174 zqhcisf3d(:,:,:) = 0.0_wp ! 3d heat content flux (W/m2) 175 zqlatisf3d(:,:,:)= 0.0_wp ! 3d ice shelf melting latent heat flux (W/m2) 176 zqhcisf2d(:,:) = fwfisf(:,:) * zt_frz * rcp ! 2d heat content flux (W/m2) 177 178 DO jj = 1,jpj 179 DO ji = 1,jpi 180 ikt = misfkt(ji,jj) 181 ikb = misfkb(ji,jj) 182 DO jk = ikt, ikb - 1 183 zfwfisf3d (ji,jj,jk) = zfwfisf3d (ji,jj,jk) + fwfisf (ji,jj) * r1_hisf_tbl(ji,jj) * e3t_n(ji,jj,jk) 184 zqhcisf3d (ji,jj,jk) = zqhcisf3d (ji,jj,jk) + zqhcisf2d(ji,jj) * r1_hisf_tbl(ji,jj) * e3t_n(ji,jj,jk) 185 zqlatisf3d(ji,jj,jk) = zqlatisf3d(ji,jj,jk) + qisf (ji,jj) * r1_hisf_tbl(ji,jj) * e3t_n(ji,jj,jk) 186 END DO 187 zfwfisf3d (ji,jj,jk) = zfwfisf3d (ji,jj,jk) + fwfisf (ji,jj) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) * e3t_n(ji,jj,jk) 188 zqhcisf3d (ji,jj,jk) = zqhcisf3d (ji,jj,jk) + zqhcisf2d(ji,jj) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) * e3t_n(ji,jj,jk) 189 zqlatisf3d(ji,jj,jk) = zqlatisf3d(ji,jj,jk) + qisf (ji,jj) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) * e3t_n(ji,jj,jk) 190 END DO 191 END DO 192 193 CALL iom_put('fwfisf3d' , zfwfisf3d (:,:,:)) 194 CALL iom_put('qlatisf3d', zqlatisf3d(:,:,:)) 195 CALL iom_put('qhcisf3d' , zqhcisf3d (:,:,:)) 196 CALL iom_put('qhcisf' , zqhcisf2d (:,: )) 197 198 CALL wrk_dealloc( jpi,jpj,jpk, zfwfisf3d, zqhcisf3d, zqlatisf3d ) 199 CALL wrk_dealloc( jpi,jpj, zqhcisf2d ) 200 END IF 201 202 ! output 203 CALL iom_put('qlatisf' , qisf) 204 CALL iom_put('fwfisf', fwfisf) 205 !============================================================================================================================================= 165 206 166 207 IF( kt == nit000 ) THEN ! set the forcing field at nit000 - 1 ! … … 177 218 END IF 178 219 ! 179 ! output180 CALL iom_put('qisf' , qisf)181 CALL iom_put('fwfisf', fwfisf)182 183 220 ! deallocation 184 221 CALL wrk_dealloc( jpi,jpj, zt_frz, zdep )
Note: See TracChangeset
for help on using the changeset viewer.