Changeset 473 for trunk/NEMO/OPA_SRC/TRA/trabbc.F90

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/TRA/trabbc.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/TRA/trabbc.F90

@@ -137 +137 @@
       !!----------------------------------------------------------------------
       !! * Modules used
-      USE ioipsl
+      USE iom
 
       !! * local declarations
-      CHARACTER (len=32) ::   clname
       INTEGER  ::   ji, jj              ! dummy loop indices
-      INTEGER  ::   inum = 11           ! temporary logical unit
-      INTEGER  ::   itime               ! temporary integers
-      REAL(wp) ::   zdate0, zdt         ! temporary scalars
-      REAL(wp), DIMENSION(1) :: zdept   ! temporary workspace
-      REAL(wp), DIMENSION(jpidta,jpjdta) ::   &
-         zlamt, zphit, zdta   ! temporary workspace
+      INTEGER  ::   inum                ! temporary logical unit
 
       NAMELIST/nambbc/ngeo_flux, ngeo_flux_const 
@@ -188 +182 @@
       CASE ( 2 )                ! variable geothermal heat flux
          ! read the geothermal fluxes in mW/m2
-         clname = 'geothermal_heating'
-#if defined key_agrif
-      if ( .NOT. Agrif_Root() ) then
-         clname = TRIM(Agrif_CFixed())//'_'//TRIM(clname)
-      endif
-#endif    
-         itime = 1
-         zlamt(:,:) = 0.
-         zphit(:,:) = 0.
-         IF(lwp) WRITE(numout,*) '             *** variable geothermal heat flux read in ', clname, ' file'
-         CALL restini( clname, jpidta, jpjdta, zlamt, zphit, 1, zdept , 'NONE',   &
-            &          itime, zdate0, zdt, inum, domain_id=nidom )
-         CALL restget( inum, 'heatflow', jpidta, jpjdta, 1, itime, .FALSE., zdta )
-         DO jj = 1, nlcj
-            DO ji = 1, nlci
-              qgh_trd(ji,jj) = zdta(mig(ji),mjg(jj))
-            END DO
-         END DO
-
-         CALL restclo( inum )
+
+         IF(lwp) WRITE(numout,*) '             *** variable geothermal heat flux'
+         CALL iom_open ( 'geothermal_heating.nc', inum )
+         CALL iom_get ( inum, jpdom_data, 'heatflow', qgh_trd )
+         CALL iom_close (inum)
+
          qgh_trd(:,:) = qgh_trd(:,:) * 1.e-3 ! conversion in W/m2
 
       CASE DEFAULT
-         IF(lwp) WRITE(numout,cform_err)
-         IF(lwp) WRITE(numout,*) '     bad flag value for ngeo_flux = ', ngeo_flux
-         nstop = nstop + 1
-
+         WRITE(ctmp1,*) '     bad flag value for ngeo_flux = ', ngeo_flux
+         CALL ctl_stop( ctmp1 )
       END SELECT