Changeset 1860 for trunk/NEMO/OPA_SRC/SBC/sbccpl.F90
- Timestamp:
- 2010-05-08T11:23:48+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbccpl.F90
r1833 r1860 23 23 USE sbc_oce ! Surface boundary condition: ocean fields 24 24 USE sbc_ice ! Surface boundary condition: ice fields 25 USE phycst ! physical constants 25 26 #if defined key_lim3 26 27 USE par_ice ! ice parameters … … 45 46 USE lib_mpp ! distribued memory computing library 46 47 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 47 USE phycst, ONLY : xlsn, rhosn, xlic, rhoic48 48 #if defined key_cpl_carbon_cycle 49 49 USE p4zflx, ONLY : oce_co2 … … 715 715 IF( srcv(jpr_qnsoce)%laction ) qns(:,:) = frcv(:,:,jpr_qnsoce) 716 716 IF( srcv(jpr_qnsmix)%laction ) qns(:,:) = frcv(:,:,jpr_qnsmix) 717 ! energy for melting solid precipitation over free ocean 718 zcoef = xlsn / rhosn 719 qns(:,:) = qns(:,:) - frcv(:,:,jpr_snow) * zcoef 717 qns(:,:) = qns(:,:) - frcv(:,:,jpr_snow) * lfus ! add the latent heat of solid precip. melting 718 720 719 ! ! solar flux over the ocean (qsr) 721 720 IF( srcv(jpr_qsroce)%laction ) qsr(:,:) = frcv(:,:,jpr_qsroce) … … 1118 1117 & + pist(:,:,1) * zicefr(:,:,1) ) ) 1119 1118 END SELECT 1120 ! ! snow melting heat flux .... 1121 ! energy for melting solid precipitation over ice-free ocean 1122 zcoef = xlsn / rhosn 1123 ztmp(:,:) = p_frld(:,:,1) * zsnow(:,:) * zcoef 1124 pqns_tot(:,:) = pqns_tot(:,:) - ztmp(:,:) 1119 ztmp(:,:) = p_frld(:,:,1) * zsnow(:,:) * lfus ! add the latent heat of solid precip. melting 1120 pqns_tot(:,:) = pqns_tot(:,:) - ztmp(:,:) ! over free ocean 1125 1121 IF( lk_diaar5 ) CALL iom_put( 'hflx_snow_cea', ztmp + zsnow(:,:) * zcptn(:,:) ) ! heat flux from snow (cell average) 1126 1122 !!gm … … 1131 1127 !! 1132 1128 !! similar job should be done for snow and precipitation temperature 1133 ! ! Iceberg melting heat flux .... 1134 ! energy for iceberg melting 1135 IF( srcv(jpr_cal)%laction ) THEN 1136 zcoef = xlic / rhoic 1137 ztmp(:,:) = frcv(:,:,jpr_cal) * zcoef 1129 ! 1130 IF( srcv(jpr_cal)%laction ) THEN ! Iceberg melting 1131 ztmp(:,:) = frcv(:,:,jpr_cal) * lfus ! add the latent heat of iceberg melting 1138 1132 pqns_tot(:,:) = pqns_tot(:,:) - ztmp(:,:) 1139 1133 IF( lk_diaar5 ) CALL iom_put( 'hflx_cal_cea', ztmp + frcv(:,:,jpr_cal) * zcptn(:,:) ) ! heat flux from calving
Note: See TracChangeset
for help on using the changeset viewer.