Changeset 8373 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/albedoice.F90

Timestamp:

2017-07-25T19:44:54+02:00 (7 years ago)

Author:

clem

Message:

remove most of wrk_alloc

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/albedoice.F90 (modified) (5 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/albedoice.F90

@@ -10 +10 @@
    !!   albedo_init   : initialisation of albedo computation
    !!----------------------------------------------------------------------
+   USE ice, ONLY : jpl
    USE phycst         ! physical constants
    USE in_out_manager ! I/O manager
@@ -82 +83 @@
       REAL(wp), INTENT(  out), DIMENSION(:,:,:) ::   pa_ice_cs           !  albedo of ice under clear    sky
       REAL(wp), INTENT(  out), DIMENSION(:,:,:) ::   pa_ice_os           !  albedo of ice under overcast sky
-      !!
+      !
       INTEGER  ::   ji, jj, jl                                           ! dummy loop indices
-      INTEGER  ::   ijpl                                                 ! number of ice categories (3rd dim of ice input arrays)
-      REAL(wp)                            ::   zswitch, z1_c1, z1_c2
-      REAL(wp)                            ::   zhref_pnd                                 
-      REAL(wp)                            ::   zalb_sm, zalb_sf, zalb_st ! albedo of snow melting, freezing, total
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zalb, zalb_it             ! intermediate variable & albedo of ice (snow free)
+      REAL(wp) ::   zswitch, z1_c1, z1_c2
+      REAL(wp) ::   zhref_pnd                                 
+      REAL(wp) ::   zalb_sm, zalb_sf, zalb_st ! albedo of snow melting, freezing, total
+      !
+      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zalb, zalb_it             ! intermediate variable & albedo of ice (snow free)
 !! MV MP
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zalb_pnd                  ! ponded sea ice albedo
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zalb_ice                  ! bare sea ice albedo
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zalb_snw                  ! snow-covered sea ice albedo
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zafrac_snw                ! relative snow fraction
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zafrac_ice                ! relative ice fraction
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zafrac_pnd                ! relative ice fraction (effective)
+      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zalb_pnd                  ! ponded sea ice albedo
+      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zalb_ice                  ! bare sea ice albedo
+      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zalb_snw                  ! snow-covered sea ice albedo
+      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zafrac_snw                ! relative snow fraction
+      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zafrac_ice                ! relative ice fraction
+      REAL(wp), DIMENSION(jpi,jpj,jpl) ::   zafrac_pnd                ! relative ice fraction (effective)
       !!
       !!---------------------------------------------------------------------
-
-      ijpl = SIZE( pt_ice, 3 )                     ! number of ice categories
-      
-      CALL wrk_alloc( jpi,jpj,ijpl, zalb, zalb_it )
-      CALL wrk_alloc( jpi,jpj,ijpl, zalb_pnd, zalb_ice, zalb_snw )
-      CALL wrk_alloc( jpi,jpj,ijpl, zalb_pnd, zafrac_snw, zafrac_ice, zafrac_pnd )
 
       IF( albd_init == 0 )   CALL albedo_init      ! initialization 
@@ -148 +143 @@
          ENDIF 
 
-         DO jl = 1, ijpl
+         DO jl = 1, jpl
             DO jj = 1, jpj
                DO ji = 1, jpi
@@ -224 +219 @@
          
          ! Overcast sky surface albedo (accounting for snow, ice melt ponds)
-         DO jl = 1, ijpl
+         DO jl = 1, jpl
             DO jj = 1, jpj
                DO ji = 1, jpi
@@ -309 +304 @@
 
       END SELECT
-      
-      CALL wrk_dealloc( jpi,jpj,ijpl, zalb, zalb_it )
-      CALL wrk_dealloc( jpi,jpj,ijpl, zalb_pnd, zalb_ice, zalb_snw )
-      CALL wrk_dealloc( jpi,jpj,ijpl, zalb_pnd, zafrac_snw, zafrac_ice, zafrac_pnd )
       !
    END SUBROUTINE albedo_ice