- Timestamp:
- Jul 24, 2024, 2:54:37 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bilan_dyn_loc.f90
r5106 r5116 41 41 ! =========== 42 42 43 integer:: ntrac44 real:: dt_app,dt_cum45 real:: ps(iip1,jjb_u:jje_u)46 real:: masse(iip1,jjb_u:jje_u,llm),pk(iip1,jjb_u:jje_u,llm)47 real:: flux_u(iip1,jjb_u:jje_u,llm)48 real:: flux_v(iip1,jjb_v:jje_v,llm)49 real:: teta(iip1,jjb_u:jje_u,llm)50 real:: phi(iip1,jjb_u:jje_u,llm)51 real:: ucov(iip1,jjb_u:jje_u,llm)52 real:: vcov(iip1,jjb_v:jje_v,llm)53 real:: trac(iip1,jjb_u:jje_u,llm,ntrac)43 INTEGER :: ntrac 44 REAL :: dt_app,dt_cum 45 REAL :: ps(iip1,jjb_u:jje_u) 46 REAL :: masse(iip1,jjb_u:jje_u,llm),pk(iip1,jjb_u:jje_u,llm) 47 REAL :: flux_u(iip1,jjb_u:jje_u,llm) 48 REAL :: flux_v(iip1,jjb_v:jje_v,llm) 49 REAL :: teta(iip1,jjb_u:jje_u,llm) 50 REAL :: phi(iip1,jjb_u:jje_u,llm) 51 REAL :: ucov(iip1,jjb_u:jje_u,llm) 52 REAL :: vcov(iip1,jjb_v:jje_v,llm) 53 REAL :: trac(iip1,jjb_u:jje_u,llm,ntrac) 54 54 55 55 ! Local : … … 61 61 !$OMP THREADPRIVATE(first) 62 62 63 real:: zz,zqy63 REAL :: zz,zqy 64 64 REAl,SAVE,ALLOCATABLE :: zfactv(:,:) 65 65 … … 69 69 !ym character*6 nom(nQ) 70 70 !ym character*6 unites(nQ) 71 character(len=6),save :: nom(nQ)72 character(len=6),save :: unites(nQ)73 74 character(len=10) file75 integer:: ifile71 CHARACTER(LEN=6),save :: nom(nQ) 72 CHARACTER(LEN=6),save :: unites(nQ) 73 74 CHARACTER(LEN=10) file 75 INTEGER :: ifile 76 76 parameter (ifile=4) 77 77 … … 84 84 !$OMP THREADPRIVATE(time,itau) 85 85 86 real:: ww86 REAL :: ww 87 87 88 88 ! variables dynamiques intermédiaires … … 110 110 111 111 ! champs de tansport en moyenne zonale 112 integer:: ntr,itr112 INTEGER :: ntr,itr 113 113 parameter (ntr=5) 114 114 … … 122 122 INTEGER,PARAMETER :: iave=1,itot=2,immc=3,itrs=4,istn=5 123 123 124 character(len=3) :: ctrs(ntr)124 CHARACTER(LEN=3) :: ctrs(ntr) 125 125 data ctrs/' ','TOT','MMC','TRS','STN'/ 126 126 … … 131 131 real,SAVE,ALLOCATABLE :: zv(:,:),psi(:,:) 132 132 133 integer:: i,j,l,iQ133 INTEGER :: i,j,l,iQ 134 134 135 135 … … 137 137 ! --------------------------------------------------------- 138 138 139 character(len=10) :: infile139 CHARACTER(LEN=10) :: infile 140 140 141 141 integer, save :: fileid 142 integer:: thoriid, zvertiid142 INTEGER :: thoriid, zvertiid 143 143 144 144 INTEGER,SAVE,ALLOCATABLE :: ndex3d(:) … … 146 146 ! Variables locales 147 147 ! 148 integer:: tau0149 real:: zjulian150 character(len=3) :: str151 character(len=10) :: ctrac152 integer:: ii,jj153 integer:: zan, dayref148 INTEGER :: tau0 149 REAL :: zjulian 150 CHARACTER(LEN=3) :: str 151 CHARACTER(LEN=10) :: ctrac 152 INTEGER :: ii,jj 153 INTEGER :: zan, dayref 154 154 ! 155 155 real,SAVE,ALLOCATABLE :: rlong(:),rlatg(:) 156 integer:: jjb,jje,jjn,ijb,ije156 INTEGER :: jjb,jje,jjn,ijb,ije 157 157 type(Request),SAVE :: Req 158 158 !$OMP THREADPRIVATE(Req) … … 178 178 itau=itau+1 179 179 180 if (first) then180 if (first) THEN 181 181 !$OMP BARRIER 182 182 !$OMP MASTER … … 223 223 ! ncum est la frequence de stokage en pas de temps 224 224 ncum=dt_cum/dt_app 225 if (abs(ncum*dt_app-dt_cum)>1.e-5*dt_app) then225 if (abs(ncum*dt_app-dt_cum)>1.e-5*dt_app) THEN 226 226 WRITE(lunout,*) & 227 227 'Pb : le pas de cumule doit etre multiple du pas' … … 296 296 do iQ=1,nQ 297 297 do itr=1,ntr 298 if(itr==1) then298 IF(itr==1) THEN 299 299 znom(itr,iQ)=nom(iQ) 300 300 znoml(itr,iQ)=nom(iQ) … … 404 404 !===================================================================== 405 405 ! 406 if(icum==0) then406 IF(icum==0) THEN 407 407 jjb=jj_begin 408 408 jje=jj_end … … 582 582 ! PAS DE TEMPS D'ECRITURE 583 583 !===================================================================== 584 if (icum==ncum) then584 if (icum==ncum) THEN 585 585 !===================================================================== 586 586 … … 776 776 ! sorties proprement dites 777 777 !$OMP MASTER 778 if (i_sortie==1) then778 if (i_sortie==1) THEN 779 779 jjb=jj_begin 780 780 jje=jj_end … … 841 841 842 842 843 end subroutinebilan_dyn_loc843 END SUBROUTINE bilan_dyn_loc
Note: See TracChangeset
for help on using the changeset viewer.