Changeset 1742 for trunk

Timestamp:

2009-11-20T09:46:06+01:00 (14 years ago)

Author:

smasson

Message:

sbccpl: bugfix of pqsr_ice + energy due to iceberg melting, see ticket:607

File:

• trunk/NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (7 diffs)

trunk/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -45 +45 @@
    USE lib_mpp         ! distribued memory computing library
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
-   USE phycst, ONLY : xlsn, rhosn
+   USE phycst, ONLY : xlsn, rhosn, xlic, rhoic
 #if defined key_cpl_carbon_cycle
    USE p4zflx, ONLY : oce_co2
@@ -1038 +1038 @@
       INTEGER ::   isec, info       ! temporary integer
       REAL(wp)::   zcoef, ztsurf    ! temporary scalar
-      REAL(wp), DIMENSION(jpi,jpj)::   zsnow    ! snow precipitation 
+      REAL(wp), DIMENSION(jpi,jpj    )::   zsnow    ! snow precipitation 
+      REAL(wp), DIMENSION(jpi,jpj,jpl)::   zicefr   ! ice fraction 
       !!----------------------------------------------------------------------
+      zicefr(:,:,1) = 1.- p_frld(:,:,1)
       !
       !                                                      ! ========================= !
@@ -1054 +1056 @@
          zsnow   (:,:) = frcv(:,:,jpr_snow)
       CASE( 'oce and ice'   )   ! received fields: jpr_sbpr, jpr_semp, jpr_oemp
-         pemp_tot(:,:) = p_frld(:,:,1) * frcv(:,:,jpr_oemp) + (1.- p_frld(:,:,1)) * frcv(:,:,jpr_sbpr) 
+         pemp_tot(:,:) = p_frld(:,:,1) * frcv(:,:,jpr_oemp) + zicefr(:,:,1) * frcv(:,:,jpr_sbpr) 
          pemp_ice(:,:) = frcv(:,:,jpr_semp)
          zsnow   (:,:) = - frcv(:,:,jpr_semp) + frcv(:,:,jpr_ievp)
@@ -1088 +1090 @@
          pqns_ice(:,:,1) = frcv(:,:,jpr_qnsice)
       CASE( 'oce and ice' )       ! the total flux is computed from ocean and ice fluxes
-         pqns_tot(:,:  ) =  p_frld(:,:,1) * frcv(:,:,jpr_qnsoce) + ( 1.- p_frld(:,:,1) ) * frcv(:,:,jpr_qnsice)
+         pqns_tot(:,:  ) =  p_frld(:,:,1) * frcv(:,:,jpr_qnsoce) + zicefr(:,:,1) * frcv(:,:,jpr_qnsice)
          pqns_ice(:,:,1) =  frcv(:,:,jpr_qnsice)
       CASE( 'mixed oce-ice' )     ! the ice flux is cumputed from the total flux, the SST and ice informations
          pqns_tot(:,:  ) = frcv(:,:,jpr_qnsmix)
          pqns_ice(:,:,1) = frcv(:,:,jpr_qnsmix)    &
-            &            + frcv(:,:,jpr_dqnsdt) * ( pist(:,:,1) - ( (rt0 + psst(:,:  ) ) *        p_frld(:,:,1)   &
-            &                                                   +          pist(:,:,1)   * ( 1. - p_frld(:,:,1) ) ) )
+            &            + frcv(:,:,jpr_dqnsdt) * ( pist(:,:,1) - ( (rt0 + psst(:,:  ) ) * p_frld(:,:,1)   &
+            &                                                   +          pist(:,:,1)   * zicefr(:,:,1) ) )
       END SELECT
       !                                                           ! snow melting heat flux ....
-      !   energy for melting solid precipitation over free ocean
+      !   energy for melting solid precipitation over ice-free ocean
       zcoef = xlsn / rhosn
       pqns_tot(:,:) = pqns_tot(:,:) - p_frld(:,:,1) * zsnow(:,:) * zcoef
@@ -1107 +1109 @@
 !!
 !! similar job should be done for snow and precipitation temperature
+      !                                                           ! Iceberg melting heat flux ....
+      !   energy for iceberg melting 
+      IF( srcv(jpr_cal)%laction ) THEN 
+         zcoef = xlic / rhoic
+         pqns_tot(:,:) = pqns_tot(:,:) - frcv(:,:,jpr_cal)(:,:) * zcoef
+      ENDIF
 
       !                                                      ! ========================= !
@@ -1115 +1123 @@
          pqsr_ice(:,:,1) = frcv(:,:,jpr_qsrice)
       CASE( 'oce and ice' )
-         pqsr_tot(:,:  ) =  p_frld(:,:,1) * frcv(:,:,jpr_qsroce) + ( 1.- p_frld(:,:,1) ) * frcv(:,:,jpr_qsrice)
+         pqsr_tot(:,:  ) =  p_frld(:,:,1) * frcv(:,:,jpr_qsroce) + zicefr(:,:,1) * frcv(:,:,jpr_qsrice)
          pqsr_ice(:,:,1) =  frcv(:,:,jpr_qsrice)
       CASE( 'mixed oce-ice' )
@@ -1122 +1130 @@
 !       ( see OASIS3 user guide, 5th edition, p39 )
          pqsr_ice(:,:,1) = frcv(:,:,jpr_qsrmix) * ( 1.- palbi(:,:,1) )   &
-            &            / (  1.- ( albedo_oce_mix(:,:  ) * ( 1.- p_frld(:,:,1) )   &
-            &                     + palbi         (:,:,1) *       p_frld(:,:,1)   )  )
+            &            / (  1.- ( albedo_oce_mix(:,:  ) * p_frld(:,:,1)   &
+            &                     + palbi         (:,:,1) * zicefr(:,:,1) ) )
       END SELECT