| [3184] | 1 | SUBROUTINE writeg1d(ngrid,nx,x,nom,titre) |
|---|
| 2 | USE comg1d_mod, ONLY: g1d_nomfich,g1d_premier,g1d_unitfich, |
|---|
| 3 | & g1d_irec,g1d_nvar,g1d_nomvar,g1d_titrevar,g1d_dimvar,g1d_nlayer |
|---|
| 4 | IMPLICIT NONE |
|---|
| 5 | |
|---|
| 6 | c....................................................................... |
|---|
| 7 | c |
|---|
| 8 | c ecriture de x pour GRADS-1D |
|---|
| 9 | c |
|---|
| 10 | c in : |
|---|
| 11 | c * ngrid ---> pour controler que l'on est bien en 1D |
|---|
| 12 | c * nx ---> taille du vecteur a stocker |
|---|
| 13 | c "1" pour une variable de surface |
|---|
| 14 | c "nlayer" pour une variable de centre de couche |
|---|
| 15 | c "nlayer+1" pour une variable d'interface |
|---|
| 16 | c * x ---> variable a stocker |
|---|
| 17 | c * nom ---> nom "pour grads" |
|---|
| 18 | c * titre ---> titre "pour grads" |
|---|
| 19 | c |
|---|
| 20 | c....................................................................... |
|---|
| 21 | c |
|---|
| 22 | |
|---|
| 23 | c |
|---|
| 24 | c....................................................................... |
|---|
| 25 | c declaration des arguments |
|---|
| 26 | c |
|---|
| 27 | INTEGER ngrid,nx,i |
|---|
| 28 | REAL*4 xr4(1000) |
|---|
| 29 | REAL x(nx) |
|---|
| 30 | CHARACTER*(*) nom,titre |
|---|
| 31 | c |
|---|
| 32 | c declaration des arguments |
|---|
| 33 | c....................................................................... |
|---|
| 34 | c declaration des variables locales |
|---|
| 35 | c |
|---|
| 36 | INTEGER ilayer,ivar |
|---|
| 37 | LOGICAL test |
|---|
| 38 | c |
|---|
| 39 | c declaration des variables locales |
|---|
| 40 | c....................................................................... |
|---|
| 41 | c controle 1D |
|---|
| 42 | c |
|---|
| 43 | c print*,'ngrid=',ngrid |
|---|
| 44 | IF (ngrid.NE.1) return |
|---|
| 45 | c |
|---|
| 46 | c controle 1D |
|---|
| 47 | c....................................................................... |
|---|
| 48 | c copy => force en reel 4 pour l'ecriture dans grads1d.dat |
|---|
| 49 | |
|---|
| 50 | do i=1,nx |
|---|
| 51 | xr4(i) = x(i) |
|---|
| 52 | enddo |
|---|
| 53 | |
|---|
| 54 | c copy => force en reel 4 pour l'ecriture dans grads1d.dat |
|---|
| 55 | c....................................................................... |
|---|
| 56 | c ouverture du fichier au premier appel |
|---|
| 57 | |
|---|
| 58 | |
|---|
| 59 | g1d_nomfich='g1d.dat' |
|---|
| 60 | |
|---|
| 61 | IF (g1d_premier) THEN |
|---|
| 62 | OPEN (g1d_unitfich,FILE=g1d_nomfich |
|---|
| 63 | & ,FORM='unformatted',ACCESS='direct',RECL=4) |
|---|
| 64 | g1d_irec=0 |
|---|
| 65 | g1d_nvar=0 |
|---|
| 66 | g1d_premier=.false. |
|---|
| 67 | ENDIF |
|---|
| 68 | |
|---|
| 69 | c ouverture du fichier au premier appel |
|---|
| 70 | c....................................................................... |
|---|
| 71 | c pour l'ecriture du fichier ctl |
|---|
| 72 | |
|---|
| 73 | test=.true. |
|---|
| 74 | DO ivar=1,g1d_nvar |
|---|
| 75 | IF (nom.EQ.g1d_nomvar(ivar)) test=.false. |
|---|
| 76 | IF (nx .GT. 1000) then |
|---|
| 77 | print*,'ERROR: nx > 1000 dans writeg1d.F' |
|---|
| 78 | print*,'Changer la dimension de xr4' |
|---|
| 79 | call exit(1) |
|---|
| 80 | ENDIF |
|---|
| 81 | ENDDO |
|---|
| 82 | IF (test) THEN |
|---|
| 83 | g1d_nvar=g1d_nvar+1 |
|---|
| 84 | g1d_nomvar(g1d_nvar)=nom |
|---|
| 85 | g1d_titrevar(g1d_nvar)=titre |
|---|
| 86 | IF (nx.EQ.1) THEN |
|---|
| 87 | g1d_dimvar(g1d_nvar)=0 |
|---|
| 88 | ELSEIF (nx.EQ.g1d_nlayer) THEN |
|---|
| 89 | g1d_dimvar(g1d_nvar)=g1d_nlayer |
|---|
| 90 | ELSEIF (nx.EQ.g1d_nlayer+1) THEN |
|---|
| 91 | g1d_dimvar(g1d_nvar)=g1d_nlayer+1 |
|---|
| 92 | ELSE |
|---|
| 93 | PRINT *,'._. probleme de dimension dans GRADS-1D ._.' |
|---|
| 94 | print*,'NX = ',nx |
|---|
| 95 | print*,'g1d_nlayer = ',g1d_nlayer |
|---|
| 96 | ENDIF |
|---|
| 97 | ENDIF |
|---|
| 98 | c |
|---|
| 99 | c pour l'ecriture du fichier ctl |
|---|
| 100 | c....................................................................... |
|---|
| 101 | c ecriture |
|---|
| 102 | c |
|---|
| 103 | IF (nx.EQ.1) THEN |
|---|
| 104 | g1d_irec=g1d_irec+1 |
|---|
| 105 | WRITE(g1d_unitfich,REC=g1d_irec) xr4(1) |
|---|
| 106 | ELSE |
|---|
| 107 | DO ilayer=1,g1d_nlayer |
|---|
| 108 | g1d_irec=g1d_irec+1 |
|---|
| 109 | WRITE(g1d_unitfich,REC=g1d_irec) xr4(ilayer) |
|---|
| 110 | ENDDO |
|---|
| 111 | ENDIF |
|---|
| 112 | c |
|---|
| 113 | c ecriture |
|---|
| 114 | c....................................................................... |
|---|
| 115 | c |
|---|
| 116 | 10001 CONTINUE |
|---|
| 117 | c |
|---|
| 118 | c....................................................................... |
|---|
| 119 | c |
|---|
| 120 | RETURN |
|---|
| 121 | END |
|---|
| 122 | |
|---|
| 123 | |
|---|
| 124 | |
|---|
| 125 | |
|---|
| 126 | |
|---|
| 127 | c ********************************************************************* |
|---|
| 128 | c ********************************************************************* |
|---|
| 129 | |
|---|
| 130 | SUBROUTINE endg1d(ngrid,nlayer,zlayer,ndt) |
|---|
| 131 | USE time_phylmdz_mod, ONLY: dtphys, daysec |
|---|
| 132 | USE comg1d_mod, ONLY: g1d_nomfich,g1d_unitfich,g1d_nvar, |
|---|
| 133 | & g1d_nomvar,g1d_titrevar,g1d_dimvar,g1d_nlayer,g1d_unitctl, |
|---|
| 134 | & g1d_nomctl,saveG1D |
|---|
| 135 | IMPLICIT NONE |
|---|
| 136 | c....................................................................... |
|---|
| 137 | c |
|---|
| 138 | c ecriture du fichier de controle pour GRADS-1D |
|---|
| 139 | c |
|---|
| 140 | c in : |
|---|
| 141 | c * ngrid ---> pour controler que l'on est bien en 1D |
|---|
| 142 | c * nlayer ---> nombre de couches |
|---|
| 143 | c * zlayer ---> altitude au centre de chaque couche (km) |
|---|
| 144 | c * ndt ---> nombre de pas de temps |
|---|
| 145 | c |
|---|
| 146 | c....................................................................... |
|---|
| 147 | c |
|---|
| 148 | |
|---|
| 149 | |
|---|
| 150 | c |
|---|
| 151 | c....................................................................... |
|---|
| 152 | c declaration des arguments |
|---|
| 153 | c |
|---|
| 154 | INTEGER ngrid,nlayer |
|---|
| 155 | REAL zlayer(nlayer) |
|---|
| 156 | INTEGER ndt |
|---|
| 157 | c |
|---|
| 158 | c declaration des arguments |
|---|
| 159 | c....................................................................... |
|---|
| 160 | c declaration des variables locales |
|---|
| 161 | c |
|---|
| 162 | INTEGER ivar,ilayer |
|---|
| 163 | c |
|---|
| 164 | |
|---|
| 165 | |
|---|
| 166 | ! integer saveG1D |
|---|
| 167 | |
|---|
| 168 | c declaration des variables locales |
|---|
| 169 | c....................................................................... |
|---|
| 170 | c contole 1D |
|---|
| 171 | c |
|---|
| 172 | IF (ngrid.NE.1) GOTO 10001 |
|---|
| 173 | c |
|---|
| 174 | c contole 1D |
|---|
| 175 | c....................................................................... |
|---|
| 176 | c |
|---|
| 177 | IF (nlayer.ne.g1d_nlayer) |
|---|
| 178 | &PRINT *,'._. probleme de dimension dans GRADS-1D (endg1d.F) ' |
|---|
| 179 | c |
|---|
| 180 | c....................................................................... |
|---|
| 181 | c |
|---|
| 182 | CLOSE (g1d_unitfich) |
|---|
| 183 | c |
|---|
| 184 | c....................................................................... |
|---|
| 185 | |
|---|
| 186 | |
|---|
| 187 | OPEN (g1d_unitctl,FILE=g1d_nomctl,FORM='formatted',RECL=4*100) |
|---|
| 188 | WRITE (g1d_unitctl,'(a4,2x,a1,a20)') 'DSET','^',g1d_nomfich |
|---|
| 189 | WRITE (g1d_unitctl,'(a5,2x,a20)') 'UNDEF ','1.E+30' |
|---|
| 190 | WRITE (g1d_unitctl,'(a11)') 'FORMAT YREV' |
|---|
| 191 | WRITE (g1d_unitctl,'(a5,2x,a30)') 'TITLE ','champs 1D' |
|---|
| 192 | WRITE (g1d_unitctl,'(a5,i4,a20)') 'XDEF ',1,' LINEAR 0 1' |
|---|
| 193 | WRITE (g1d_unitctl,'(a5,i4,a20)') 'YDEF ',1,' LINEAR 0 1' |
|---|
| 194 | WRITE (g1d_unitctl,'(a5,i4,a20)') 'ZDEF ',g1d_nlayer,' LEVELS' |
|---|
| 195 | WRITE (g1d_unitctl,'(5(1x,f13.5))') |
|---|
| 196 | & (zlayer(ilayer),ilayer=1,g1d_nlayer) |
|---|
| 197 | |
|---|
| 198 | c Writing true timestep in g1d.ctl (in planet "minutes"= sol/(60*24) |
|---|
| 199 | ! ivar =min( max(int(1440.*dtphys/daysec +0.5),1) , 99) |
|---|
| 200 | ! WRITE (g1d_unitctl,'(a4,2x,i10,a19,i2,a2)') |
|---|
| 201 | ! & 'TDEF ',ndt,' LINEAR 01JAN2000 ', ivar,'MN ' |
|---|
| 202 | |
|---|
| 203 | ivar =min( max(int(1440.*dtphys*saveG1D/daysec +0.5),1) , 99) |
|---|
| 204 | ! not sure ivar is right, but it doesnt matter |
|---|
| 205 | WRITE (g1d_unitctl,'(a4,2x,i10,a19,i2,a2)') |
|---|
| 206 | & 'TDEF ',ndt/saveG1D,' LINEAR 01JAN2000 ', ivar,'MN ' |
|---|
| 207 | |
|---|
| 208 | WRITE (g1d_unitctl,'(a5,i5)') 'VARS ',g1d_nvar |
|---|
| 209 | DO ivar=1,g1d_nvar |
|---|
| 210 | WRITE (g1d_unitctl,'(a9,3x,i4,i3,1x,a39)') |
|---|
| 211 | & g1d_nomvar(ivar),g1d_dimvar(ivar),99,g1d_titrevar(ivar) |
|---|
| 212 | ENDDO |
|---|
| 213 | WRITE (g1d_unitctl,'(a7)') 'ENDVARS' |
|---|
| 214 | CLOSE (g1d_unitctl) |
|---|
| 215 | c |
|---|
| 216 | c....................................................................... |
|---|
| 217 | c |
|---|
| 218 | 10001 CONTINUE |
|---|
| 219 | c |
|---|
| 220 | c....................................................................... |
|---|
| 221 | c |
|---|
| 222 | RETURN |
|---|
| 223 | END |
|---|
