Changeset 631 for trunk/NEMO/OPA_SRC/DIA/diadimg.F90
- Timestamp:
- 2007-03-02T17:40:54+01:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DIA/diadimg.F90
r460 r631 64 64 REAL(sp),DIMENSION(jpk) :: z4dep ! vertical level (sp) 65 65 66 CHARACTER(LEN=4) :: clver='@!01' 66 CHARACTER(LEN=80) :: clname ! name of file in case of dimgnnn 67 CHARACTER(LEN=4) :: clver='@!01' ! dimg string identifier 67 68 !!--------------------------------------------------------------------------- 68 69 … … 101 102 END SELECT 102 103 103 !! * Open file 104 OPEN (inum, FILE=cd_name, FORM='UNFORMATTED', ACCESS='DIRECT', RECL=irecl4 ) 104 IF ( ln_dimgnnn ) THEN 105 WRITE(clname,'(a,a,i3.3)') TRIM(cd_name),'.',narea 106 OPEN (inum, FILE=clname, FORM='UNFORMATTED', ACCESS='DIRECT', RECL=irecl4 ) 107 WRITE(inum,REC=1 ) clver, cd_text, irecl4, & 108 & jpi,jpj, klev, 1 , 1 , & 109 & zwest, zsouth, zdx, zdy, zspval, & 110 & z4dep(1:klev), & 111 & ztimm, & 112 & narea, jpnij,jpiglo,jpjglo,jpizoom, jpjzoom, & ! extension to dimg for mpp output 113 & nlcit,nlcjt, nldit, nldjt, nleit, nlejt, nimppt, njmppt ! 105 114 106 !! * Write header on record #1 107 IF(lwp) WRITE(inum,REC=1 ) clver, cd_text, irecl4, & 108 & jpi,jpj, klev*jpnij, 1 , 1 , & 109 & zwest, zsouth, zdx, zdy, zspval, & 110 & (z4dep(1:klev),jn=1,jpnij), & 111 & ztimm, & 112 & narea, jpnij,jpiglo,jpjglo,jpizoom, jpjzoom, & ! extension to dimg for mpp output 113 & nlcit,nlcjt, nldit, nldjt, nleit, nlejt, nimppt, njmppt ! 115 !! * Write klev levels 116 IF ( cd_type == 'I' ) THEN 114 117 115 !! * Write klev levels 116 IF ( cd_type == 'I' ) THEN 118 DO jk = 1, klev 119 irec =1 + jk 120 z42d(:,:) = ptab(:,:,ksubi(jk)) 121 WRITE(inum,REC=irec) z42d(:,:) 122 END DO 123 ELSE 124 DO jk = 1, klev 125 irec =1 + jk 126 z42d(:,:) = ptab(:,:,jk) 127 WRITE(inum,REC=irec) z42d(:,:) 128 END DO 129 ENDIF 130 ELSE 131 !! Standard dimgproc (1 file per variable, all procs. write to this file ) 132 !! * Open file 133 OPEN (inum, FILE=cd_name, FORM='UNFORMATTED', ACCESS='DIRECT', RECL=irecl4 ) 117 134 118 DO jk = 1, klev 119 irec =1 + klev * (narea -1) + jk 120 z42d(:,:) = ptab(:,:,ksubi(jk)) 121 WRITE(inum,REC=irec) z42d(:,:) 122 END DO 123 ELSE 124 DO jk = 1, klev 125 irec =1 + klev * (narea -1) + jk 126 z42d(:,:) = ptab(:,:,jk) 127 WRITE(inum,REC=irec) z42d(:,:) 128 END DO 135 !! * Write header on record #1 136 IF(lwp) WRITE(inum,REC=1 ) clver, cd_text, irecl4, & 137 & jpi,jpj, klev*jpnij, 1 , 1 , & 138 & zwest, zsouth, zdx, zdy, zspval, & 139 & (z4dep(1:klev),jn=1,jpnij), & 140 & ztimm, & 141 & narea, jpnij,jpiglo,jpjglo,jpizoom, jpjzoom, & ! extension to dimg for mpp output 142 & nlcit,nlcjt, nldit, nldjt, nleit, nlejt, nimppt, njmppt ! 143 144 !! * Write klev levels 145 IF ( cd_type == 'I' ) THEN 146 147 DO jk = 1, klev 148 irec =1 + klev * (narea -1) + jk 149 z42d(:,:) = ptab(:,:,ksubi(jk)) 150 WRITE(inum,REC=irec) z42d(:,:) 151 END DO 152 ELSE 153 DO jk = 1, klev 154 irec =1 + klev * (narea -1) + jk 155 z42d(:,:) = ptab(:,:,jk) 156 WRITE(inum,REC=irec) z42d(:,:) 157 END DO 158 ENDIF 129 159 ENDIF 130 160
Note: See TracChangeset
for help on using the changeset viewer.