<p><b>dwj07@fsu.edu</b> 2012-04-04 15:58:42 -0600 (Wed, 04 Apr 2012)</p><p><br>
        -- BRANCH COMMIT --<br>
<br>
        Fixing configuration of idealized acc.<br>
        Some tuning still has to be done.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/basin_src/basin-template.F
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/basin_src/basin-template.F        2012-04-04 21:12:48 UTC (rev 1749)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/basin_src/basin-template.F        2012-04-04 21:58:42 UTC (rev 1750)
@@ -89,7 +89,7 @@
! Step 3: Specify some Parameters
real (kind=8), parameter :: &
- h_total_max = 2000.0, &
+ h_total_max = 5000.0, &
u_max = 0.0, &
u_src_max = 0.1, & ! max wind stress, N/m2
beta = 1.4e-11, &
@@ -1000,7 +1000,7 @@
real (kind=4), allocatable, dimension(:) :: x,y, work_kmt
real (kind=4), allocatable, dimension(:,:) :: ztopo
integer :: nx, ny, inx, iny, ix, iy, kmt_neighbor_max
-real :: pi, dtr, zdata, rlon, rlat, r, ymin, ymax, xmin, xmax
+real :: pi, dtr, zdata, rlon, rlat, r, ymin, ymax, xmin, xmax, xmid, sigma
real :: latmin, latmax, lonmin, lonmax
logical :: flag, kmt_flag
pi = 4.0*atan(1.0)
@@ -1020,23 +1020,29 @@
where(yCell.lt.1.001*ymin) kmt = 0
where(yCell.gt.0.999*ymax) kmt = 0
+ xmin = minval(xCell)
+ xmax = maxval(xCell)
+ xmid = (xmax-xmin)/2.0
+ sigma = (xmin-xmax)/10.0
+
do iCell=1,nCells
- if ( xCell(iCell).gt.1000.0e3 &
- .and.xCell(iCell).lt.1900.0e3) then
+! if ( xCell(iCell).gt.xmid-4*sigma &
+! .and.xCell(iCell).lt.xmid+4*sigma) then
- zdata = 5000.0-*RIDGE*exp(-2.0*((xCell(iCell)-1450.0e3)/350.0e3)**2)
+ zdata = 5000.0-*RIDGE*exp(-2.0*((xCell(iCell)-xmid)/sigma)**2)
r = 0
kmt_flag=.false.
do k=1,nVertLevelsMod
if(.not.kmt_flag) then
- r = r + dz(k)
+! r = r + dz(k)
+ r = r + h_total_max / nVertLevelsMOD
if(r.gt.zdata) then
kmt(iCell) = k
kmt_flag = .true.
endif
endif
enddo
- endif
+! endif
enddo
allocate(work_kmt(nCells))
Modified: branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/makeMeshes.sh
===================================================================
--- branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/makeMeshes.sh        2012-04-04 21:12:48 UTC (rev 1749)
+++ branches/ocean_projects/ocean_test_cases_staging/ocean/idealized_acc/makeMeshes.sh        2012-04-04 21:58:42 UTC (rev 1750)
@@ -5,8 +5,9 @@
## To start, grids are periodic in the x direction. This may change later                         ##
########################################################################################
-SPACINGS="80000 40000 20000 10000"
-VERTLEVS="25 40"
+SPACINGS="8000 4000 2000 1000"
+SPACINGS="8000"
+VERTLEVS="34"
RIDGESIZES="2000"
TCNAME="idealized_acc"
@@ -29,8 +30,8 @@
        PROCS="$3"
fi
-X_EXTENT="320000"
-Y_EXTENT="160000"
+X_EXTENT="3200000"
+Y_EXTENT="1600000"
############################################################
## Setup Default namelist keys and values to update later ##
</font>
</pre>