Changeset 1756 for trunk/NEMO/OPA_SRC/SBC/sbccpl.F90

Timestamp:

2009-11-25T15:15:20+01:00 (14 years ago)

Author:

smasson

Message:

implement AR5 diagnostics, see ticket:610

File:

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

trunk/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -49 +49 @@
    USE p4zflx, ONLY : oce_co2
 #endif
+   USE diaar5, ONLY :   lk_diaar5
    IMPLICIT NONE
    PRIVATE
@@ -1038 +1039 @@
       INTEGER ::   isec, info       ! temporary integer
       REAL(wp)::   zcoef, ztsurf    ! temporary scalar
+      REAL(wp), DIMENSION(jpi,jpj    )::   zcptn    ! rcp * tn(:,:,1)
+      REAL(wp), DIMENSION(jpi,jpj    )::   ztmp     ! temporary array
       REAL(wp), DIMENSION(jpi,jpj    )::   zsnow    ! snow precipitation 
       REAL(wp), DIMENSION(jpi,jpj,jpl)::   zicefr   ! ice fraction 
       !!----------------------------------------------------------------------
       zicefr(:,:,1) = 1.- p_frld(:,:,1)
+      IF( lk_diaar5 )   zcptn(:,:) = rcp * tn(:,:,1)
       !
       !                                                      ! ========================= !
@@ -1055 +1059 @@
          pemp_ice(:,:) = frcv(:,:,jpr_ievp) - frcv(:,:,jpr_snow)
          zsnow   (:,:) = frcv(:,:,jpr_snow)
+                           CALL iom_put( 'rain'         , frcv(:,:,jpr_rain)              )   ! liquid precipitation 
+         IF( lk_diaar5 )   CALL iom_put( 'hflx_rain_cea', frcv(:,:,jpr_rain) * zcptn(:,:) )   ! heat flux from liq. precip. 
+         ztmp(:,:) = frcv(:,:,jpr_tevp) - frcv(:,:,jpr_ievp) * zicefr(:,:,1)
+                           CALL iom_put( 'evap_ao_cea'  , ztmp                            )   ! ice-free oce evap (cell average)
+         IF( lk_diaar5 )   CALL iom_put( 'hflx_evap_cea', ztmp(:,:         ) * zcptn(:,:) )   ! heat flux from from evap (cell ave)
       CASE( 'oce and ice'   )   ! received fields: jpr_sbpr, jpr_semp, jpr_oemp
          pemp_tot(:,:) = p_frld(:,:,1) * frcv(:,:,jpr_oemp) + zicefr(:,:,1) * frcv(:,:,jpr_sbpr) 
@@ -1061 +1070 @@
       END SELECT
       psprecip(:,:) = - pemp_ice(:,:)
-      CALL iom_put( 'snowpre', psprecip )     ! Snow precipitation 
+      CALL iom_put( 'snowpre'    , zsnow                               )   ! Snow
+      CALL iom_put( 'snow_ao_cea', zsnow(:,:         ) * p_frld(:,:,1) )   ! Snow        over ice-free ocean  (cell average)
+      CALL iom_put( 'snow_ai_cea', zsnow(:,:         ) * zicefr(:,:,1) )   ! Snow        over sea-ice         (cell average)
+      CALL iom_put( 'subl_ai_cea', frcv (:,:,jpr_ievp) * zicefr(:,:,1) )   ! Sublimation over sea-ice         (cell average)
       !   
       !                                                           ! runoffs and calving (put in emp_tot)
-      IF( srcv(jpr_rnf)%laction )   pemp_tot(:,:) = pemp_tot(:,:) -      frcv(:,:,jpr_rnf)
-      IF( srcv(jpr_cal)%laction )   pemp_tot(:,:) = pemp_tot(:,:) - ABS( frcv(:,:,jpr_cal) )
+      IF( srcv(jpr_rnf)%laction ) THEN 
+         pemp_tot(:,:) = pemp_tot(:,:) - frcv(:,:,jpr_rnf)
+                           CALL iom_put( 'runoffs'      , frcv(:,:,jpr_rnf )              )   ! rivers
+         IF( lk_diaar5 )   CALL iom_put( 'hflx_rnf_cea' , frcv(:,:,jpr_rnf ) * zcptn(:,:) )   ! heat flux from rivers
+      ENDIF
+      IF( srcv(jpr_cal)%laction ) THEN 
+         pemp_tot(:,:) = pemp_tot(:,:) - frcv(:,:,jpr_cal)
+         CALL iom_put( 'calving', frcv(:,:,jpr_cal) )
+      ENDIF
       !
 !!gm :  this seems to be internal cooking, not sure to need that in a generic interface 
@@ -1101 +1120 @@
       !   energy for melting solid precipitation over ice-free ocean
       zcoef = xlsn / rhosn
-      pqns_tot(:,:) = pqns_tot(:,:) - p_frld(:,:,1) * zsnow(:,:) * zcoef
+      ztmp(:,:) = p_frld(:,:,1) * zsnow(:,:) * zcoef
+      pqns_tot(:,:) = pqns_tot(:,:) - ztmp(:,:)
+      IF( lk_diaar5 )   CALL iom_put( 'hflx_snow_cea', ztmp + zsnow(:,:) * zcptn(:,:) )   ! heat flux from snow (cell average)
 !!gm
 !!    currently it is taken into account in leads budget but not in the qns_tot, and thus not in 
@@ -1113 +1134 @@
       IF( srcv(jpr_cal)%laction ) THEN 
          zcoef = xlic / rhoic
-         pqns_tot(:,:) = pqns_tot(:,:) - frcv(:,:,jpr_cal) * zcoef
+         ztmp(:,:) = frcv(:,:,jpr_cal) * zcoef
+         pqns_tot(:,:) = pqns_tot(:,:) - ztmp(:,:)
+         IF( lk_diaar5 )   CALL iom_put( 'hflx_cal_cea', ztmp + frcv(:,:,jpr_cal) * zcptn(:,:) )   ! heat flux from calving
       ENDIF