Skip to content
Snippets Groups Projects
Commit e2ecdd3d authored by Victor Stepanenko's avatar Victor Stepanenko
Browse files

Control point I/O is now operational for the case when carbon_model%par /= 1

parent 57e59221
No related branches found
No related tags found
No related merge requests found
......@@ -209,11 +209,16 @@ if (allocated(hx1ml_2d)) then
call SHIFTVI(gs%M+1,dnx,dny,hy2ml_2d,work,1_iintegers)
call WRFILD_3D_LAKE(work, coords3d, recn, unitout, .false., .true., parparams, filen)
endif
if (allocated(POCL_2d)) then
call SHIFTVI(gs%M+1,dnx,dny,POCL_2d,work,1_iintegers)
call WRFILD_3D_LAKE(work, coords3d, recn, unitout, .false., .true., parparams, filen)
endif
if (allocated(POCD_2d)) then
call SHIFTVI(gs%M+1,dnx,dny,POCD_2d,work,1_iintegers)
call WRFILD_3D_LAKE(work, coords3d, recn, unitout, .false., .true., parparams, filen)
endif
deallocate(work)
if (allocated(DOC_2d)) then
coords3d%k0 = 1
coords3d%k1 = 2*(gs%M+1)
coords3d%k00 = 1
......@@ -230,6 +235,7 @@ do j = 1, 2
enddo
call WRFILD_3D_LAKE(work, coords3d, recn, unitout, .false., .true., parparams, filen)
deallocate(work)
endif
if (allocated(itherm_2d)) then
coords3d%k0 = 1
coords3d%k1 = gs%M+2
......@@ -570,11 +576,16 @@ if (allocated(hx1ml_2d)) then
call RDFILD_3D_LAKE(coords3d, recn, unitout, .false., .true., parparams, filen, work)
call SHIFTVI(gs%M+1,dnx,dny,hy2ml_2d,work,2_iintegers)
endif
if (allocated(POCL_2d)) then
call RDFILD_3D_LAKE(coords3d, recn, unitout, .false., .true., parparams, filen, work)
call SHIFTVI(gs%M+1,dnx,dny,POCL_2d,work,2_iintegers)
endif
if (allocated(POCD_2d)) then
call RDFILD_3D_LAKE(coords3d, recn, unitout, .false., .true., parparams, filen, work)
call SHIFTVI(gs%M+1,dnx,dny,POCD_2d,work,2_iintegers)
endif
deallocate(work)
if (allocated(DOC_2d)) then
coords3d%k0 = 1
coords3d%k1 = 2*(gs%M+1)
coords3d%k00 = 1
......@@ -591,6 +602,7 @@ do j = 1, 2
enddo
enddo
deallocate(work)
endif
if (allocated(itherm_2d)) then
coords3d%k0 = 1
coords3d%k1 = gs%M+2
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment