[cdftools] 57/228: JMM cosmetics in cdfspice.f90
Alastair McKinstry
mckinstry at moszumanska.debian.org
Fri Jun 12 08:21:29 UTC 2015
This is an automated email from the git hooks/post-receive script.
mckinstry pushed a commit to branch master
in repository cdftools.
commit 3e83b1be052ac50a8b32beda9a07aed54c77f8ef
Author: molines <molines at 1055176f-818a-41d9-83e1-73fbe5b947c5>
Date: Wed Jun 16 10:05:21 2010 +0000
JMM cosmetics in cdfspice.f90
git-svn-id: http://servforge.legi.grenoble-inp.fr/svn/CDFTOOLS/trunk@333 1055176f-818a-41d9-83e1-73fbe5b947c5
---
cdfspice.f90 | 50 +++++++++++++++++++++++++-------------------------
1 file changed, 25 insertions(+), 25 deletions(-)
diff --git a/cdfspice.f90 b/cdfspice.f90
index 17ae9bd..313d39d 100644
--- a/cdfspice.f90
+++ b/cdfspice.f90
@@ -103,31 +103,31 @@ PROGRAM cdfspice
! Compute spiciness
DO jt=1,npt
- PRINT *,' TIME = ', jt, tim(jt)/86400.,' days'
- DO jk = 1, npk
- zmask(:,:)=1.
-
- ztemp(:,:)= getvar(cfilet, 'votemper', jk ,npiglo, npjglo,ktime=jt)
- zsal(:,:) = getvar(cfilet, 'vosaline', jk ,npiglo, npjglo,ktime=jt)
-
- WHERE(zsal == 0 ) zmask = 0
-
- ! spiciness at time jt, at level jk
- zspi(:,:) = 0
- zsalref(:,:) = zsal(:,:) - 35.
- ztempt(:,:) = 1.
- DO ji=1,6
- zsalt(:,:) = 1.
- DO jj=1,5
- zspi(:,:) = zspi(:,:) + beta(ji,jj)*ztempt(:,:)*zsalt(:,:)
- zsalt(:,:) = zsalt(:,:)*zsalref(:,:)
- END DO
- ztempt(:,:) = ztempt(:,:)*ztemp(:,:)
- END DO
-
- ierr = putvar(ncout, id_varout(1) ,REAL(zspi)*zmask, jk,npiglo, npjglo,ktime=jt)
-
- END DO ! loop to next level
+ PRINT *,' TIME = ', jt, tim(jt)/86400.,' days'
+ DO jk = 1, npk
+ zmask(:,:)=1.
+
+ ztemp(:,:)= getvar(cfilet, 'votemper', jk ,npiglo, npjglo,ktime=jt)
+ zsal(:,:) = getvar(cfilet, 'vosaline', jk ,npiglo, npjglo,ktime=jt)
+
+ WHERE(zsal == 0 ) zmask = 0
+
+ ! spiciness at time jt, at level jk
+ zspi(:,:) = 0
+ zsalref(:,:) = zsal(:,:) - 35.
+ ztempt(:,:) = 1.
+ DO ji=1,6
+ zsalt(:,:) = 1.
+ DO jj=1,5
+ zspi(:,:) = zspi(:,:) + beta(ji,jj)*ztempt(:,:)*zsalt(:,:)
+ zsalt(:,:) = zsalt(:,:)*zsalref(:,:)
+ END DO
+ ztempt(:,:) = ztempt(:,:)*ztemp(:,:)
+ END DO
+
+ ierr = putvar(ncout, id_varout(1) ,REAL(zspi)*zmask, jk,npiglo, npjglo,ktime=jt)
+
+ END DO ! loop to next level
END DO ! next time frame
istatus = closeout(ncout)
--
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-science/packages/cdftools.git
More information about the debian-science-commits
mailing list