Changeset 473 for trunk/NEMO/OPA_SRC/SBC/tau_forced_monthly.h90

Timestamp:

2006-05-11T17:04:37+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_060: SM: IOM + 301 levels + CORE + begining of ctl_stop

File:

• trunk/NEMO/OPA_SRC/SBC/tau_forced_monthly.h90 (modified) (6 diffs)

trunk/NEMO/OPA_SRC/SBC/tau_forced_monthly.h90

@@ -9 +9 @@
    !! * local modules variables
    INTEGER ::   &
-      numtau,   &  ! logical unit for the i-component of the wind data
-      numtav,   &  ! logical unit for the j-component of the wind data
+      numtau,   &   ! logical unit for the i-component of the wind data
+      numtav,   &   ! logical unit for the j-component of the wind data
       ntau1, ntau2  ! index of the first and second record used
-
-   CHARACTER (len=34) ::   &      !!! * monthly climatology/interanual fields
-      cl_taux,  & ! generic name of the i-component monthly NetCDF file
-      cl_tauy     ! generic name of the j-component monthly NetCDF file
-
    REAL(wp), DIMENSION(jpi,jpj,2) ::   &
       taux_dta,    &  ! i- and j-components of the surface stress (Pascal)
@@ -56 +51 @@
       !!----------------------------------------------------------------------
       !! * Modules used
-      USE ioipsl       ! NetCDF library
+      USE iom      
       !! * Arguments
       INTEGER, INTENT( in  ) ::   kt   ! ocean time step
 
       !! * Local declarations
-      INTEGER, PARAMETER ::   jpmonth = 12
-      INTEGER ::   &
-         imois, iman, itime,   &
-         i15,   &
-         ipi, ipj, ipk
-      INTEGER, DIMENSION(jpmonth) ::   istep
-      REAL(wp) , DIMENSION(jpi,jpj)::  &
-         zlon  , &
-         zlat
-      REAL(wp) , DIMENSION(jpk)::  &
-         zlev
-      REAL(wp) ::   &
-         zsecond,   & ! ???
-         zdate0,    & ! ???
-         zxy          ! coefficient of the linear time interpolation
+      INTEGER :: imois, iman, i15
+      REAL(wp) :: zxy          ! coefficient of the linear time interpolation
       !!---------------------------------------------------------------------
-      cl_taux = 'taux_1m.nc'
-      cl_tauy = 'tauy_1m.nc'
 
       ! -------------- !
@@ -86 +66 @@
       ! iman=number of dates in data file (12 for a year of monthly values)
       iman  = INT( raamo )
-      itime = jpmonth
-      ipi   = jpiglo
-      ipj   = jpjglo
-      ipk   = jpk
-
       i15 = INT( 2*FLOAT( nday ) / ( FLOAT( nobis(nmonth) ) + 0.5 ) )
-
       imois = nmonth + i15 - 1
       IF( imois == 0 ) imois = iman
-
 
       ! -------------------- !
@@ -102 +75 @@
 
       IF( kt == nit000 ) THEN
-         ntau1 = 0
-         IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*) ' tau    : MONTHLY climatological wind stress (NetCDF files)'
-         IF(lwp) WRITE(numout,*) ' ~~~    '
          
-         ! title, dimensions and tests
-
-#if defined key_agrif
-      if ( .NOT. Agrif_Root() ) then
-         cl_taux = TRIM(Agrif_CFixed())//'_'//TRIM(cl_taux)
-      endif
-#endif
+         ntau1 = 0   ! initialization
+         IF(lwp) THEN
+            WRITE(numout,*)
+            WRITE(numout,*) ' tau    : MONTHLY climatological wind stress (NetCDF files)'
+         ENDIF
+         CALL iom_open ( 'taux_1m.nc', numtau )
+         CALL iom_open ( 'tauy_1m.nc', numtav )
          
