Changeset 3335 for branches/2012/dev_r3322_NOCS09_SAS/NEMOGCM/NEMO/SAS_SRC/sbcsas.F90

Timestamp:

2012-03-30T09:33:25+02:00 (12 years ago)

Author:

sga

Message:

Branch dev_r3322_NOCS09_SAS: remove interior ocean allocation

add allocation of surface fields used by LIM codes.

File:

• branches/2012/dev_r3322_NOCS09_SAS/NEMOGCM/NEMO/SAS_SRC/sbcsas.F90 (modified) (3 diffs)

branches/2012/dev_r3322_NOCS09_SAS/NEMOGCM/NEMO/SAS_SRC/sbcsas.F90

@@ -98 +98 @@
       ssu_m(:,:) = sf_sas(jf_uwd)%fnow(:,:,1) * umask(:,:,1)    ! u-velocity
       ssv_m(:,:) = sf_sas(jf_vwd)%fnow(:,:,1) * vmask(:,:,1)    ! v-velocity 
+      !
+      tsn(:,:,1,jp_tem) = sst_m(:,:)
+      tsn(:,:,1,jp_sal) = sss_m(:,:)
+      ub (:,:,1       ) = ssu_m(:,:)
+      vb (:,:,1       ) = ssv_m(:,:)
 
       IF(ln_ctl) THEN                  ! print control
@@ -123 +128 @@
       INTEGER  :: ifpr                               ! dummy loop indice
       INTEGER  :: jfld                               ! dummy loop arguments
-      INTEGER  :: inum, idv, idimv                   ! local integer
+      INTEGER  :: inum, idv, idimv, jpm              ! local integer
       !!
       CHARACTER(len=100)            ::  cn_dir   !   Root directory for location of core files
@@ -177 +182 @@
       CALL fld_fill( sf_sas, slf_d, cn_dir, 'sbc_sas_init', 'Data in file', 'namsbc_sas' )
       !
+      ! lim code currently uses surface temperature and salinity in tsn array for initialisation
+      ! and ub, vb arrays in ice dynamics
+      ! so allocate enough of arrays to use
+      !
+      jpm = MAX(jp_tem, jp_sal)
+      ALLOCATE( tsn(jpi,jpj,1,jpm), STAT=ierr0 )
+      ALLOCATE( ub(jpi,jpj,1)     , STAT=ierr1 )
+      ALLOCATE( vb(jpi,jpj,1)     , STAT=ierr2 )
+      ierr = ierr0 + ierr1 + ierr2
+      IF( ierr > 0 ) THEN
+         CALL ctl_stop('sbc_sas_init: unable to allocate surface arrays')
+      ENDIF
+      !
    END SUBROUTINE sbc_sas_init