Changeset 15440 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/ICE/icewri.F90

Timestamp:

2021-10-23T12:18:24+02:00 (3 years ago)

Author:

cetlod

Message:

dev_PISCO : merge with trunk@15439

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/ICE/icewri.F90 (modified) (4 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/ICE/icewri.F90

@@ -20 +20 @@
    USE ice            ! sea-ice: variables
    USE icevar         ! sea-ice: operations
+   USE icealb , ONLY : rn_alb_oce
    !
    USE ioipsl         !
@@ -53 +54 @@
       REAL(wp) ::   z2da, z2db, zrho1, zrho2
       REAL(wp) ::   zmiss_val       ! missing value retrieved from xios
-      REAL(wp), DIMENSION(jpi,jpj)     ::   z2d, zfast                     ! 2D workspace
+      REAL(wp), DIMENSION(jpi,jpj)     ::   z2d                            ! 2D workspace
       REAL(wp), DIMENSION(jpi,jpj)     ::   zmsk00, zmsk05, zmsk15, zmsksn ! O%, 5% and 15% concentration mask and snow mask
       REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zmsk00l, zmsksnl               ! cat masks
+      REAL(wp), DIMENSION(:,:), ALLOCATABLE ::   zfast, zalb, zmskalb      ! 2D workspace
       !
       ! Global ice diagnostics (SIMIP)
@@ -131 +133 @@
       IF( iom_use('vice'    ) )   CALL iom_put( 'vice'   , v_ice    )                                                       ! ice velocity v
       !
-      IF( iom_use('icevel') .OR. iom_use('fasticepres') ) THEN                                                              ! module of ice velocity
+      IF( iom_use('icevel') .OR. iom_use('fasticepres') ) THEN                                                              ! module of ice velocity & fast ice
+         ALLOCATE( zfast(jpi,jpj) )
          DO_2D( 0, 0, 0, 0 )
             z2da  = u_ice(ji,jj) + u_ice(ji-1,jj)
@@ -144 +147 @@
          END WHERE
          CALL iom_put( 'fasticepres', zfast )
-      ENDIF
-
+         DEALLOCATE( zfast )
+      ENDIF
+      !
+      IF( iom_use('icealb') .OR. iom_use('albedo') ) THEN                                                                   ! ice albedo and surface albedo
+         ALLOCATE( zalb(jpi,jpj), zmskalb(jpi,jpj) )
+         ! ice albedo
+         WHERE( at_i_b < 1.e-03 )
+            zmskalb(:,:) = 0._wp
+            zalb   (:,:) = rn_alb_oce
+         ELSEWHERE
+            zmskalb(:,:) = 1._wp
+            zalb   (:,:) = SUM( alb_ice * a_i_b, dim=3 ) / at_i_b
+         END WHERE
+         CALL iom_put( 'icealb' , zalb * zmskalb + zmiss_val * ( 1._wp - zmskalb ) )
+         ! ice+ocean albedo
+         zalb(:,:) = SUM( alb_ice * a_i_b, dim=3 ) + rn_alb_oce * ( 1._wp - at_i_b )
+         CALL iom_put( 'albedo' , zalb )
+         DEALLOCATE( zalb, zmskalb )
+      ENDIF
+      !
       ! --- category-dependent fields --- !
       IF( iom_use('icemask_cat' ) )   CALL iom_put( 'icemask_cat' ,                  zmsk00l                                   ) ! ice mask 0%