-         CALL flinopen( cl_taux, mig(1), nlci, mjg(1), nlcj,   &   ! taux on U-grid
-                        .FALSE., ipi   , ipj, ipk   ,        &
-                        zlon , zlat  , zlev   , itime,       &
-                        istep, zdate0, zsecond, numtau )
-         
-         IF( itime /= jpmonth ) THEN
-            IF(lwp) WRITE(numout,cform_err)
-            IF(lwp) WRITE(numout,*) '   problem with time coordinates in file ', cl_taux
-            IF(lwp) WRITE(numout,*) '   itime = ', itime,' jpmonth = ',jpmonth
-            nstop = nstop + 1
-         ENDIF
-         IF( ipi /= jpidta .AND. ipj /= jpjdta .AND. ipk /= 1 ) THEN
-            IF(lwp) WRITE(numout,cform_err)
-            IF(lwp) WRITE(numout,*) '   problem with size read in file ', cl_taux
-            IF(lwp) WRITE(numout,*) '   ipi = ',ipi,' jpidta = ',jpidta
-            IF(lwp) WRITE(numout,*) '   ipj = ',ipj,' jpjdta = ',jpjdta
-            IF(lwp) WRITE(numout,*) '   ipk = ',ipk,' must be 1'
-            nstop = nstop + 1
-         ENDIF
-#if defined key_agrif
-      if ( .NOT. Agrif_Root() ) then
-         cl_tauy = TRIM(Agrif_CFixed())//'_'//TRIM(cl_tauy)
-      endif
-#endif
-         CALL flinopen( cl_tauy, mig(1), nlci, mjg(1), nlcj,   &   ! tauy on V-grid
-                        .FALSE., ipi   , ipj, ipk   ,        &
-                        zlon , zlat  , zlev   , itime,       &
-                        istep, zdate0, zsecond, numtav )
-
-         IF( itime /= jpmonth ) THEN          
-            IF(lwp) WRITE(numout,cform_err)
-            IF(lwp) WRITE(numout,*) '   problem with time coordinates in file ', cl_tauy
-            IF(lwp) WRITE(numout,*) '   itime = ', itime,' jpmonth = ',jpmonth
-            nstop = nstop + 1
-         ENDIF
-         IF( ipi /= jpidta .AND. ipj /= jpjdta .AND. ipk /= 1) THEN
-            IF(lwp) WRITE(numout,cform_err)
-            IF(lwp) WRITE(numout,*) '   problem with size read in file ', cl_tauy
-            IF(lwp) WRITE(numout,*) '   ipi = ',ipi,' jpidta = ',jpidta
-            IF(lwp) WRITE(numout,*) '   ipj = ',ipj,' jpjdta = ',jpjdta
-            IF(lwp) WRITE(numout,*) '   ipk = ',ipk,' must be 1'
-            nstop = nstop + 1
-         ENDIF
       ENDIF
       
@@ -178 +104 @@
          ! Read the corresponding 2 monthly stress data
          ! ntau1
-         CALL flinget( numtau,'sozotaux',    &               ! i-component at U-pt
-            jpidta,jpjdta,1,jpmonth,ntau1,   &
-            ntau1,mig(1),nlci,mjg(1),nlcj,taux_dta(1:nlci,1:nlcj,1) )
-         CALL flinget( numtav,'sometauy',    &               ! j-component at V-pt
-            jpidta,jpjdta,1,jpmonth,ntau1,   &
-            ntau1,mig(1),nlci,mjg(1),nlcj,tauy_dta(1:nlci,1:nlcj,1) )
-         ! ntau2
-         CALL flinget( numtau,'sozotaux',    &               ! i-component at U-pt
-            jpidta,jpjdta,1,jpmonth,ntau2,   &
-            ntau2,mig(1),nlci,mjg(1),nlcj,taux_dta(1:nlci,1:nlcj,2) )
-         CALL flinget( numtav,'sometauy',    &               ! j-component at V-pt
-            jpidta,jpjdta,1,jpmonth,ntau2,   &
-            ntau2,mig(1),nlci,mjg(1),nlcj,tauy_dta(1:nlci,1:nlcj,2) )
+         CALL iom_get ( numtau, jpdom_data, 'sozotaux', taux_dta(:,:,1), ntau1 )
+         CALL iom_get ( numtav, jpdom_data, 'sometauy', tauy_dta(:,:,1), ntau1 )
+
+         CALL iom_get ( numtau, jpdom_data, 'sozotaux', taux_dta(:,:,2), ntau2 )
+         CALL iom_get ( numtav, jpdom_data, 'sometauy', tauy_dta(:,:,2), ntau2 )
          
          IF(lwp .AND. nitend-nit000 <= 100 ) THEN
@@ -230 +148 @@
       ! Closing of the 2 files (required in mpp)
       IF( kt == nitend ) THEN
-          CALL flinclo(numtau)
-          CALL flinclo(numtav)
+          CALL iom_close(numtau)
+          CALL iom_close(numtav)
       ENDIF