Changeset 3152 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

Timestamp:

2011-11-18T10:19:26+01:00 (12 years ago)

Author:

smasson

Message:

dev_NEMO_MERGE_2011: new dynamical allocation in IOM and SBC

File:

• branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (15 diffs)

branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -47 +47 @@
    USE iom             ! NetCDF library
    USE lib_mpp         ! distribued memory computing library
+   USE wrk_nemo_2      ! work arrays
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
 #if defined key_cpl_carbon_cycle
@@ -217 +218 @@
       !!              * initialise the OASIS coupler
       !!----------------------------------------------------------------------
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   zacs => wrk_2d_3 , zaos => wrk_2d_4   ! clear & overcast sky albedos
-      !!
       INTEGER, INTENT(in) ::   k_ice    ! ice management in the sbc (=0/1/2/3)
       !!
       INTEGER ::   jn   ! dummy loop index
+      REAL(wp), POINTER, DIMENSION(:,:) ::   zacs, zaos
       !!
       NAMELIST/namsbc_cpl/  sn_snd_temp, sn_snd_alb   , sn_snd_thick, sn_snd_crt   , sn_snd_co2,   &
@@ -229 +228 @@
       !!---------------------------------------------------------------------
 
-      IF( wrk_in_use(2, 3,4) ) THEN
-         CALL ctl_stop('sbc_cpl_init: requested workspace arrays unavailable')   ;   RETURN
-      ENDIF
+      CALL wrk_alloc( jpi,jpj, zacs, zaos )
 
       ! ================================ !
@@ -618 +615 @@
          &   CALL ctl_stop( 'sbc_cpl_init: diurnal cycle reconstruction (ln_dm2dc) needs daily couping for solar radiation' )
 
-      IF( wrk_not_released(2, 3,4) )   CALL ctl_stop('sbc_cpl_init: failed to release workspace arrays')
+      CALL wrk_dealloc( jpi,jpj, zacs, zaos )
       !
    END SUBROUTINE sbc_cpl_init
