Changeset 80 in PSPA for parmelaPSPA
- Timestamp:
- Nov 15, 2012, 12:36:13 PM (12 years ago)
- Location:
- parmelaPSPA/trunk
- Files:
-
- 1 added
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
parmelaPSPA/trunk/cellfeld.f
r12 r80 1 subroutine cellfeld(file super,nc,ne,ns)1 subroutine cellfeld(fileelec,nc,ne,ns) 2 2 c---------------------------------------------------------------------- 3 c modifie glm : suppression de filesuper (var_char.h) remplace par fileelec 3 4 c 4 5 include 'param_sz.h' … … 7 8 include 'misccom.h' 8 9 include 'syscom.h' 9 include 'var_char.h'10 c include 'var_char.h' 10 11 include 'ucom.h' 11 12 c … … 15 16 c-------------------------------------------------------------------------- 16 17 c* 18 character*256 fileelec 17 19 data rcurflag/.true./ 18 20 itest=0 19 21 iorsay=0 20 22 c open (unit=nfld,file=filesuper//'.sffld',status='old') 21 nlast = index(filesuper, ' ') 22 if (nlast .eq. 0) nlast = len(filesuper) + 1 23 open (unit=nfld,file=filesuper(1:nlast-1)//'.sffld', 23 nlast = index(fileelec, ' ') 24 print *, ' cellfeld recoit : ', fileelec 25 if (nlast .eq. 0) nlast = len(fileelec) + 1 26 open (unit=nfld,file=fileelec(1:nlast-1)//'.sffld', 24 27 1 status='old') 25 28 do 100 i=nc,ne,ns -
parmelaPSPA/trunk/file.f
r12 r80 2 2 c system dependant * 3 3 c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* 4 subroutine ouvrir 4 subroutine ouvrir(lonnom, work_path) 5 6 C ********************************************** 7 C BUT : glm 8 C ---- 9 C work_path: PATH des fichiers a ouvrir 10 C lonnom : longueur de la chaine precedente 11 C ********************************************** 12 13 14 5 15 c 6 16 include 'ucom.h' 7 17 include 'debug.h' 8 18 include 'flagcom.h' 19 20 character*256 nomfic 21 character*256 work_path 22 character*256 nomf 23 integer lonnom 24 9 25 c 10 26 c------------------------- open -------------------------------------- … … 27 43 open(unit=nav,access='sequential', 28 44 * form='formatted',status='scratch') 29 open(unit=ndes1,file='parmdesz',form='formatted', 45 46 nomf = nomfic(work_path(1:lonnom),'parmdesz') 47 print *, ' parmela file : fichier desz : ', nomf 48 49 50 open(unit=ndes1,file=nomf,form='formatted', 30 51 * status='unknown') 31 52 open(unit=ndes2,file='/tmp/parmdest',form='unformatted', -
parmelaPSPA/trunk/main.c
r73 r80 1 #include <string.h> 2 #include <stdio.h> 3 #include <stdlib.h> 1 4 2 #include <stdio.h> 5 static int maxLongPath=132; 3 6 4 void parmela_( );7 void parmela_(char*, int*); 5 8 6 9 int main(int argc,char* argv[]) 7 10 { 8 parmela_(); 11 12 if ( argc != 2 ) 13 { 14 printf(" usage : parmela complete_path_of_files \n"); 15 return 0; 16 } 17 18 char workingPath[maxLongPath]; 19 20 21 int taille = strlen(argv[1]); 22 if ( taille < maxLongPath - 1) 23 { 24 strcpy(workingPath, argv[1]); 25 } 26 else 27 { 28 printf(" parmela ERROR : path name of working area too long (should be < %d \n",maxLongPath); 29 return 0; 30 } 31 32 parmela_(workingPath, &taille); 9 33 return 0; 10 34 } -
parmelaPSPA/trunk/makefile
r78 r80 2 2 #FFLAGS= -O4 -m32 -fno-second-underscore -funroll-loops -I. 3 3 #FFLAGS= -O4 -fno-second-underscore -funroll-loops -I. 4 FFLAGS= -O4 -fno-second-underscore-funroll-loops -I.4 FFLAGS= -O4 -funroll-loops -I. 5 5 6 6 #VER=macOS 7 7 VER= 8 8 9 SRC= adjust.f alpham.f bend.f bfld.f bfldlal.f blockdata.f buncher.f byfield.f card.f cathod.f cathode.f cellfeld.f cellimp.f cfield.f change.f cread.f datime.f dawfc.f ddate.f design.f distriphase.f distrir.f distrirhom.f dpout.f drift.f eh.f eint.f elispdp.f emit90dp.f erf.f field.f fieldlal.f file.f fldring.f flds.f fsch.f gammln.f gammp.f gaus.f gcf.f geners.f genfc.f genfld.f gpindp.f graphic.f gser.f indx.f input.f interp.f interv.f irschfgt.f lodgau.f loduni.f luf.f no rran.f out1.f out2.f out3.f out4.f out5.f out6dp.f outlaldp.f output.f pardyn.f parmelamainv503.f pdist.f poisson.f prch.f prchold.f pvoldp.f quad.f ranf.f rannor.f rn32.f rndm.f rnorth.f rotate.f rtmfc.f rtrans.f sandom.f sbload.f sccfc.f scheff1.f scheff2.f scheff3.f scwall.f search.f setwall.f sextupole.f solnoid.f sort.f sortie.f spl1d1.f spl1d2.f spl2d1.f spl2d2.f spln1.f spln2.f srchfd.f sscal.f swap.f tank.f trwave.f trwfield.f trwimp.f wiggler.f xi0.f xi10x.f zout.f zoutlal.f9 SRC= adjust.f alpham.f bend.f bfld.f bfldlal.f blockdata.f buncher.f byfield.f card.f cathod.f cathode.f cellfeld.f cellimp.f cfield.f change.f cread.f datime.f dawfc.f ddate.f design.f distriphase.f distrir.f distrirhom.f dpout.f drift.f eh.f eint.f elispdp.f emit90dp.f erf.f field.f fieldlal.f file.f fldring.f flds.f fsch.f gammln.f gammp.f gaus.f gcf.f geners.f genfc.f genfld.f gpindp.f graphic.f gser.f indx.f input.f interp.f interv.f irschfgt.f lodgau.f loduni.f luf.f nomfic.f norran.f out1.f out2.f out3.f out4.f out5.f out6dp.f outlaldp.f output.f pardyn.f parmelamainv503.f pdist.f poisson.f prch.f prchold.f pvoldp.f quad.f ranf.f rannor.f rn32.f rndm.f rnorth.f rotate.f rtmfc.f rtrans.f sandom.f sbload.f sccfc.f scheff1.f scheff2.f scheff3.f scwall.f search.f setwall.f sextupole.f solnoid.f sort.f sortie.f spl1d1.f spl1d2.f spl2d1.f spl2d2.f spln1.f spln2.f srchfd.f sscal.f swap.f tank.f trwave.f trwfield.f trwimp.f wiggler.f xi0.f xi10x.f zout.f zoutlal.f 10 10 11 11 OBJ=$(SRC:.f=.o) … … 14 14 15 15 parmela : main.o $(OBJ) 16 gfortran main.o $(OBJ) -o parmela16 gfortran main.o $(OBJ) -o parmela 17 17 18 18 main.o:main.c … … 181 181 $(FF) -c $(FFLAGS) luf.f 182 182 183 nomfic.o:nomfic.f 184 $(FF) -c $(FFLAGS) nomfic.f 185 183 186 norran.o:norran.f 184 187 $(FF) -c $(FFLAGS) norran.f -
parmelaPSPA/trunk/parmelamainv503.f
r73 r80 233 233 c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* 234 234 c----------------------------------------------------------------------- 235 236 c modif F. Touze, G. Le Meur 15/11/2012 235 237 c---->program parmela 236 subroutine parmela( )238 subroutine parmela(work_path, lonnom) 237 239 c----------------------------------------------------------------------- 238 240 c label n --> n*100 … … 259 261 include 'ucom.h' 260 262 include 'upawcom.h' 263 264 character*132 work_path 265 character*256 nom_input 266 character*256 nomfic 267 character*256 fileelec 268 character*256 ficpoiss 269 integer lonnom 270 261 271 character*25 kvers 262 272 common/aswap/acor(7,imaa) … … 283 293 common /wtstep/wtstep,nstep ! mcd 284 294 c-------------------------------------------------------------------------- 295 285 296 character Istart_date*9,Iend_date*9,Istart_time*8,Iend_time*8 286 297 real*4 time_begin … … 292 303 293 304 write(*,*) ' Go Parmela' 305 306 294 307 c_for ibm 295 308 c time_begin = secnds(0.0) … … 345 358 c--- 346 359 c--- 347 open(unit=nin,file='parmin',status='old', 360 nom_input = nomfic(work_path(1:lonnom),'parmin') 361 print *, ' parmela : fichier input : ', nom_input 362 open(unit=nin,file=nom_input,status='old', 348 363 * access='sequential',form='formatted') 364 365 366 c$$$ open(unit=nin,file='parmin',status='old', 367 c$$$ * access='sequential',form='formatted') 349 368 open(unit=nnout,file='parmdebug',status='unknown', 350 369 * access='sequential',form='formatted') 351 call ouvrir 370 call ouvrir(lonnom, work_path) 352 371 c--- 353 372 call setwall … … 821 840 C modif glm 29/10/2012 read(nin,'(a12)')filesuper 822 841 read(nin,*) filesuper 842 print *, ' j''ai lu le fichier champ : ', filesuper 823 843 write(nnout,'(1x,a12)')filesuper 824 844 if(nn.lt.1)then … … 839 859 ns=vv(3) 840 860 endif 841 call cellfeld(filesuper,nc,ne,ns) 861 862 c$$$ call cellfeld(filesuper,nc,ne,ns) 863 fileelec = nomfic(work_path(1:lonnom),filesuper) 864 print *, ' j''envoie lu le fichier champ : ', fileelec 865 call cellfeld(fileelec,nc,ne,ns) ! modif glm 16/11/2012 842 866 go to 99 843 867 c---after design … … 961 985 rmult=1. 962 986 endif 963 call poisson(filepois,zoffset,rmult) 987 ficpoiss = nomfic(work_path(1:lonnom),filepois) 988 call poisson(ficpoiss,zoffset,rmult) 964 989 go to 99 965 990 c -
parmelaPSPA/trunk/poisson.f
r43 r80 7 7 include 'ucom.h' 8 8 c 9 character filepois* 129 character filepois*256 ! dimension modifiee glm (avant : 12) 10 10 c----------------------------------------------------------------------- 11 print *, ' poisson recoit : ', filepois 11 12 nlast = index(filepois, ' ') 12 13 if (nlast .eq. 0) nlast = len(filepois) + 1
Note: See TracChangeset
for help on using the changeset viewer.