- Timestamp:
- 2010-05-08T11:23:48+02:00 (14 years ago)
- Location:
- trunk/NEMO/OPA_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DOM/phycst.F90
r1739 r1860 4 4 !! Definition of of both ocean and ice parameters used in the code 5 5 !!===================================================================== 6 !! History : ! 90-10 (C. Levy - G. Madec) Original code 7 !! ! 91-11 (G. Madec) 8 !! ! 91-12 (M. Imbard) 9 !! 8.5 ! 02-08 (G. Madec, C. Ethe) F90, add ice constants 10 !! 9.0 ! 06-08 (G. Madec) style 6 !! History : OPA ! 1990-10 (C. Levy, G. Madec) Original code 7 !! ! 1991-11 (G. Madec, M. Imbard) 8 !! NEMO 1.0 ! 2002-08 (G. Madec, C. Ethe) F90, add ice constants 9 !! - ! 2006-08 (G. Madec) style 11 10 !!---------------------------------------------------------------------- 12 11 … … 59 58 rcdic = 2.034396_wp , & !: thermal conductivity of fresh ice 60 59 cpic = 2067.0 , & 61 ! add the following lines62 60 lsub = 2.834e+6 , & !: pure ice latent heat of sublimation (J.kg-1) 63 61 lfus = 0.334e+6 , & !: latent heat of fusion of fresh ice (J.kg-1) … … 69 67 rcpsn = 6.9069e+5_wp, & !: density times specific heat for snow 70 68 rcpic = 1.8837e+6_wp, & !: volumetric latent heat fusion of sea ice 71 xlsn = 110.121e+6_wp , & !: volumetric latent heat fusion of snow 72 xlic = 300.33e+6_wp , & !: volumetric latent heat fusion of ice 69 lfus = 0.3337e+6 , & !: latent heat of fusion of fresh ice (J.kg-1) 70 xlsn = 110.121e+6_wp , & !: = lfus * rhosn, volumetric latent heat fusion of snow 71 xlic = 300.33e+6_wp , & !: = lfus * rhosn, volumetric latent heat fusion of ice 73 72 xsn = 2.8e+6 , & !: latent heat of sublimation of snow 74 73 rhoic = 900._wp , & !: volumic mass of sea ice (kg/m3) … … 83 82 stefan = 5.67e-8_wp !: Stefan-Boltzmann constant 84 83 !!---------------------------------------------------------------------- 85 !! OPA 9.0 , LOCEAN-IPSL (2005)84 !! NEMO/OPA 3.2 , LOCEAN-IPSL (2010) 86 85 !! $Id$ 87 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt86 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 88 87 !!---------------------------------------------------------------------- 89 88 … … 160 159 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) = ro0cpr = ', ro0cpr 161 160 161 #if defined key_lim3 162 xlsn = lfus * rhosn ! volumetric latent heat fusion of snow [J/m3] 163 #else 164 lfus = xlsn / rhosn ! latent heat of fusion of fresh ice 165 #endif 166 162 167 IF(lwp) THEN 163 168 WRITE(numout,*) -
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.