@@ -665 +662 @@
       !!                        emp = emps   evap. - precip. (- runoffs) (- calving) ('ocean only case)
       !!----------------------------------------------------------------------
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   ztx => wrk_2d_1 , zty => wrk_2d_2
-      !!
       INTEGER, INTENT(in) ::   kt       ! ocean model time step index
       INTEGER, INTENT(in) ::   k_fsbc   ! frequency of sbc (-> ice model) computation 
@@ -680 +674 @@
       REAL(wp) ::   zcdrag = 1.5e-3        ! drag coefficient
       REAL(wp) ::   zzx, zzy               ! temporary variables
-      !!----------------------------------------------------------------------
-
-      IF( wrk_in_use(2, 1,2) ) THEN
-         CALL ctl_stop('sbc_cpl_rcv: requested workspace arrays unavailable')   ;   RETURN
-      ENDIF
+      REAL(wp), POINTER, DIMENSION(:,:) ::   ztx, zty 
+      !!----------------------------------------------------------------------
+
+      CALL wrk_alloc( jpi,jpj, ztx, zty )
 
       IF( kt == nit000 )   CALL sbc_cpl_init( k_ice )          ! initialisation
@@ -867 +860 @@
       ENDIF
       !
-      IF( wrk_not_released(2, 1,2) )   CALL ctl_stop('sbc_cpl_rcv: failed to release workspace arrays')
+      CALL wrk_dealloc( jpi,jpj, ztx, zty )
       !
    END SUBROUTINE sbc_cpl_rcv
@@ -905 +898 @@
       !! ** Action  :   return ptau_i, ptau_j, the stress over the ice at cp_ice_msh point
       !!----------------------------------------------------------------------
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   ztx => wrk_2d_1 , zty => wrk_2d_2
-      !!
       REAL(wp), INTENT(out), DIMENSION(:,:) ::   p_taui   ! i- & j-components of atmos-ice stress [N/m2]
       REAL(wp), INTENT(out), DIMENSION(:,:) ::   p_tauj   ! at I-point (B-grid) or U & V-point (C-grid)
@@ -913 +903 @@
       INTEGER ::   ji, jj                          ! dummy loop indices
       INTEGER ::   itx                             ! index of taux over ice
-      !!----------------------------------------------------------------------
-
-      IF( wrk_in_use(2, 1,2) ) THEN
-         CALL ctl_stop('sbc_cpl_ice_tau: requested workspace arrays unavailable')   ;   RETURN
-      ENDIF
+      REAL(wp), POINTER, DIMENSION(:,:) ::   ztx, zty 
+      !!----------------------------------------------------------------------
+      
+      CALL wrk_alloc( jpi,jpj, ztx, zty )
 
       IF( srcv(jpr_itx1)%laction ) THEN   ;   itx =  jpr_itx1   
@@ -1079 +1068 @@
       ENDIF
       !   
-      IF( wrk_not_released(2, 1,2) )   CALL ctl_stop('sbc_cpl_ice_tau: failed to release workspace arrays')
+      CALL wrk_dealloc( jpi,jpj, ztx, zty )
       !
    END SUBROUTINE sbc_cpl_ice_tau
@@ -1124 +1113 @@
       !!                   sprecip             solid precipitation over the ocean  
       !!----------------------------------------------------------------------
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   zcptn  => wrk_2d_2   ! rcp * tsn(:,:,1,jp_tem)
-      USE wrk_nemo, ONLY:   ztmp   => wrk_2d_3   ! temporary array
-      USE wrk_nemo, ONLY:   zicefr => wrk_2d_4   ! ice fraction 
-      !!
       REAL(wp), INTENT(in   ), DIMENSION(:,:)   ::   p_frld     ! lead fraction                [0 to 1]
       ! optional arguments, used only in 'mixed oce-ice' case
@@ -1136 +1120 @@
       !
       INTEGER ::   jl   ! dummy loop index
-      !!----------------------------------------------------------------------
-
-      IF( wrk_in_use(2, 2,3,4) ) THEN
-         CALL ctl_stop('sbc_cpl_ice_flx: requested workspace arrays unavailable')   ;   RETURN
-      ENDIF
+      REAL(wp), POINTER, DIMENSION(:,:) ::   zcptn, ztmp, zicefr
+      !!----------------------------------------------------------------------
+
+      CALL wrk_alloc( jpi,jpj, zcptn, ztmp, zicefr )
 
       zicefr(:,:) = 1.- p_frld(:,:)
@@ -1316 +1299 @@
       END SELECT
 
-      IF( wrk_not_released(2, 2,3,4) ) CALL ctl_stop('sbc_cpl_ice_flx: failed to release workspace arrays')
+      CALL wrk_dealloc( jpi,jpj, zcptn, ztmp, zicefr )
       !
    END SUBROUTINE sbc_cpl_ice_flx
@@ -1330 +1313 @@
       !!              all the needed fields (as defined in sbc_cpl_init)
       !!----------------------------------------------------------------------
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   zfr_l => wrk_2d_1   ! 1. - fr_i(:,:)
-      USE wrk_nemo, ONLY:   ztmp1 => wrk_2d_2 , ztmp2 => wrk_2d_3
-      USE wrk_nemo, ONLY:   ztmp3 => wrk_3d_1 , ztmp4 => wrk_3d_2
-      USE wrk_nemo, ONLY:   zotx1 => wrk_2d_4 , zoty1 => wrk_2d_5 , zotz1 => wrk_2d_6
-      USE wrk_nemo, ONLY:   zitx1 => wrk_2d_7 , zity1 => wrk_2d_8 , zitz1 => wrk_2d_9
-      !
       INTEGER, INTENT(in) ::   kt
       !
       INTEGER ::   ji, jj, jl   ! dummy loop indices
       INTEGER ::   isec, info   ! local integer
-      !!----------------------------------------------------------------------
-
-      IF( wrk_in_use(2, 1,2,3,4,5,6,7,8,9) .OR. wrk_in_use(3, 1,2)  ) THEN
-         CALL ctl_stop('sbc_cpl_snd: requested workspace arrays are unavailable')   ;   RETURN
-      ENDIF
+      REAL(wp), POINTER, DIMENSION(:,:)   ::   ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmp3, ztmp4   
+      !!----------------------------------------------------------------------
+
+      CALL wrk_alloc( jpi,jpj, zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1 )
+      CALL wrk_alloc( jpi,jpj,jpl, ztmp3, ztmp4 )
 
       isec = ( kt - nit000 ) * NINT(rdttra(1))        ! date of exchanges
@@ -1565 +1542 @@
       ENDIF
       !
-      IF( wrk_not_released(2, 1,2,3,4,5,6,7,8,9) .OR. wrk_not_released(3, 1,2) )   CALL ctl_stop('sbc_cpl_snd: failed to release workspace arrays')
+      CALL wrk_dealloc( jpi,jpj, zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1 )
+      CALL wrk_dealloc( jpi,jpj,jpl, ztmp3, ztmp4 )
       !
    END SUBROUTINE sbc_cpl_snd