New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3184 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90 – NEMO

Ignore:
Timestamp:
2011-11-25T16:42:54+01:00 (12 years ago)
Author:
cetlod
Message:

New dynamic allocation in NST modules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90

    r2715 r3184  
    1010   USE trc 
    1111   USE lib_mpp 
     12   USE wrk_nemo_2 
    1213 
    1314   IMPLICIT NONE 
     
    2930      !!--------------------------------------------- 
    3031#include "domzgr_substitute.h90" 
    31       USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 
    32       USE wrk_nemo, ONLY: wrk_2d_1 
    33       USE wrk_nemo, ONLY: wrk_4d_1, wrk_4d_2, wrk_4d_3, wrk_4d_4 
    3432      !!  
    3533      INTEGER :: ji,jj,jk,jl 
     
    4139 
    4240#if defined SPONGE_TOP 
    43       IF( wrk_in_use(4, 1,2,3,4) .OR. wrk_in_use(2, 1) ) THEN 
    44          CALL ctl_stop('Agrif_Sponge_trc : requested workspace arrays unavailable') 
    45          RETURN 
    46       ENDIF 
    47       localviscsponge => wrk_2d_1 
    48       trbdiff(:,:,:,:) => wrk_4d_1(:,:,:,1:jptra) 
    49       ztru   (:,:,:,:) => wrk_4d_2(:,:,:,1:jptra) 
    50       ztrv   (:,:,:,:) => wrk_4d_3(:,:,:,1:jptra) 
    51       ztab   (:,:,:,:) => wrk_4d_4(:,:,:,1:jptra) 
     41      CALL wrk_alloc( jpi, jpj, localviscsponge ) 
     42      CALL wrk_alloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab ) 
    5243 
    5344      timecoeff = REAL(Agrif_NbStepint(),wp)/Agrif_rhot() 
     
    155146      ENDDO 
    156147  
    157       IF( wrk_not_released(4, 1,2,3,4) .OR. wrk_not_released(2, 1) ) THEN 
    158          CALL ctl_stop('Agrif_Sponge_trc : failed to release workspace arrays.') 
    159          RETURN 
    160       ENDIF 
     148      CALL wrk_dealloc( jpi, jpj, localviscsponge ) 
     149      CALL wrk_dealloc( jpi, jpj, jpk, jptra, trbdiff, ztru, ztrv, ztab ) 
    161150 
    162151#endif 
Note: See TracChangeset for help on using the changeset viewer.