mode=0 working version; need to test mode=3

Author: simon

makefile

@@ -14,13 +14,9 @@ CLEAN = rm
 NETCDFINC = -I/aracbox/lib/netcdf/4.1.2/intel_12/include
 NETCDFLIB = -L/aracbox/lib/netcdf/4.1.2/intel_12/lib
 
-#FFLAGS = -check bounds -warn interface -g -traceback -fpe0 -mcmodel large -shared-intel -fp-model precise -align commons $(NETCDFINC)
-#FFLAGS = -O3 -mcmodel large -shared-intel -fp-model precise -align commons $(NETCDFINC)
 #FFLAGS = -O3 -mcmodel large -shared-intel -fp-model precise -assume byterecl $(NETCDFINC)
 FFLAGS = -check bounds -warn interface -g -traceback -mcmodel large -shared-intel -fp-model precise -assume byterecl $(NETCDFINC)
 
-#LIBS = -check bounds -warn interface -g -traceback -fpe0 -mcmodel large -shared-intel -fp-model precise -align commons $(NETCDFLIB) -lnetcdf -lnetcdff -limf -lm 
-#LIBS = -O3 -mcmodel large -shared-intel -fp-model precise -align commons $(NETCDFLIB) -lnetcdf -lnetcdff -limf -lm 
 #LIBS = -O3 -mcmodel large -shared-intel -fp-model precise -assume byterecl $(NETCDFLIB) -lnetcdf -lnetcdff -limf -lm 
 LIBS = -check bounds -warn interface -g -traceback -mcmodel large -shared-intel -fp-model precise -assume byterecl $(NETCDFLIB) -lnetcdf -lnetcdff -limf -lm 
 #-----------------------------------------------------------------------
@@ -39,6 +35,7 @@ SRCDIR = pom
 #-----------------------------------------------------------------------
 OBJS = pom08_iosub.o    \
        module_time.o    \
+       interp.o         \
        gridgen.o        \
        wave.o           \
        wind.o           \
@@ -74,4 +71,4 @@ $(BIN): $(OBJS)
 # Cleaning target
 #-----------------------------------------------------------------------
 clean:
-	@rm -f *.o *.mod *.il *.exe 
+	@rm -f *.o *.mod *.il *.exe  *genmod.f90

pom/gridgen.f

@@ -109,7 +109,7 @@ subroutine gridgen
 c        h(1,j)=h(2,j)
 c        h(im,j)=h(im-1,j)
 c      enddo
-c      call areas_masks         ! obtain fsm,dum,dvm
+      call areas_masks         ! obtain fsm,dum,dvm
 C
 C     Adjust bottom topography so that cell to cell variations
 C     in h do not exceed parameter slmax:
@@ -159,7 +159,7 @@ subroutine gridgen
 C     Set initial conditions:
 C
 c!sc:vel is used for ub
-      vel=0.2e0 ! current velocity
+      vel=0.1e0 ! current velocity
 
       do k=1,kbm1
 !sc:add dz and dzz calculation
@@ -171,7 +171,7 @@ subroutine gridgen
             sb(i,j,k)=35.e0-sbias
             tclim(i,j,k)=tb(i,j,k)
             sclim(i,j,k)=sb(i,j,k)
-!            ub(i,j,k)=vel*dum(i,j)
+            ub(i,j,k)=vel*dum(i,j)
           end do
         end do
       end do
@@ -187,7 +187,7 @@ subroutine gridgen
 C
       do j=1,jm
         do i=1,im
-!          uab(i,j)=vel*dum(i,j)
+          uab(i,j)=vel*dum(i,j)
         end do
       end do
 C