/[MITgcm]/MITgcm_contrib/mlosch/optim_m1qn3/testbed/driver.F
ViewVC logotype

Diff of /MITgcm_contrib/mlosch/optim_m1qn3/testbed/driver.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by mlosch, Wed May 2 12:27:42 2012 UTC revision 1.2 by mlosch, Mon May 7 12:09:27 2012 UTC
# Line 77  c     first guess of x,f,g (izs=-1) Line 77  c     first guess of x,f,g (izs=-1)
77          call m1qn3 (simul_rc,euclid,ctonbe,ctcabe,n,x,f,g,dxmin,df1,          call m1qn3 (simul_rc,euclid,ctonbe,ctcabe,n,x,f,g,dxmin,df1,
78       &       epsg,normtype,impres,io,imode,omode,niter,nsim,       &       epsg,normtype,impres,io,imode,omode,niter,nsim,
79       &       iz,dz,ndz,reverse,indic,izs,rzs,dzs)       &       iz,dz,ndz,reverse,indic,izs,rzs,dzs)
80          write(fmt,'(A,I2,A)') '(A,I4,I8,',n,'E12.4)'          write(fmt,'(A,I4,A)') '(A,I4,I8,',n,'E12.4)'
81          write(*,fmt) 'ml-driver iter = ',iiter,reverse, (1.-x(i),i=1,n)          write(*,fmt) 'ml-driver iter = ',iiter,reverse, (1.-x(i),i=1,n)
82          call mysimul(indic,n,f,x,g,izs,rzs,dzs)          call mysimul(indic,n,f,x,g,izs,rzs,dzs)
83          write(fmt,'(A,I2,A)') '(A,',n+1,'E12.4)'          write(fmt,'(A,I4,A)') '(A,',n+1,'E12.4)'
84          write(*,fmt) 'ml-driver f    = ', f, (g(i),i=1,n)          write(*,fmt) 'ml-driver f    = ', f, (g(i),i=1,n)
85  c        print *, 'ml-driver: niter = ', niter, nsim  c        print *, 'ml-driver: niter = ', niter, nsim
86          if ( reverse .lt. 0 ) goto 100          if ( reverse .lt. 0 ) goto 100
# Line 99  c Line 99  c
99         call m1qn3 (simul_rc,euclid,ctonbe,ctcabe,n,x,f,g,dxmin,df1,         call m1qn3 (simul_rc,euclid,ctonbe,ctcabe,n,x,f,g,dxmin,df1,
100       &      epsg,normtype,impres,io,imode,omode,niter,nsim,       &      epsg,normtype,impres,io,imode,omode,niter,nsim,
101       &      iz,dz,ndz,reverse,indic,izs,rzs,dzs)       &      iz,dz,ndz,reverse,indic,izs,rzs,dzs)
102         write(fmt,'(A,I2,A)') '(A,I4,I8,',n,'E12.4)'         write(fmt,'(A,I4,A)') '(A,I4,I8,',n,'E12.4)'
103         write(*,fmt) 'ml-driver iter = ',iiter,reverse, (1.-x(i),i=1,n)         write(*,fmt) 'ml-driver iter = ',iiter,reverse, (1.-x(i),i=1,n)
104         call mysimul(indic,n,f,x,g,izs,rzs,dzs)         call mysimul(indic,n,f,x,g,izs,rzs,dzs)
105         write(fmt,'(A,I2,A)') '(A,',n+1,'E12.4)'         write(fmt,'(A,I4,A)') '(A,',n+1,'E12.4)'
106         write(*,fmt) 'ml-driver f    = ', f, (g(i),i=1,n)         write(*,fmt) 'ml-driver f    = ', f, (g(i),i=1,n)
107  c       print *, 'ml-driver niter =', niter, nsim  c       print *, 'ml-driver niter =', niter, nsim
108        enddo        enddo
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22