Changeset 41


Ignore:
Timestamp:
04/26/11 15:45:50 (13 years ago)
Author:
jbrlod
Message:

run ref with batch

Location:
trunk
Files:
3 added
5 deleted
6 edited

Legend:

Unmodified
Added
Removed

  • trunk/scripts/run_batch_sinobad.sh

    r40 r41  
    1 #!/bin/tcsh 
     1#!/bin/sh 
    22#PBS -e ../scripts/script_outputs/ 
    33#PBS -o ../scripts/script_outputs/ 
    … …  
    88#Command to submit the job: 
    99#qsub run_batch_sinobad.sh 
    10 cd /usr/home/jbrlod/these/postdoc/collaborations/NEMO/nemoyao/svn/trunk/src/ 
     10cd $PBS_O_WORKDIR 
     11cd ../src/ 
     12 
     13echo "****************************************"  
     14echo "*       run reference experiment       *" 
     15echo "*             (sequential)             *" 
     16echo "*              26/04/2011              *" 
     17echo "*  julien.brajard@locean-ipsl.upmc.fr  *" 
     18echo "****************************************" 
    1119 
    1220$YAODIR/etc/YaoI -p +O3 sinobad 
    13 ./sinobad -i ../scripts/sinobad.i 
     21time ./sinobad -i ../scripts/sinobad.i 
    1422 

  • trunk/src/Yworkdir/Y2sinobad.h

    r38 r41  
    38733873         Yting[0]=YS1_tb( Yi, Yj, 0, YTemps-1); 
    38743874         Ysbcgyre_qsr_qns(Yi, Yj)->forward( Yting[0] ); 
    3875          Ydbg_nanf ("F: ", YNBI_sbcgyre_qsr_qns, "sbcgyre_qsr_qns"); 
    38763875         Ysbcgyre_emp1(Yi, Yj)->forward(); 
    38773876         if (1==0 || Yi-1<0) 
    … …  
    38933892         Yting[4]=YS1_sbcgyre_emp1( Yi, Yj, YTemps); 
    38943893         Ysbcgyre_zsumemp_zsurf(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    3895          Ydbg_nanf ("F: ", YNBI_sbcgyre_zsumemp_zsurf, "sbcgyre_zsumemp_zsurf"); 
    38963894 
    38973895} 
    … …  
    39033901         Yting[2]=YS2_sbcgyre_zsumemp_zsurf( 31, 21, YTemps); 
    39043902         Ysbcgyre_emp(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2] ); 
    3905          Ydbg_nanf ("F: ", YNBI_sbcgyre_emp, "sbcgyre_emp"); 
    39063903         Yting[0]=YS1_sbcgyre_emp( Yi, Yj, YTemps); 
    39073904         Ysbcgyre_emps(Yi, Yj)->forward( Yting[0] ); 
    3908          Ydbg_nanf ("F: ", YNBI_sbcgyre_emps, "sbcgyre_emps"); 
    39093905         Ysbcgyre_utau_vtau(Yi, Yj)->forward(); 
    39103906         if (1==0 || Yi-1<0) 
    … …  
    39193915         Yting[3]=YS2_sbcgyre_utau_vtau( Yi, Yj, YTemps); 
    39203916         Ysbcgyre_wndm(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    3921          Ydbg_nanf ("F: ", YNBI_sbcgyre_wndm, "sbcgyre_wndm"); 
    39223917 
    39233918} 
    … …  
    39393934         Yting[3]=YS1_sb( Yi, Yj, Yk-1, YTemps-1); 
    39403935         Ybn2_rn2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    3941          Ydbg_nanf ("F: ", YNBI_bn2_rn2, "bn2_rn2"); 
    39423936 
    39433937} 
    … …  
    40114005         Yting[18]=YS1_va( Yi+1, Yj, 1, YTemps-1); 
    40124006         Yzdf_avt_avmu_avmv(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17], Yting[18] ); 
    4013          Ydbg_nanf ("F: ", YNBI_zdf_avt_avmu_avmv, "zdf_avt_avmu_avmv"); 
    40144007 
    40154008} 
    … …  
    40524045         Yting[30]=YS1_zdf_avt_avmu_avmv( Yi, Yj, 30, YTemps); 
    40534046         Yzdfmxl_hmld(Yi, Yj)->forward(); 
    4054          Ydbg_nanf ("F: ", YNBI_zdfmxl_hmld, "zdfmxl_hmld"); 
    40554047         Yting[0]=YS1_rhop( Yi, Yj, 0, YTemps-1); 
    40564048         Yting[1]=YS1_rhop( Yi, Yj, 1, YTemps-1); 
    … …  
    40854077         Yting[30]=YS1_rhop( Yi, Yj, 30, YTemps-1); 
    40864078         Yzdfmxl_nmln_hmlp_hmlpt(Yi, Yj)->forward(); 
    4087          Ydbg_nanf ("F: ", YNBI_zdfmxl_nmln_hmlp_hmlpt, "zdfmxl_nmln_hmlp_hmlpt"); 
    40884079 
    40894080} 
    … …  
    40984089         Yting[2]=YS1_sa( Yi, Yj, 0, YTemps-1); 
    40994090         Yta_sa_trasbc(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2] ); 
    4100          Ydbg_nanf ("F: ", YNBI_ta_sa_trasbc, "ta_sa_trasbc"); 
    41014091 
    41024092} 
    … …  
    41114101         Yting[1]=YS1_ta_sa_trasbc( Yi, Yj, YTemps); 
    41124102         Yta_traqsr(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4113          Ydbg_nanf ("F: ", YNBI_ta_traqsr, "ta_traqsr"); 
    41144103 
    41154104} 
    … …  
    41344123         Yting[4]=YS1_ta( Yi+1, Yj, Yk, YTemps-1); 
    41354124         Yzcent_hor_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4136          Ydbg_nanf ("F: ", YNBI_zcent_hor_traadv_cen2, "zcent_hor_traadv_cen2"); 
    41374125         Yting[0]=YS1_ua( Yi, Yj, Yk, YTemps-1); 
    41384126         Yting[1]=YS1_va( Yi, Yj, Yk, YTemps-1); 
    … …  
    41474135         Yting[4]=YS1_sa( Yi+1, Yj, Yk, YTemps-1); 
    41484136         Yzcens_hor_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4149          Ydbg_nanf ("F: ", YNBI_zcens_hor_traadv_cen2, "zcens_hor_traadv_cen2"); 
    41504137 
    41514138} 
    … …  
    41654152         Yting[4]=YS2_zcent_hor_traadv_cen2( Yi, Yj, Yk, YTemps); 
    41664153         Yta_hor_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4167          Ydbg_nanf ("F: ", YNBI_ta_hor_traadv_cen2, "ta_hor_traadv_cen2"); 
    41684154         Yting[0]=YS2_ta_sa_trasbc( Yi, Yj, YTemps); 
    41694155         if (1==0 || Yi-1<0) 
    … …  
    41784164         Yting[4]=YS2_zcens_hor_traadv_cen2( Yi, Yj, Yk, YTemps); 
    41794165         Ysa_hor_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4180          Ydbg_nanf ("F: ", YNBI_sa_hor_traadv_cen2, "sa_hor_traadv_cen2"); 
    41814166 
    41824167} 
    … …  
    41974182         Yting[2]=YS1_ta( Yi, Yj, Yk-1, YTemps-1); 
    41984183         Yzcent_ver_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4199          Ydbg_nanf ("F: ", YNBI_zcent_ver_traadv_cen2, "zcent_ver_traadv_cen2"); 
    42004184         Yting[0]=YS1_wa( Yi, Yj, Yk, YTemps-1); 
    42014185         Yting[1]=YS1_sa( Yi, Yj, Yk, YTemps-1); 
    … …  
    42054189         Yting[2]=YS1_sa( Yi, Yj, Yk-1, YTemps-1); 
    42064190         Yzcens_ver_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4207          Ydbg_nanf ("F: ", YNBI_zcens_ver_traadv_cen2, "zcens_ver_traadv_cen2"); 
    42084191 
    42094192} 
    … …  
    42224205         Yting[2]=YS1_zcent_ver_traadv_cen2( Yi, Yj, Yk+1, YTemps); 
    42234206         Yta_ver_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4224          Ydbg_nanf ("F: ", YNBI_ta_ver_traadv_cen2, "ta_ver_traadv_cen2"); 
    42254207         Yting[0]=YS1_sa_hor_traadv_cen2( Yi, Yj, Yk, YTemps); 
    42264208         Yting[1]=YS1_zcens_ver_traadv_cen2( Yi, Yj, Yk, YTemps); 
    … …  
    42304212         Yting[2]=YS1_zcens_ver_traadv_cen2( Yi, Yj, Yk+1, YTemps); 
    42314213         Ysa_ver_traadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4232          Ydbg_nanf ("F: ", YNBI_sa_ver_traadv_cen2, "sa_ver_traadv_cen2"); 
    42334214 
    42344215} 
    … …  
    42504231         Yting[2]=YS1_tb( Yi, Yj+1, Yk, YTemps-1); 
    42514232         Yztuv_traldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4252          Ydbg_nanf ("F: ", YNBI_ztuv_traldf_lap, "ztuv_traldf_lap"); 
    42534233         Yting[0]=YS1_sb( Yi, Yj, Yk, YTemps-1); 
    42544234         if (1==0 || Yi+1>YA1_S3dt-1) 
    … …  
    42614241         Yting[2]=YS1_sb( Yi, Yj+1, Yk, YTemps-1); 
    42624242         Yzsuv_traldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4263          Ydbg_nanf ("F: ", YNBI_zsuv_traldf_lap, "zsuv_traldf_lap"); 
    42644243         Yting[0]=YS1_ta_ver_traadv_cen2( Yi, Yj, Yk, YTemps); 
    42654244         Yting[1]=YS1_ztuv_traldf_lap( Yi, Yj, Yk, YTemps); 
    … …  
    42744253         Yting[4]=YS2_ztuv_traldf_lap( Yi, Yj-1, Yk, YTemps); 
    42754254         Yta_traldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4276          Ydbg_nanf ("F: ", YNBI_ta_traldf_lap, "ta_traldf_lap"); 
    42774255         Yting[0]=YS1_sa_ver_traadv_cen2( Yi, Yj, Yk, YTemps); 
    42784256         Yting[1]=YS1_zsuv_traldf_lap( Yi, Yj, Yk, YTemps); 
    … …  
    42874265         Yting[4]=YS2_zsuv_traldf_lap( Yi, Yj-1, Yk, YTemps); 
    42884266         Ysa_traldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4289          Ydbg_nanf ("F: ", YNBI_sa_traldf_lap, "sa_traldf_lap"); 
    42904267 
    42914268} 
    … …  
    43034280         Yting[1]=YS1_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps); 
    43044281         Yzwisd_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4305          Ydbg_nanf ("F: ", YNBI_zwisd_trazdf_imp, "zwisd_trazdf_imp"); 
    43064282         Yting[0]=YS1_ta_traldf_lap( Yi, Yj, Yk, YTemps); 
    43074283         Yting[1]=YS1_tb( Yi, Yj, Yk, YTemps-1); 
    43084284         Yzrhs_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4309          Ydbg_nanf ("F: ", YNBI_zrhs_trazdf_imp, "zrhs_trazdf_imp"); 
    43104285         Yting[0]=YS1_sa_traldf_lap( Yi, Yj, Yk, YTemps); 
    43114286         Yting[1]=YS1_sb( Yi, Yj, Yk, YTemps-1); 
    43124287         Yzrhs_sa_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4313          Ydbg_nanf ("F: ", YNBI_zrhs_sa_trazdf_imp, "zrhs_sa_trazdf_imp"); 
    43144288 
    43154289} 
    … …  
    43304304         Yting[3]=YS1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps); 
    43314305         Yzwt_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4332          Ydbg_nanf ("F: ", YNBI_zwt_trazdf_imp, "zwt_trazdf_imp"); 
    43334306         Yting[0]=YS1_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps); 
    43344307         if (1==0 || Yk-1<0) 
    … …  
    43424315         Yting[3]=YS1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps); 
    43434316         Yta_de_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4344          Ydbg_nanf ("F: ", YNBI_ta_de_trazdf_imp, "ta_de_trazdf_imp"); 
    43454317         Yting[0]=YS1_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps); 
    43464318         if (1==0 || Yk-1<0) 
    … …  
    43544326         Yting[3]=YS1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps); 
    43554327         Ysa_de_trazdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4356          Ydbg_nanf ("F: ", YNBI_sa_de_trazdf_imp, "sa_de_trazdf_imp"); 
    43574328 
    43584329} 
    … …  
    43704341         Yting[4]=YS1_ta_c( Yi, Yj, Yk); 
    43714342         Yta(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4372          Ydbg_nanf ("F: ", YNBI_ta, "ta"); 
    43734343         Yting[0]=YS2_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps); 
    43744344         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    43804350         Yting[4]=YS1_sa_c( Yi, Yj, Yk); 
    43814351         Ysa(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4382          Ydbg_nanf ("F: ", YNBI_sa, "sa"); 
    43834352 
    43844353} 
    … …  
    43944363         Yting[3]=YS1_ta_c( Yi, Yj, Yk); 
    43954364         Ytb(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4396          Ydbg_nanf ("F: ", YNBI_tb, "tb"); 
    43974365         Yting[0]=YS1_sb( Yi, Yj, Yk, YTemps-1); 
    43984366         Yting[1]=YS1_sa( Yi, Yj, Yk, YTemps-1); 
    … …  
    44004368         Yting[3]=YS1_sa_c( Yi, Yj, Yk); 
    44014369         Ysb(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4402          Ydbg_nanf ("F: ", YNBI_sb, "sb"); 
    44034370 
    44044371} 
    … …  
    44134380         Yting[1]=YS1_sb( Yi, Yj, Yk, YTemps); 
    44144381         Yrhd(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4415          Ydbg_nanf ("F: ", YNBI_rhd, "rhd"); 
    44164382         Yting[0]=YS1_rhd( Yi, Yj, Yk, YTemps); 
    44174383         Yrhop(Yi, Yj, Yk)->forward( Yting[0] ); 
    4418          Ydbg_nanf ("F: ", YNBI_rhop, "rhop"); 
    44194384 
    44204385} 
    … …  
    44464411         Yting[5]=YS1_va( Yi, Yj+1, Yk, YTemps-1); 
    44474412         Yhmf_TF_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5] ); 
    4448          Ydbg_nanf ("F: ", YNBI_hmf_TF_dynadv_cen2, "hmf_TF_dynadv_cen2"); 
    44494413 
    44504414} 
    … …  
    44634427         Yting[3]=YS2_hmf_TF_dynadv_cen2( Yi, Yj-1, Yk, YTemps); 
    44644428         Yzua_dhmf_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4465          Ydbg_nanf ("F: ", YNBI_zua_dhmf_dynadv_cen2, "zua_dhmf_dynadv_cen2"); 
    44664429         Yting[0]=YS3_hmf_TF_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    44674430         if (1==0 || Yi-1<0) 
    … …  
    44754438         Yting[3]=YS4_hmf_TF_dynadv_cen2( Yi, Yj-1, Yk, YTemps); 
    44764439         Yzva_dhmf_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4477          Ydbg_nanf ("F: ", YNBI_zva_dhmf_dynadv_cen2, "zva_dhmf_dynadv_cen2"); 
    44784440 
    44794441} 
    … …  
    45074469         Yting[6]=YS1_va( Yi, Yj, Yk-1, YTemps-1); 
    45084470         Yzfuv_UW_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6] ); 
    4509          Ydbg_nanf ("F: ", YNBI_zfuv_UW_dynadv_cen2, "zfuv_UW_dynadv_cen2"); 
    45104471 
    45114472} 
    … …  
    45244485         Yting[2]=YS1_zua_dhmf_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    45254486         Yua_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4526          Ydbg_nanf ("F: ", YNBI_ua_dynadv_cen2, "ua_dynadv_cen2"); 
    45274487         Yting[0]=YS2_zfuv_UW_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    45284488         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    45324492         Yting[2]=YS1_zva_dhmf_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    45334493         Yva_dynadv_cen2(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4534          Ydbg_nanf ("F: ", YNBI_va_dynadv_cen2, "va_dynadv_cen2"); 
    45354494 
    45364495} 
    … …  
    45544513         Yting[3]=YS1_va( Yi+1, Yj, Yk, YTemps-1); 
    45554514         Yzwxyz_vor_ene(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4556          Ydbg_nanf ("F: ", YNBI_zwxyz_vor_ene, "zwxyz_vor_ene"); 
    45574515 
    45584516} 
    … …  
    45804538         Yting[6]=YS2_zwxyz_vor_ene( Yi+1, Yj, Yk, YTemps); 
    45814539         Yua_vor_ene(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6] ); 
    4582          Ydbg_nanf ("F: ", YNBI_ua_vor_ene, "ua_vor_ene"); 
    45834540         Yting[0]=YS1_va_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    45844541         if (1==0 || Yi-1<0) 
    … …  
    46014558         Yting[6]=YS1_zwxyz_vor_ene( Yi, Yj+1, Yk, YTemps); 
    46024559         Yva_vor_ene(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6] ); 
    4603          Ydbg_nanf ("F: ", YNBI_va_vor_ene, "va_vor_ene"); 
    46044560         Yting[0]=YS1_ua_vor_ene( Yi, Yj, Yk, YTemps); 
    46054561         if (1==0 || Yj-1<0) 
    … …  
    46144570         Yting[4]=YS1_hdivn( Yi+1, Yj, Yk, YTemps-2); 
    46154571         Yua_dyn_ldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4616          Ydbg_nanf ("F: ", YNBI_ua_dyn_ldf_lap, "ua_dyn_ldf_lap"); 
    46174572         Yting[0]=YS1_va_vor_ene( Yi, Yj, Yk, YTemps); 
    46184573         if (1==0 || Yi-1<0) 
    … …  
    46274582         Yting[4]=YS1_hdivn( Yi, Yj+1, Yk, YTemps-2); 
    46284583         Yva_dyn_ldf_lap(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    4629          Ydbg_nanf ("F: ", YNBI_va_dyn_ldf_lap, "va_dyn_ldf_lap"); 
    46304584 
    46314585} 
    … …  
    46674621         Yting[7]=YS2_zhpij_dyn_hpg_zco( Yi, Yj, Yk-1, YTemps); 
    46684622         Yzhpij_dyn_hpg_zco(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7] ); 
    4669          Ydbg_nanf ("F: ", YNBI_zhpij_dyn_hpg_zco, "zhpij_dyn_hpg_zco"); 
    46704623         Yting[0]=YS1_ua_dyn_ldf_lap( Yi, Yj, Yk, YTemps); 
    46714624         Yting[1]=YS1_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps); 
    46724625         Yua_dyn_hpg_zco(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4673          Ydbg_nanf ("F: ", YNBI_ua_dyn_hpg_zco, "ua_dyn_hpg_zco"); 
    46744626         Yting[0]=YS1_va_dyn_ldf_lap( Yi, Yj, Yk, YTemps); 
    46754627         Yting[1]=YS2_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps); 
    46764628         Yva_dyn_hpg_zco(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4677          Ydbg_nanf ("F: ", YNBI_va_dyn_hpg_zco, "va_dyn_hpg_zco"); 
    46784629 
    46794630} 
    … …  
    46904641         Yting[1]=YS2_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps); 
    46914642         Yzwisd_ua_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4692          Ydbg_nanf ("F: ", YNBI_zwisd_ua_dynzdf_imp, "zwisd_ua_dynzdf_imp"); 
    46934643         Yting[0]=YS3_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps); 
    46944644         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    46974647         Yting[1]=YS3_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps); 
    46984648         Yzwisd_va_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4699          Ydbg_nanf ("F: ", YNBI_zwisd_va_dynzdf_imp, "zwisd_va_dynzdf_imp"); 
    47004649         Yting[0]=YS1_ua_dyn_hpg_zco( Yi, Yj, Yk, YTemps); 
    47014650         Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1); 
    47024651         Yting[2]=YS1_sbcgyre_utau_vtau( Yi, Yj, YTemps); 
    47034652         Yzrhs_ua_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4704          Ydbg_nanf ("F: ", YNBI_zrhs_ua_dynzdf_imp, "zrhs_ua_dynzdf_imp"); 
    47054653         Yting[0]=YS1_va_dyn_hpg_zco( Yi, Yj, Yk, YTemps); 
    47064654         Yting[1]=YS1_vb( Yi, Yj, Yk, YTemps-1); 
    47074655         Yting[2]=YS2_sbcgyre_utau_vtau( Yi, Yj, YTemps); 
    47084656         Yzrhs_va_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4709          Ydbg_nanf ("F: ", YNBI_zrhs_va_dynzdf_imp, "zrhs_va_dynzdf_imp"); 
    47104657 
    47114658} 
    … …  
    47264673         Yting[3]=YS1_zwt_ua_dynzdf_imp( Yi, Yj, Yk-1, YTemps); 
    47274674         Yzwt_ua_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4728          Ydbg_nanf ("F: ", YNBI_zwt_ua_dynzdf_imp, "zwt_ua_dynzdf_imp"); 
    47294675         Yting[0]=YS1_zwisd_ua_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    47304676         if (1==0 || Yk-1<0) 
    … …  
    47384684         Yting[3]=YS1_zwt_ua_dynzdf_imp( Yi, Yj, Yk-1, YTemps); 
    47394685         Yua_de_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4740          Ydbg_nanf ("F: ", YNBI_ua_de_dynzdf_imp, "ua_de_dynzdf_imp"); 
    47414686         Yting[0]=YS1_zwisd_va_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    47424687         if (1==0 || Yk-1<0) 
    … …  
    47504695         Yting[3]=YS1_zwt_va_dynzdf_imp( Yi, Yj, Yk-1, YTemps); 
    47514696         Yzwt_va_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4752          Ydbg_nanf ("F: ", YNBI_zwt_va_dynzdf_imp, "zwt_va_dynzdf_imp"); 
    47534697         Yting[0]=YS1_zwisd_va_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    47544698         if (1==0 || Yk-1<0) 
    … …  
    47624706         Yting[3]=YS1_zwt_va_dynzdf_imp( Yi, Yj, Yk-1, YTemps); 
    47634707         Yva_de_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4764          Ydbg_nanf ("F: ", YNBI_va_de_dynzdf_imp, "va_de_dynzdf_imp"); 
    47654708 
    47664709} 
    … …  
    47774720         Yting[3]=YS1_zwt_ua_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    47784721         Yua_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4779          Ydbg_nanf ("F: ", YNBI_ua_dynzdf_imp, "ua_dynzdf_imp"); 
    47804722         Yting[0]=YS2_zwisd_va_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    47814723         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    47864728         Yting[3]=YS1_zwt_va_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    47874729         Yva_dynzdf_imp(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4788          Ydbg_nanf ("F: ", YNBI_va_dynzdf_imp, "va_dynzdf_imp"); 
    47894730 
    47904731} 
    … …  
    47984739         Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1); 
    47994740         Yua_norm_dynzdf(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4800          Ydbg_nanf ("F: ", YNBI_ua_norm_dynzdf, "ua_norm_dynzdf"); 
    48014741         Yting[0]=YS1_va_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    48024742         Yting[1]=YS1_vb( Yi, Yj, Yk, YTemps-1); 
    48034743         Yva_norm_dynzdf(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    4804          Ydbg_nanf ("F: ", YNBI_va_norm_dynzdf, "va_norm_dynzdf"); 
    48054744 
    48064745} 
    … …  
    48214760         Yting[2]=YS1_sshn( Yi+1, Yj, YTemps-1); 
    48224761         Yspguv_now_dynspg(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2] ); 
    4823          Ydbg_nanf ("F: ", YNBI_spguv_now_dynspg, "spguv_now_dynspg"); 
    48244762 
    48254763} 
    … …  
    48354773         Yting[2]=YS1_spguv_now_dynspg( Yi, Yj, YTemps); 
    48364774         Yua_1_dynspg_flt(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4837          Ydbg_nanf ("F: ", YNBI_ua_1_dynspg_flt, "ua_1_dynspg_flt"); 
    48384775         Yting[0]=YS1_va_norm_dynzdf( Yi, Yj, Yk, YTemps); 
    48394776         Yting[1]=YS1_vb( Yi, Yj, Yk, YTemps-1); 
    48404777         Yting[2]=YS2_spguv_now_dynspg( Yi, Yj, YTemps); 
    48414778         Yva_1_dynspg_flt(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    4842          Ydbg_nanf ("F: ", YNBI_va_1_dynspg_flt, "va_1_dynspg_flt"); 
    48434779 
    48444780} 
    … …  
    48814817         Yting[30]=YS1_ua_1_dynspg_flt( Yi, Yj, 30, YTemps); 
    48824818         Ysum_spgu_dynspg_flt(Yi, Yj)->forward(); 
    4883          Ydbg_nanf ("F: ", YNBI_sum_spgu_dynspg_flt, "sum_spgu_dynspg_flt"); 
    48844819         Yting[0]=YS1_va_1_dynspg_flt( Yi, Yj, 0, YTemps); 
    48854820         Yting[1]=YS1_va_1_dynspg_flt( Yi, Yj, 1, YTemps); 
    … …  
    49144849         Yting[30]=YS1_va_1_dynspg_flt( Yi, Yj, 30, YTemps); 
    49154850         Ysum_spgv_dynspg_flt(Yi, Yj)->forward(); 
    4916          Ydbg_nanf ("F: ", YNBI_sum_spgv_dynspg_flt, "sum_spgv_dynspg_flt"); 
    49174851 
    49184852} 
    … …  
    49344868         Yting[3]=YS1_sum_spgv_dynspg_flt( Yi, Yj-1, YTemps); 
    49354869         Ygcb_dynspg_flt(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    4936          Ydbg_nanf ("F: ", YNBI_gcb_dynspg_flt, "gcb_dynspg_flt"); 
    49374870         Yting[0]=YS1_gcx2( Yi, Yj, YTemps-1); 
    49384871         Yting[1]=YS1_gcx2( Yi, Yj, YTemps-2); 
    49394872         Ygcx_dynspg_flt(Yi, Yj)->forward( Yting[0], Yting[1] ); 
    4940          Ydbg_nanf ("F: ", YNBI_gcx_dynspg_flt, "gcx_dynspg_flt"); 
    49414873 
    49424874} 
    … …  
    50164948         Yting[18]=YS1_solsor_dynspg_flt( Yi-1, Yj-1, YTemps); 
    50174949         Ysolsor_dynspg_flt(Yi, Yj)->forward(); 
    5018          Ydbg_nanf ("F: ", YNBI_solsor_dynspg_flt, "solsor_dynspg_flt"); 
    50194950 
    50204951} 
    … …  
    50264957         Yting[0]=YS1_solsor_dynspg_flt( Yi, Yj, YTemps); 
    50274958         Ygcx2(Yi, Yj)->forward(); 
    5028          Ydbg_nanf ("F: ", YNBI_gcx2, "gcx2"); 
    50294959 
    50304960} 
    … …  
    50454975         Yting[2]=YS1_gcx2( Yi, Yj+1, YTemps); 
    50464976         Yspguv_dynspg_flt(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2] ); 
    5047          Ydbg_nanf ("F: ", YNBI_spguv_dynspg_flt, "spguv_dynspg_flt"); 
    50484977 
    50494978} 
    … …  
    50594988         Yting[2]=YS1_ua_c( Yi, Yj, Yk); 
    50604989         Yua(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    5061          Ydbg_nanf ("F: ", YNBI_ua, "ua"); 
    50624990         Yting[0]=YS1_va_1_dynspg_flt( Yi, Yj, Yk, YTemps); 
    50634991         Yting[1]=YS2_spguv_dynspg_flt( Yi, Yj, YTemps); 
    50644992         Yting[2]=YS1_va_c( Yi, Yj, Yk); 
    50654993         Yva(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2] ); 
    5066          Ydbg_nanf ("F: ", YNBI_va, "va"); 
    50674994 
    50684995} 
    … …  
    50785005         Yting[3]=YS1_sshn_c( Yi, Yj); 
    50795006         Ysshn(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    5080          Ydbg_nanf ("F: ", YNBI_sshn, "sshn"); 
    50815007         Yting[0]=YS1_sshb( Yi, Yj, YTemps-1); 
    50825008         Yting[1]=YS1_sshn( Yi, Yj, YTemps-1); 
    … …  
    50845010         Yting[3]=YS1_sshn_c( Yi, Yj); 
    50855011         Ysshb(Yi, Yj)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    5086          Ydbg_nanf ("F: ", YNBI_sshb, "sshb"); 
    50875012 
    50885013} 
    … …  
    50995024         Yting[3]=YS1_ua_c( Yi, Yj, Yk); 
    51005025         Yub(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    5101          Ydbg_nanf ("F: ", YNBI_ub, "ub"); 
    51025026         Yting[0]=YS1_vb( Yi, Yj, Yk, YTemps-1); 
    51035027         Yting[1]=YS1_va( Yi, Yj, Yk, YTemps-1); 
    … …  
    51055029         Yting[3]=YS1_va_c( Yi, Yj, Yk); 
    51065030         Yvb(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    5107          Ydbg_nanf ("F: ", YNBI_vb, "vb"); 
    51085031 
    51095032} 
    … …  
    51265049         Yting[3]=YS1_va( Yi, Yj-1, Yk, YTemps); 
    51275050         Yhdivn(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    5128          Ydbg_nanf ("F: ", YNBI_hdivn, "hdivn"); 
    51295051         if (1==0 || Yi+1>YA1_S3dt-1) 
    51305052         Yting[0]=0;  
    … …  
    51385060         Yting[3]=YS1_ua( Yi, Yj, Yk, YTemps); 
    51395061         Yrotn(Yi, Yj, Yk)->forward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    5140          Ydbg_nanf ("F: ", YNBI_rotn, "rotn"); 
    51415062 
    51425063} 
    … …  
    51525073         Yting[1]=YS1_hdivn( Yi, Yj, Yk, YTemps); 
    51535074         Ywa(Yi, Yj, Yk)->forward( Yting[0], Yting[1] ); 
    5154          Ydbg_nanf ("F: ", YNBI_wa, "wa"); 
    51555075 
    51565076} 
    … …  
    51795099         YG1_wa( Yi, Yj, Yk+1, YTemps) += Ytbeta[0]; 
    51805100         YG1_hdivn( Yi, Yj, Yk, YTemps) += Ytbeta[1]; 
    5181          Ydbg_nanf ("B: ", YNBI_wa, "wa"); 
    51825101 
    51835102} 
    … …  
    52095128         YG1_ua( Yi, Yj+1, Yk, YTemps) += Ytbeta[2]; 
    52105129         YG1_ua( Yi, Yj, Yk, YTemps) += Ytbeta[3]; 
    5211          Ydbg_nanf ("B: ", YNBI_rotn, "rotn"); 
    52125130         Yting[0]=YS1_ua( Yi, Yj, Yk, YTemps); 
    52135131         if (1==0 || Yi-1<0) 
    … …  
    52315149         else  
    52325150         YG1_va( Yi, Yj-1, Yk, YTemps) += Ytbeta[3]; 
    5233          Ydbg_nanf ("B: ", YNBI_hdivn, "hdivn"); 
    52345151 
    52355152} 
    … …  
    52535170         YG1_va( Yi, Yj, Yk, YTemps) += Ytbeta[2]; 
    52545171         YG1_va_c( Yi, Yj, Yk) += Ytbeta[3]; 
    5255          Ydbg_nanf ("B: ", YNBI_vb, "vb"); 
    52565172         Yting[0]=YS1_ub( Yi, Yj, Yk, YTemps-1); 
    52575173         Yting[1]=YS1_ua( Yi, Yj, Yk, YTemps-1); 
    … …  
    52665182         YG1_ua( Yi, Yj, Yk, YTemps) += Ytbeta[2]; 
    52675183         YG1_ua_c( Yi, Yj, Yk) += Ytbeta[3]; 
    5268          Ydbg_nanf ("B: ", YNBI_ub, "ub"); 
    52695184 
    52705185} 
    … …  
    52875202         YG1_sshn( Yi, Yj, YTemps) += Ytbeta[2]; 
    52885203         YG1_sshn_c( Yi, Yj) += Ytbeta[3]; 
    5289          Ydbg_nanf ("B: ", YNBI_sshb, "sshb"); 
    52905204         Yting[0]=YS1_sshb( Yi, Yj, YTemps-1); 
    52915205         Yting[1]=YS1_wa( Yi, Yj, 0, YTemps-1); 
    … …  
    53005214         YG1_sbcgyre_emp( Yi, Yj, YTemps) += Ytbeta[2]; 
    53015215         YG1_sshn_c( Yi, Yj) += Ytbeta[3]; 
    5302          Ydbg_nanf ("B: ", YNBI_sshn, "sshn"); 
    53035216 
    53045217} 
    … …  
    53205233         YG2_spguv_dynspg_flt( Yi, Yj, YTemps) += Ytbeta[1]; 
    53215234         YG1_va_c( Yi, Yj, Yk) += Ytbeta[2]; 
    5322          Ydbg_nanf ("B: ", YNBI_va, "va"); 
    53235235         Yting[0]=YS1_ua_1_dynspg_flt( Yi, Yj, Yk, YTemps); 
    53245236         Yting[1]=YS1_spguv_dynspg_flt( Yi, Yj, YTemps); 
    … …  
    53315243         YG1_spguv_dynspg_flt( Yi, Yj, YTemps) += Ytbeta[1]; 
    53325244         YG1_ua_c( Yi, Yj, Yk) += Ytbeta[2]; 
    5333          Ydbg_nanf ("B: ", YNBI_ua, "ua"); 
    53345245 
    53355246} 
    … …  
    53625273         #pragma omp atomic 
    53635274         YG1_gcx2( Yi, Yj+1, YTemps) += Ytbeta[2]; 
    5364          Ydbg_nanf ("B: ", YNBI_spguv_dynspg_flt, "spguv_dynspg_flt"); 
    53655275 
    53665276} 
    … …  
    53775287         #pragma omp atomic 
    53785288         YG1_solsor_dynspg_flt( Yi, Yj, YTemps) += Ytbeta[0]; 
    5379          Ydbg_nanf ("B: ", YNBI_gcx2, "gcx2"); 
    53805289 
    53815290} 
    … …  
    55105419         else  
    55115420         YG1_solsor_dynspg_flt( Yi-1, Yj-1, YTemps) += Ytbeta[18]; 
    5512          Ydbg_nanf ("B: ", YNBI_solsor_dynspg_flt, "solsor_dynspg_flt"); 
    55135421 
    55145422} 
    … …  
    55255433         YG1_gcx2( Yi, Yj, YTemps-1) += Ytbeta[0]; 
    55265434         YG1_gcx2( Yi, Yj, YTemps-2) += Ytbeta[1]; 
    5527          Ydbg_nanf ("B: ", YNBI_gcx_dynspg_flt, "gcx_dynspg_flt"); 
    55285435         Yting[0]=YS1_sum_spgu_dynspg_flt( Yi, Yj, YTemps); 
    55295436         if (1==0 || Yi-1<0) 
    … …  
    55495456         #pragma omp atomic 
    55505457         YG1_sum_spgv_dynspg_flt( Yi, Yj-1, YTemps) += Ytbeta[3]; 
    5551          Ydbg_nanf ("B: ", YNBI_gcb_dynspg_flt, "gcb_dynspg_flt"); 
    55525458 
    55535459} 
    … …  
    56235529         YG1_va_1_dynspg_flt( Yi, Yj, 29, YTemps) += Ytbeta[29]; 
    56245530         YG1_va_1_dynspg_flt( Yi, Yj, 30, YTemps) += Ytbeta[30]; 
    5625          Ydbg_nanf ("B: ", YNBI_sum_spgv_dynspg_flt, "sum_spgv_dynspg_flt"); 
    56265531         Yting[0]=YS1_ua_1_dynspg_flt( Yi, Yj, 0, YTemps); 
    56275532         Yting[1]=YS1_ua_1_dynspg_flt( Yi, Yj, 1, YTemps); 
    … …  
    56895594         YG1_ua_1_dynspg_flt( Yi, Yj, 29, YTemps) += Ytbeta[29]; 
    56905595         YG1_ua_1_dynspg_flt( Yi, Yj, 30, YTemps) += Ytbeta[30]; 
    5691          Ydbg_nanf ("B: ", YNBI_sum_spgu_dynspg_flt, "sum_spgu_dynspg_flt"); 
    5692  
    5693 } 
    5694 } 
    5695 //#pragma omp parallel for num_threads(16) 
     5596 
     5597} 
     5598} 
     5599#pragma omp parallel for num_threads(16) 
    56965600for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    56975601  Yk=YY; 
    … …  
    57085612         YG1_vb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    57095613         YG2_spguv_now_dynspg( Yi, Yj, YTemps) += Ytbeta[2]; 
    5710          Ydbg_nanf ("B: ", YNBI_va_1_dynspg_flt, "va_1_dynspg_flt"); 
    57115614         Yting[0]=YS1_ua_norm_dynzdf( Yi, Yj, Yk, YTemps); 
    57125615         Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1); 
    … …  
    57185621         YG1_ub( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    57195622         YG1_spguv_now_dynspg( Yi, Yj, YTemps) += Ytbeta[2]; 
    5720          Ydbg_nanf ("B: ", YNBI_ua_1_dynspg_flt, "ua_1_dynspg_flt"); 
    5721  
    5722 } 
    5723 } 
    5724 //#pragma omp parallel for num_threads(16) 
     5623 
     5624} 
     5625} 
     5626#pragma omp parallel for num_threads(16) 
    57255627for(YY=YA2_S3dt-1; YY>=0; --YY){ 
    57265628  Yj=YY; 
    … …  
    57395641         Yspguv_now_dynspg(Yi, Yj)->backward( Yting[0], Yting[1], Yting[2] ); 
    57405642         Yvsmatt (YNBS_spguv_now_dynspg, YNBI_spguv_now_dynspg, YMAX_JAC_NBI, &YG1Y_spguv_now_dynspg, Yjac[0], Ytbeta); 
    5741          //      #pragma omp atomic 
     5643         #pragma omp atomic 
    57425644         YG1_sshn( Yi, Yj, YTemps-1) += Ytbeta[0]; 
    57435645         if (1==0 || Yj+1>YA2_S3dt-1){}  
    57445646         else  
    5745            //    #pragma omp atomic 
     5647         #pragma omp atomic 
    57465648         YG1_sshn( Yi, Yj+1, YTemps-1) += Ytbeta[1]; 
    57475649         if (1==0 || Yi+1>YA1_S3dt-1){}  
    57485650         else  
    5749            //    #pragma omp atomic 
     5651         #pragma omp atomic 
    57505652         YG1_sshn( Yi+1, Yj, YTemps-1) += Ytbeta[2]; 
    5751          Ydbg_nanf ("B: ", YNBI_spguv_now_dynspg, "spguv_now_dynspg"); 
    5752  
    5753 } 
    5754 } 
    5755 //#pragma omp parallel for num_threads(16) 
     5653 
     5654} 
     5655} 
     5656#pragma omp parallel for num_threads(16) 
    57565657for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    57575658  Yk=YY; 
    … …  
    57645665         Yva_norm_dynzdf(Yi, Yj, Yk)->backward( Yting[0], Yting[1] ); 
    57655666         Yvsmatt (YNBS_va_norm_dynzdf, YNBI_va_norm_dynzdf, YMAX_JAC_NBI, &YG1Y_va_norm_dynzdf, Yjac[0], Ytbeta); 
    5766          //      #pragma omp atomic 
     5667         #pragma omp atomic 
    57675668         YG1_va_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    57685669         YG1_vb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    5769          Ydbg_nanf ("B: ", YNBI_va_norm_dynzdf, "va_norm_dynzdf"); 
    57705670         Yting[0]=YS1_ua_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    57715671         Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1); 
    … …  
    57735673         Yua_norm_dynzdf(Yi, Yj, Yk)->backward( Yting[0], Yting[1] ); 
    57745674         Yvsmatt (YNBS_ua_norm_dynzdf, YNBI_ua_norm_dynzdf, YMAX_JAC_NBI, &YG1Y_ua_norm_dynzdf, Yjac[0], Ytbeta); 
    5775          //      #pragma omp atomic 
     5675         #pragma omp atomic 
    57765676         YG1_ua_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    57775677         YG1_ub( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    5778          Ydbg_nanf ("B: ", YNBI_ua_norm_dynzdf, "ua_norm_dynzdf"); 
    57795678 
    57805679} 
    … …  
    58005699         YG1_va_de_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[2]; 
    58015700         YG1_zwt_va_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[3]; 
    5802          Ydbg_nanf ("B: ", YNBI_va_dynzdf_imp, "va_dynzdf_imp"); 
    58035701         Yting[0]=YS2_zwisd_ua_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    58045702         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    58175715         YG1_ua_de_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[2]; 
    58185716         YG1_zwt_ua_dynzdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[3]; 
    5819          Ydbg_nanf ("B: ", YNBI_ua_dynzdf_imp, "ua_dynzdf_imp"); 
    58205717 
    58215718} 
    … …  
    58455742         else  
    58465743         YG1_zwt_va_dynzdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3]; 
    5847          Ydbg_nanf ("B: ", YNBI_va_de_dynzdf_imp, "va_de_dynzdf_imp"); 
    58485744         Yting[0]=YS1_zwisd_va_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    58495745         if (1==0 || Yk-1<0) 
    … …  
    58675763         else  
    58685764         YG1_zwt_va_dynzdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3]; 
    5869          Ydbg_nanf ("B: ", YNBI_zwt_va_dynzdf_imp, "zwt_va_dynzdf_imp"); 
    58705765         Yting[0]=YS1_zwisd_ua_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    58715766         if (1==0 || Yk-1<0) 
    … …  
    58895784         else  
    58905785         YG1_zwt_ua_dynzdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3]; 
    5891          Ydbg_nanf ("B: ", YNBI_ua_de_dynzdf_imp, "ua_de_dynzdf_imp"); 
    58925786         Yting[0]=YS1_zwisd_ua_dynzdf_imp( Yi, Yj, Yk, YTemps); 
    58935787         if (1==0 || Yk-1<0) 
    … …  
    59115805         else  
    59125806         YG1_zwt_ua_dynzdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3]; 
    5913          Ydbg_nanf ("B: ", YNBI_zwt_ua_dynzdf_imp, "zwt_ua_dynzdf_imp"); 
    5914  
    5915 } 
    5916 //#pragma omp parallel for num_threads(16) 
     5807 
     5808} 
     5809#pragma omp parallel for num_threads(16) 
    59175810for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    59185811  Yk=YY; 
    … …  
    59295822         YG1_vb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    59305823         YG2_sbcgyre_utau_vtau( Yi, Yj, YTemps) += Ytbeta[2]; 
    5931          Ydbg_nanf ("B: ", YNBI_zrhs_va_dynzdf_imp, "zrhs_va_dynzdf_imp"); 
    59325824         Yting[0]=YS1_ua_dyn_hpg_zco( Yi, Yj, Yk, YTemps); 
    59335825         Yting[1]=YS1_ub( Yi, Yj, Yk, YTemps-1); 
    … …  
    59395831         YG1_ub( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    59405832         YG1_sbcgyre_utau_vtau( Yi, Yj, YTemps) += Ytbeta[2]; 
    5941          Ydbg_nanf ("B: ", YNBI_zrhs_ua_dynzdf_imp, "zrhs_ua_dynzdf_imp"); 
    59425833         Yting[0]=YS3_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps); 
    59435834         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    59485839         Yzwisd_va_dynzdf_imp(Yi, Yj, Yk)->backward( Yting[0], Yting[1] ); 
    59495840         Yvsmatt (YNBS_zwisd_va_dynzdf_imp, YNBI_zwisd_va_dynzdf_imp, YMAX_JAC_NBI, &YG1Y_zwisd_va_dynzdf_imp, Yjac[0], Ytbeta); 
    5950          //      #pragma omp atomic 
     5841         #pragma omp atomic 
    59515842         YG3_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    59525843         if (1==0 || Yk+1>YA3_S3dt-1){}  
    59535844         else  
    5954            //    #pragma omp atomic 
     5845         #pragma omp atomic 
    59555846         YG3_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps) += Ytbeta[1]; 
    5956          Ydbg_nanf ("B: ", YNBI_zwisd_va_dynzdf_imp, "zwisd_va_dynzdf_imp"); 
    59575847         Yting[0]=YS2_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps); 
    59585848         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    59635853         Yzwisd_ua_dynzdf_imp(Yi, Yj, Yk)->backward( Yting[0], Yting[1] ); 
    59645854         Yvsmatt (YNBS_zwisd_ua_dynzdf_imp, YNBI_zwisd_ua_dynzdf_imp, YMAX_JAC_NBI, &YG1Y_zwisd_ua_dynzdf_imp, Yjac[0], Ytbeta); 
    5965          //      #pragma omp atomic 
     5855         #pragma omp atomic 
    59665856         YG2_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    59675857         if (1==0 || Yk+1>YA3_S3dt-1){}  
    59685858         else  
    5969            //    #pragma omp atomic 
     5859         #pragma omp atomic 
    59705860         YG2_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps) += Ytbeta[1]; 
    5971          Ydbg_nanf ("B: ", YNBI_zwisd_ua_dynzdf_imp, "zwisd_ua_dynzdf_imp"); 
    59725861 
    59735862} 
    … …  
    59845873         YG1_va_dyn_ldf_lap( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    59855874         YG2_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps) += Ytbeta[1]; 
    5986          Ydbg_nanf ("B: ", YNBI_va_dyn_hpg_zco, "va_dyn_hpg_zco"); 
    59875875         Yting[0]=YS1_ua_dyn_ldf_lap( Yi, Yj, Yk, YTemps); 
    59885876         Yting[1]=YS1_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps); 
    … …  
    59925880         YG1_ua_dyn_ldf_lap( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    59935881         YG1_zhpij_dyn_hpg_zco( Yi, Yj, Yk, YTemps) += Ytbeta[1]; 
    5994          Ydbg_nanf ("B: ", YNBI_ua_dyn_hpg_zco, "ua_dyn_hpg_zco"); 
    59955882         if (1==0 || Yi+1>YA1_S3dt-1) 
    59965883         Yting[0]=0;  
    … …  
    60475934         else  
    60485935         YG2_zhpij_dyn_hpg_zco( Yi, Yj, Yk-1, YTemps) += Ytbeta[7]; 
    6049          Ydbg_nanf ("B: ", YNBI_zhpij_dyn_hpg_zco, "zhpij_dyn_hpg_zco"); 
    6050  
    6051 } 
    6052 //#pragma omp parallel for num_threads(16) 
     5936 
     5937} 
     5938#pragma omp parallel for num_threads(16) 
    60535939for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    60545940  Yk=YY; 
    … …  
    60805966         else  
    60815967         YG1_hdivn( Yi, Yj+1, Yk, YTemps-2) += Ytbeta[4]; 
    6082          Ydbg_nanf ("B: ", YNBI_va_dyn_ldf_lap, "va_dyn_ldf_lap"); 
    60835968         Yting[0]=YS1_ua_vor_ene( Yi, Yj, Yk, YTemps); 
    60845969         if (1==0 || Yj-1<0) 
    … …  
    61045989         else  
    61055990         YG1_hdivn( Yi+1, Yj, Yk, YTemps-2) += Ytbeta[4]; 
    6106          Ydbg_nanf ("B: ", YNBI_ua_dyn_ldf_lap, "ua_dyn_ldf_lap"); 
    61075991         Yting[0]=YS1_va_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    61085992         if (1==0 || Yi-1<0) 
    … …  
    61426026         else  
    61436027         YG1_zwxyz_vor_ene( Yi, Yj+1, Yk, YTemps) += Ytbeta[6]; 
    6144          Ydbg_nanf ("B: ", YNBI_va_vor_ene, "va_vor_ene"); 
    61456028         Yting[0]=YS1_ua_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    61466029         if (1==0 || Yj-1<0) 
    … …  
    61806063         else  
    61816064         YG2_zwxyz_vor_ene( Yi+1, Yj, Yk, YTemps) += Ytbeta[6]; 
    6182          Ydbg_nanf ("B: ", YNBI_ua_vor_ene, "ua_vor_ene"); 
    61836065 
    61846066} 
    … …  
    61996081         Yzwxyz_vor_ene(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    62006082         Yvsmatt (YNBS_zwxyz_vor_ene, YNBI_zwxyz_vor_ene, YMAX_JAC_NBI, &YG1Y_zwxyz_vor_ene, Yjac[0], Ytbeta); 
    6201          //      #pragma omp atomic 
     6083         #pragma omp atomic 
    62026084         YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    62036085         if (1==0 || Yj+1>YA2_S3dt-1){}  
    62046086         else  
    6205            //    #pragma omp atomic 
     6087         #pragma omp atomic 
    62066088         YG1_ua( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[1]; 
    6207          //      #pragma omp atomic 
     6089         #pragma omp atomic 
    62086090         YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[2]; 
    62096091         if (1==0 || Yi+1>YA1_S3dt-1){}  
    62106092         else  
    6211            //    #pragma omp atomic 
     6093         #pragma omp atomic 
    62126094         YG1_va( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[3]; 
    6213          Ydbg_nanf ("B: ", YNBI_zwxyz_vor_ene, "zwxyz_vor_ene"); 
    6214  
    6215 } 
    6216  
    6217 } 
    6218 } 
    6219 //#pragma omp parallel for num_threads(16) 
     6095 
     6096} 
     6097 
     6098} 
     6099} 
     6100#pragma omp parallel for num_threads(16) 
    62206101for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    62216102  Yk=YY; 
    … …  
    62326113         Yva_dynadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] ); 
    62336114         Yvsmatt (YNBS_va_dynadv_cen2, YNBI_va_dynadv_cen2, YMAX_JAC_NBI, &YG1Y_va_dynadv_cen2, Yjac[0], Ytbeta); 
    6234          //      #pragma omp atomic 
     6115         #pragma omp atomic 
    62356116         YG2_zfuv_UW_dynadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    62366117         if (1==0 || Yk+1>YA3_S3dt-1){}  
    62376118         else  
    6238            //    #pragma omp atomic 
     6119         #pragma omp atomic 
    62396120         YG2_zfuv_UW_dynadv_cen2( Yi, Yj, Yk+1, YTemps) += Ytbeta[1]; 
    62406121         YG1_zva_dhmf_dynadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[2]; 
    6241          Ydbg_nanf ("B: ", YNBI_va_dynadv_cen2, "va_dynadv_cen2"); 
    62426122         Yting[0]=YS1_zfuv_UW_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    62436123         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    62496129         Yua_dynadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] ); 
    62506130         Yvsmatt (YNBS_ua_dynadv_cen2, YNBI_ua_dynadv_cen2, YMAX_JAC_NBI, &YG1Y_ua_dynadv_cen2, Yjac[0], Ytbeta); 
    6251          //      #pragma omp atomic 
     6131         #pragma omp atomic 
    62526132         YG1_zfuv_UW_dynadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    62536133         if (1==0 || Yk+1>YA3_S3dt-1){}  
    62546134         else  
    6255            //    #pragma omp atomic 
     6135         #pragma omp atomic 
    62566136         YG1_zfuv_UW_dynadv_cen2( Yi, Yj, Yk+1, YTemps) += Ytbeta[1]; 
    62576137         YG1_zua_dhmf_dynadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[2]; 
    6258          Ydbg_nanf ("B: ", YNBI_ua_dynadv_cen2, "ua_dynadv_cen2"); 
    6259  
    6260 } 
    6261 } 
    6262 //#pragma omp parallel for num_threads(16) 
     6138 
     6139} 
     6140} 
     6141#pragma omp parallel for num_threads(16) 
    62636142for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    62646143  Yk=YY; 
    … …  
    62886167         Yzfuv_UW_dynadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6] ); 
    62896168         Yvsmatt (YNBS_zfuv_UW_dynadv_cen2, YNBI_zfuv_UW_dynadv_cen2, YMAX_JAC_NBI, &YG1Y_zfuv_UW_dynadv_cen2, Yjac[0], Ytbeta); 
    6290          //      #pragma omp atomic 
     6169         #pragma omp atomic 
    62916170         YG1_wa( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    62926171         if (1==0 || Yi+1>YA1_S3dt-1){}  
    62936172         else  
    6294            //    #pragma omp atomic 
     6173         #pragma omp atomic 
    62956174         YG1_wa( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    62966175         if (1==0 || Yj+1>YA2_S3dt-1){}  
    62976176         else  
    6298            //    #pragma omp atomic 
     6177         #pragma omp atomic 
    62996178         YG1_wa( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[2]; 
    6300          //      #pragma omp atomic 
     6179         #pragma omp atomic 
    63016180         YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[3]; 
    63026181         if (1==0 || Yk-1<0){}  
    63036182         else  
    6304            //    #pragma omp atomic 
     6183         #pragma omp atomic 
    63056184         YG1_ua( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[4]; 
    6306          //      #pragma omp atomic 
     6185         #pragma omp atomic 
    63076186         YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[5]; 
    63086187         if (1==0 || Yk-1<0){}  
    63096188         else  
    6310            //    #pragma omp atomic 
     6189         #pragma omp atomic 
    63116190         YG1_va( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[6]; 
    6312          Ydbg_nanf ("B: ", YNBI_zfuv_UW_dynadv_cen2, "zfuv_UW_dynadv_cen2"); 
    6313  
    6314 } 
    6315 } 
    6316 //#pragma omp parallel for num_threads(16) 
     6191 
     6192} 
     6193} 
     6194#pragma omp parallel for num_threads(16) 
    63176195for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    63186196  Yk=YY; 
    … …  
    63426220         else  
    63436221         YG4_hmf_TF_dynadv_cen2( Yi, Yj-1, Yk, YTemps) += Ytbeta[3]; 
    6344          Ydbg_nanf ("B: ", YNBI_zva_dhmf_dynadv_cen2, "zva_dhmf_dynadv_cen2"); 
    63456222         Yting[0]=YS1_hmf_TF_dynadv_cen2( Yi, Yj, Yk, YTemps); 
    63466223         if (1==0 || Yi-1<0) 
    … …  
    63646241         else  
    63656242         YG2_hmf_TF_dynadv_cen2( Yi, Yj-1, Yk, YTemps) += Ytbeta[3]; 
    6366          Ydbg_nanf ("B: ", YNBI_zua_dhmf_dynadv_cen2, "zua_dhmf_dynadv_cen2"); 
    63676243 
    63686244} 
    … …  
    63916267         Yhmf_TF_dynadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5] ); 
    63926268         Yvsmatt (YNBS_hmf_TF_dynadv_cen2, YNBI_hmf_TF_dynadv_cen2, YMAX_JAC_NBI, &YG1Y_hmf_TF_dynadv_cen2, Yjac[0], Ytbeta); 
    6393          //      #pragma omp atomic 
     6269         #pragma omp atomic 
    63946270         YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    63956271         if (1==0 || Yi+1>YA1_S3dt-1){}  
    63966272         else  
    6397            //    #pragma omp atomic 
     6273         #pragma omp atomic 
    63986274         YG1_ua( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    63996275         if (1==0 || Yj+1>YA2_S3dt-1){}  
    64006276         else  
    6401            //    #pragma omp atomic 
     6277         #pragma omp atomic 
    64026278         YG1_ua( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[2]; 
    6403          //      #pragma omp atomic 
     6279         #pragma omp atomic 
    64046280         YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[3]; 
    64056281         if (1==0 || Yi+1>YA1_S3dt-1){}  
    64066282         else  
    6407            //    #pragma omp atomic 
     6283         #pragma omp atomic 
    64086284         YG1_va( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[4]; 
    64096285         if (1==0 || Yj+1>YA2_S3dt-1){}  
    64106286         else  
    6411            //    #pragma omp atomic 
     6287         #pragma omp atomic 
    64126288         YG1_va( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[5]; 
    6413          Ydbg_nanf ("B: ", YNBI_hmf_TF_dynadv_cen2, "hmf_TF_dynadv_cen2"); 
    6414  
    6415 } 
    6416  
    6417 } 
    6418 } 
    6419 //#pragma omp parallel for num_threads(16) 
     6289 
     6290} 
     6291 
     6292} 
     6293} 
     6294#pragma omp parallel for num_threads(16) 
    64206295for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    64216296  Yk=YY; 
    … …  
    64276302         Yrhop(Yi, Yj, Yk)->backward( Yting[0] ); 
    64286303         Yvsmatt (YNBS_rhop, YNBI_rhop, YMAX_JAC_NBI, &YG1Y_rhop, Yjac[0], Ytbeta); 
    6429          //      #pragma omp atomic 
     6304         #pragma omp atomic 
    64306305         YG1_rhd( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    6431          Ydbg_nanf ("B: ", YNBI_rhop, "rhop"); 
    64326306         Yting[0]=YS1_tb( Yi, Yj, Yk, YTemps); 
    64336307         Yting[1]=YS1_sb( Yi, Yj, Yk, YTemps); 
    … …  
    64376311         YG1_tb( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    64386312         YG1_sb( Yi, Yj, Yk, YTemps) += Ytbeta[1]; 
    6439          Ydbg_nanf ("B: ", YNBI_rhd, "rhd"); 
    6440  
    6441 } 
    6442 } 
    6443 //#pragma omp parallel for num_threads(16) 
     6313 
     6314} 
     6315} 
     6316#pragma omp parallel for num_threads(16) 
    64446317for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    64456318  Yk=YY; 
    … …  
    64546327         Ysb(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    64556328         Yvsmatt (YNBS_sb, YNBI_sb, YMAX_JAC_NBI, &YG1Y_sb, Yjac[0], Ytbeta); 
    6456          //      #pragma omp atomic 
     6329         #pragma omp atomic 
    64576330         YG1_sb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    6458          //      #pragma omp atomic 
     6331         #pragma omp atomic 
    64596332         YG1_sa( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    6460          //      #pragma omp atomic 
     6333         #pragma omp atomic 
    64616334         YG1_sa( Yi, Yj, Yk, YTemps) += Ytbeta[2]; 
    64626335         YG1_sa_c( Yi, Yj, Yk) += Ytbeta[3]; 
    6463          Ydbg_nanf ("B: ", YNBI_sb, "sb"); 
    64646336         Yting[0]=YS1_tb( Yi, Yj, Yk, YTemps-1); 
    64656337         Yting[1]=YS1_ta( Yi, Yj, Yk, YTemps-1); 
    … …  
    64696341         Ytb(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    64706342         Yvsmatt (YNBS_tb, YNBI_tb, YMAX_JAC_NBI, &YG1Y_tb, Yjac[0], Ytbeta); 
    6471          //      #pragma omp atomic 
     6343         #pragma omp atomic 
    64726344         YG1_tb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    6473          //      #pragma omp atomic 
     6345         #pragma omp atomic 
    64746346         YG1_ta( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    6475          //      #pragma omp atomic 
     6347         #pragma omp atomic 
    64766348         YG1_ta( Yi, Yj, Yk, YTemps) += Ytbeta[2]; 
    64776349         YG1_ta_c( Yi, Yj, Yk) += Ytbeta[3]; 
    6478          Ydbg_nanf ("B: ", YNBI_tb, "tb"); 
    64796350 
    64806351} 
    … …  
    65026373         YG1_zwt_trazdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[3]; 
    65036374         YG1_sa_c( Yi, Yj, Yk) += Ytbeta[4]; 
    6504          Ydbg_nanf ("B: ", YNBI_sa, "sa"); 
    65056375         Yting[0]=YS2_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps); 
    65066376         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    65216391         YG1_zwt_trazdf_imp( Yi, Yj, Yk, YTemps) += Ytbeta[3]; 
    65226392         YG1_ta_c( Yi, Yj, Yk) += Ytbeta[4]; 
    6523          Ydbg_nanf ("B: ", YNBI_ta, "ta"); 
    65246393 
    65256394} 
    … …  
    65496418         else  
    65506419         YG1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3]; 
    6551          Ydbg_nanf ("B: ", YNBI_sa_de_trazdf_imp, "sa_de_trazdf_imp"); 
    65526420         Yting[0]=YS1_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps); 
    65536421         if (1==0 || Yk-1<0) 
    … …  
    65716439         else  
    65726440         YG1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3]; 
    6573          Ydbg_nanf ("B: ", YNBI_ta_de_trazdf_imp, "ta_de_trazdf_imp"); 
    65746441         Yting[0]=YS1_zwisd_trazdf_imp( Yi, Yj, Yk, YTemps); 
    65756442         if (1==0 || Yk-1<0) 
    … …  
    65936460         else  
    65946461         YG1_zwt_trazdf_imp( Yi, Yj, Yk-1, YTemps) += Ytbeta[3]; 
    6595          Ydbg_nanf ("B: ", YNBI_zwt_trazdf_imp, "zwt_trazdf_imp"); 
    6596  
    6597 } 
    6598 //#pragma omp parallel for num_threads(16) 
     6462 
     6463} 
     6464#pragma omp parallel for num_threads(16) 
    65996465for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    66006466  Yk=YY; 
    … …  
    66086474         Yvsmatt (YNBS_zrhs_sa_trazdf_imp, YNBI_zrhs_sa_trazdf_imp, YMAX_JAC_NBI, &YG1Y_zrhs_sa_trazdf_imp, Yjac[0], Ytbeta); 
    66096475         YG1_sa_traldf_lap( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    6610          //      #pragma omp atomic 
     6476         #pragma omp atomic 
    66116477         YG1_sb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    6612          Ydbg_nanf ("B: ", YNBI_zrhs_sa_trazdf_imp, "zrhs_sa_trazdf_imp"); 
    66136478         Yting[0]=YS1_ta_traldf_lap( Yi, Yj, Yk, YTemps); 
    66146479         Yting[1]=YS1_tb( Yi, Yj, Yk, YTemps-1); 
    … …  
    66176482         Yvsmatt (YNBS_zrhs_trazdf_imp, YNBI_zrhs_trazdf_imp, YMAX_JAC_NBI, &YG1Y_zrhs_trazdf_imp, Yjac[0], Ytbeta); 
    66186483         YG1_ta_traldf_lap( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    6619          //      #pragma omp atomic 
     6484         #pragma omp atomic 
    66206485         YG1_tb( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    6621          Ydbg_nanf ("B: ", YNBI_zrhs_trazdf_imp, "zrhs_trazdf_imp"); 
    66226486         Yting[0]=YS1_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps); 
    66236487         if (1==0 || Yk+1>YA3_S3dt-1) 
    … …  
    66286492         Yzwisd_trazdf_imp(Yi, Yj, Yk)->backward( Yting[0], Yting[1] ); 
    66296493         Yvsmatt (YNBS_zwisd_trazdf_imp, YNBI_zwisd_trazdf_imp, YMAX_JAC_NBI, &YG1Y_zwisd_trazdf_imp, Yjac[0], Ytbeta); 
    6630          //      #pragma omp atomic 
     6494         #pragma omp atomic 
    66316495         YG1_zdf_avt_avmu_avmv( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    66326496         if (1==0 || Yk+1>YA3_S3dt-1){}  
    66336497         else  
    6634            //    #pragma omp atomic 
     6498         #pragma omp atomic 
    66356499         YG1_zdf_avt_avmu_avmv( Yi, Yj, Yk+1, YTemps) += Ytbeta[1]; 
    6636          Ydbg_nanf ("B: ", YNBI_zwisd_trazdf_imp, "zwisd_trazdf_imp"); 
    6637  
    6638 } 
    6639 } 
    6640 //#pragma omp parallel for num_threads(16) 
     6500 
     6501} 
     6502} 
     6503#pragma omp parallel for num_threads(16) 
    66416504for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    66426505  Yk=YY; 
    … …  
    66676530         else  
    66686531         YG2_zsuv_traldf_lap( Yi, Yj-1, Yk, YTemps) += Ytbeta[4]; 
    6669          Ydbg_nanf ("B: ", YNBI_sa_traldf_lap, "sa_traldf_lap"); 
    66706532         Yting[0]=YS1_ta_ver_traadv_cen2( Yi, Yj, Yk, YTemps); 
    66716533         Yting[1]=YS1_ztuv_traldf_lap( Yi, Yj, Yk, YTemps); 
    … …  
    66916553         else  
    66926554         YG2_ztuv_traldf_lap( Yi, Yj-1, Yk, YTemps) += Ytbeta[4]; 
    6693          Ydbg_nanf ("B: ", YNBI_ta_traldf_lap, "ta_traldf_lap"); 
    66946555         Yting[0]=YS1_sb( Yi, Yj, Yk, YTemps-1); 
    66956556         if (1==0 || Yi+1>YA1_S3dt-1) 
    … …  
    67046565         Yzsuv_traldf_lap(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] ); 
    67056566         Yvsmatt (YNBS_zsuv_traldf_lap, YNBI_zsuv_traldf_lap, YMAX_JAC_NBI, &YG1Y_zsuv_traldf_lap, Yjac[0], Ytbeta); 
    6706          //      #pragma omp atomic 
     6567         #pragma omp atomic 
    67076568         YG1_sb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    67086569         if (1==0 || Yi+1>YA1_S3dt-1){}  
    67096570         else  
    6710            //    #pragma omp atomic 
     6571         #pragma omp atomic 
    67116572         YG1_sb( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    67126573         if (1==0 || Yj+1>YA2_S3dt-1){}  
    67136574         else  
    6714            //    #pragma omp atomic 
     6575         #pragma omp atomic 
    67156576         YG1_sb( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[2]; 
    6716          Ydbg_nanf ("B: ", YNBI_zsuv_traldf_lap, "zsuv_traldf_lap"); 
    67176577         Yting[0]=YS1_tb( Yi, Yj, Yk, YTemps-1); 
    67186578         if (1==0 || Yi+1>YA1_S3dt-1) 
    … …  
    67276587         Yztuv_traldf_lap(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] ); 
    67286588         Yvsmatt (YNBS_ztuv_traldf_lap, YNBI_ztuv_traldf_lap, YMAX_JAC_NBI, &YG1Y_ztuv_traldf_lap, Yjac[0], Ytbeta); 
    6729          //      #pragma omp atomic 
     6589         #pragma omp atomic 
    67306590         YG1_tb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    67316591         if (1==0 || Yi+1>YA1_S3dt-1){}  
    67326592         else  
    6733            //    #pragma omp atomic 
     6593         #pragma omp atomic 
    67346594         YG1_tb( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    67356595         if (1==0 || Yj+1>YA2_S3dt-1){}  
    67366596         else  
    6737            //    #pragma omp atomic 
     6597         #pragma omp atomic 
    67386598         YG1_tb( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[2]; 
    6739          Ydbg_nanf ("B: ", YNBI_ztuv_traldf_lap, "ztuv_traldf_lap"); 
    6740  
    6741 } 
    6742 } 
    6743 //#pragma omp parallel for num_threads(16) 
     6599 
     6600} 
     6601} 
     6602#pragma omp parallel for num_threads(16) 
    67446603for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    67456604  Yk=YY; 
    … …  
    67576616         Yvsmatt (YNBS_sa_ver_traadv_cen2, YNBI_sa_ver_traadv_cen2, YMAX_JAC_NBI, &YG1Y_sa_ver_traadv_cen2, Yjac[0], Ytbeta); 
    67586617         YG1_sa_hor_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    6759          //      #pragma omp atomic 
     6618         #pragma omp atomic 
    67606619         YG1_zcens_ver_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[1]; 
    67616620         if (1==0 || Yk+1>YA3_S3dt-1){}  
    67626621         else  
    6763            //    #pragma omp atomic 
     6622         #pragma omp atomic 
    67646623         YG1_zcens_ver_traadv_cen2( Yi, Yj, Yk+1, YTemps) += Ytbeta[2]; 
    6765          Ydbg_nanf ("B: ", YNBI_sa_ver_traadv_cen2, "sa_ver_traadv_cen2"); 
    67666624         Yting[0]=YS1_ta_hor_traadv_cen2( Yi, Yj, Yk, YTemps); 
    67676625         Yting[1]=YS1_zcent_ver_traadv_cen2( Yi, Yj, Yk, YTemps); 
    … …  
    67746632         Yvsmatt (YNBS_ta_ver_traadv_cen2, YNBI_ta_ver_traadv_cen2, YMAX_JAC_NBI, &YG1Y_ta_ver_traadv_cen2, Yjac[0], Ytbeta); 
    67756633         YG1_ta_hor_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[0]; 
    6776          //      #pragma omp atomic 
     6634         #pragma omp atomic 
    67776635         YG1_zcent_ver_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[1]; 
    67786636         if (1==0 || Yk+1>YA3_S3dt-1){}  
    67796637         else  
    6780            //    #pragma omp atomic 
     6638         #pragma omp atomic 
    67816639         YG1_zcent_ver_traadv_cen2( Yi, Yj, Yk+1, YTemps) += Ytbeta[2]; 
    6782          Ydbg_nanf ("B: ", YNBI_ta_ver_traadv_cen2, "ta_ver_traadv_cen2"); 
    6783  
    6784 } 
    6785 } 
    6786 //#pragma omp parallel for num_threads(16) 
     6640 
     6641} 
     6642} 
     6643#pragma omp parallel for num_threads(16) 
    67876644for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    67886645  Yk=YY; 
    … …  
    67996656         Yzcens_ver_traadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] ); 
    68006657         Yvsmatt (YNBS_zcens_ver_traadv_cen2, YNBI_zcens_ver_traadv_cen2, YMAX_JAC_NBI, &YG1Y_zcens_ver_traadv_cen2, Yjac[0], Ytbeta); 
    6801          //      #pragma omp atomic 
     6658         #pragma omp atomic 
    68026659         YG1_wa( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    6803          //      #pragma omp atomic 
     6660         #pragma omp atomic 
    68046661         YG1_sa( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    68056662         if (1==0 || Yk-1<0){}  
    68066663         else  
    6807            //    #pragma omp atomic 
     6664         #pragma omp atomic 
    68086665         YG1_sa( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[2]; 
    6809          Ydbg_nanf ("B: ", YNBI_zcens_ver_traadv_cen2, "zcens_ver_traadv_cen2"); 
    68106666         Yting[0]=YS1_wa( Yi, Yj, Yk, YTemps-1); 
    68116667         Yting[1]=YS1_ta( Yi, Yj, Yk, YTemps-1); 
    … …  
    68176673         Yzcent_ver_traadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2] ); 
    68186674         Yvsmatt (YNBS_zcent_ver_traadv_cen2, YNBI_zcent_ver_traadv_cen2, YMAX_JAC_NBI, &YG1Y_zcent_ver_traadv_cen2, Yjac[0], Ytbeta); 
    6819          //      #pragma omp atomic 
     6675         #pragma omp atomic 
    68206676         YG1_wa( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    6821          //      #pragma omp atomic 
     6677         #pragma omp atomic 
    68226678         YG1_ta( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    68236679         if (1==0 || Yk-1<0){}  
    68246680         else  
    6825            //    #pragma omp atomic 
     6681         #pragma omp atomic 
    68266682         YG1_ta( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[2]; 
    6827          Ydbg_nanf ("B: ", YNBI_zcent_ver_traadv_cen2, "zcent_ver_traadv_cen2"); 
    6828  
    6829 } 
    6830 } 
    6831 //#pragma omp parallel for num_threads(16) 
     6683 
     6684} 
     6685} 
     6686#pragma omp parallel for num_threads(16) 
    68326687for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    68336688  Yk=YY; 
    … …  
    68596714         YG2_zcens_hor_traadv_cen2( Yi, Yj-1, Yk, YTemps) += Ytbeta[3]; 
    68606715         YG2_zcens_hor_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[4]; 
    6861          Ydbg_nanf ("B: ", YNBI_sa_hor_traadv_cen2, "sa_hor_traadv_cen2"); 
    68626716         Yting[0]=YS1_ta_traqsr( Yi, Yj, Yk, YTemps); 
    68636717         if (1==0 || Yi-1<0) 
    … …  
    68836737         YG2_zcent_hor_traadv_cen2( Yi, Yj-1, Yk, YTemps) += Ytbeta[3]; 
    68846738         YG2_zcent_hor_traadv_cen2( Yi, Yj, Yk, YTemps) += Ytbeta[4]; 
    6885          Ydbg_nanf ("B: ", YNBI_ta_hor_traadv_cen2, "ta_hor_traadv_cen2"); 
    68866739 
    68876740} 
    … …  
    69036756         Yzcens_hor_traadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    69046757         Yvsmatt (YNBS_zcens_hor_traadv_cen2, YNBI_zcens_hor_traadv_cen2, YMAX_JAC_NBI, &YG1Y_zcens_hor_traadv_cen2, Yjac[0], Ytbeta); 
    6905          //      #pragma omp atomic 
     6758         #pragma omp atomic 
    69066759         YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    6907          //      #pragma omp atomic 
     6760         #pragma omp atomic 
    69086761         YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    6909          //      #pragma omp atomic 
     6762         #pragma omp atomic 
    69106763         YG1_sa( Yi, Yj, Yk, YTemps-1) += Ytbeta[2]; 
    69116764         if (1==0 || Yj+1>YA2_S3dt-1){}  
    69126765         else  
    6913            //    #pragma omp atomic 
     6766         #pragma omp atomic 
    69146767         YG1_sa( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[3]; 
    69156768         if (1==0 || Yi+1>YA1_S3dt-1){}  
    69166769         else  
    6917            //    #pragma omp atomic 
     6770         #pragma omp atomic 
    69186771         YG1_sa( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[4]; 
    6919          Ydbg_nanf ("B: ", YNBI_zcens_hor_traadv_cen2, "zcens_hor_traadv_cen2"); 
    69206772         Yting[0]=YS1_ua( Yi, Yj, Yk, YTemps-1); 
    69216773         Yting[1]=YS1_va( Yi, Yj, Yk, YTemps-1); 
    … …  
    69326784         Yzcent_hor_traadv_cen2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4] ); 
    69336785         Yvsmatt (YNBS_zcent_hor_traadv_cen2, YNBI_zcent_hor_traadv_cen2, YMAX_JAC_NBI, &YG1Y_zcent_hor_traadv_cen2, Yjac[0], Ytbeta); 
    6934          //      #pragma omp atomic 
     6786         #pragma omp atomic 
    69356787         YG1_ua( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    6936          //      #pragma omp atomic 
     6788         #pragma omp atomic 
    69376789         YG1_va( Yi, Yj, Yk, YTemps-1) += Ytbeta[1]; 
    6938          //      #pragma omp atomic 
     6790         #pragma omp atomic 
    69396791         YG1_ta( Yi, Yj, Yk, YTemps-1) += Ytbeta[2]; 
    69406792         if (1==0 || Yj+1>YA2_S3dt-1){}  
    69416793         else  
    6942            //    #pragma omp atomic 
     6794         #pragma omp atomic 
    69436795         YG1_ta( Yi, Yj+1, Yk, YTemps-1) += Ytbeta[3]; 
    69446796         if (1==0 || Yi+1>YA1_S3dt-1){}  
    69456797         else  
    6946            //    #pragma omp atomic 
     6798         #pragma omp atomic 
    69476799         YG1_ta( Yi+1, Yj, Yk, YTemps-1) += Ytbeta[4]; 
    6948          Ydbg_nanf ("B: ", YNBI_zcent_hor_traadv_cen2, "zcent_hor_traadv_cen2"); 
    69496800 
    69506801} 
    … …  
    69656816         YG1_sbcgyre_qsr_qns( Yi, Yj, YTemps) += Ytbeta[0]; 
    69666817         YG1_ta_sa_trasbc( Yi, Yj, YTemps) += Ytbeta[1]; 
    6967          Ydbg_nanf ("B: ", YNBI_ta_traqsr, "ta_traqsr"); 
    69686818 
    69696819} 
    … …  
    69846834         #pragma omp atomic 
    69856835         YG1_sa( Yi, Yj, 0, YTemps-1) += Ytbeta[2]; 
    6986          Ydbg_nanf ("B: ", YNBI_ta_sa_trasbc, "ta_sa_trasbc"); 
    69876836 
    69886837} 
    … …  
    70586907         YG1_rhop( Yi, Yj, 29, YTemps-1) += Ytbeta[29]; 
    70596908         YG1_rhop( Yi, Yj, 30, YTemps-1) += Ytbeta[30]; 
    7060          Ydbg_nanf ("B: ", YNBI_zdfmxl_nmln_hmlp_hmlpt, "zdfmxl_nmln_hmlp_hmlpt"); 
    70616909         Yting[0]=YS1_zdf_avt_avmu_avmv( Yi, Yj, 0, YTemps); 
    70626910         Yting[1]=YS1_zdf_avt_avmu_avmv( Yi, Yj, 1, YTemps); 
    … …  
    71246972         YG1_zdf_avt_avmu_avmv( Yi, Yj, 29, YTemps) += Ytbeta[29]; 
    71256973         YG1_zdf_avt_avmu_avmv( Yi, Yj, 30, YTemps) += Ytbeta[30]; 
    7126          Ydbg_nanf ("B: ", YNBI_zdfmxl_hmld, "zdfmxl_hmld"); 
    7127  
    7128 } 
    7129 } 
    7130 //#pragma omp parallel for num_threads(16) 
     6974 
     6975} 
     6976} 
     6977#pragma omp parallel for num_threads(16) 
    71316978for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    71326979  Yk=YY; 
    … …  
    72077054         if (1==0 || Yi-1<0){}  
    72087055         else  
    7209            //    #pragma omp atomic 
     7056         #pragma omp atomic 
    72107057         YG1_ua( Yi-1, Yj, 29, YTemps-1) += Ytbeta[3]; 
    72117058         if (1==0 || Yi-1<0 || Yj+1>YA2_S3dt-1){}  
    72127059         else  
    7213            //    #pragma omp atomic 
     7060         #pragma omp atomic 
    72147061         YG1_ua( Yi-1, Yj+1, 29, YTemps-1) += Ytbeta[4]; 
    7215          //      #pragma omp atomic 
     7062         #pragma omp atomic 
    72167063         YG1_ua( Yi, Yj, 29, YTemps-1) += Ytbeta[5]; 
    72177064         if (1==0 || Yj+1>YA2_S3dt-1){}  
    72187065         else  
    7219            //    #pragma omp atomic 
     7066         #pragma omp atomic 
    72207067         YG1_ua( Yi, Yj+1, 29, YTemps-1) += Ytbeta[6]; 
    72217068         if (1==0 || Yi-1<0){}  
    72227069         else  
    7223            //    #pragma omp atomic 
     7070         #pragma omp atomic 
    72247071         YG1_ua( Yi-1, Yj, 1, YTemps-1) += Ytbeta[7]; 
    72257072         if (1==0 || Yi-1<0 || Yj+1>YA2_S3dt-1){}  
    72267073         else  
    7227            //    #pragma omp atomic 
     7074         #pragma omp atomic 
    72287075         YG1_ua( Yi-1, Yj+1, 1, YTemps-1) += Ytbeta[8]; 
    7229          //      #pragma omp atomic 
     7076         #pragma omp atomic 
    72307077         YG1_ua( Yi, Yj, 1, YTemps-1) += Ytbeta[9]; 
    72317078         if (1==0 || Yj+1>YA2_S3dt-1){}  
    72327079         else  
    7233            //    #pragma omp atomic 
     7080         #pragma omp atomic 
    72347081         YG1_ua( Yi, Yj+1, 1, YTemps-1) += Ytbeta[10]; 
    72357082         if (1==0 || Yj-1<0){}  
    72367083         else  
    7237            //    #pragma omp atomic 
     7084         #pragma omp atomic 
    72387085         YG1_va( Yi, Yj-1, 29, YTemps-1) += Ytbeta[11]; 
    7239          //      #pragma omp atomic 
     7086         #pragma omp atomic 
    72407087         YG1_va( Yi, Yj, 29, YTemps-1) += Ytbeta[12]; 
    72417088         if (1==0 || Yi+1>YA1_S3dt-1 || Yj-1<0){}  
    72427089         else  
    7243            //    #pragma omp atomic 
     7090         #pragma omp atomic 
    72447091         YG1_va( Yi+1, Yj-1, 29, YTemps-1) += Ytbeta[13]; 
    72457092         if (1==0 || Yi+1>YA1_S3dt-1){}  
    72467093         else  
    7247            //    #pragma omp atomic 
     7094         #pragma omp atomic 
    72487095         YG1_va( Yi+1, Yj, 29, YTemps-1) += Ytbeta[14]; 
    72497096         if (1==0 || Yj-1<0){}  
    72507097         else  
    7251            //    #pragma omp atomic 
     7098         #pragma omp atomic 
    72527099         YG1_va( Yi, Yj-1, 1, YTemps-1) += Ytbeta[15]; 
    7253          //      #pragma omp atomic 
     7100         #pragma omp atomic 
    72547101         YG1_va( Yi, Yj, 1, YTemps-1) += Ytbeta[16]; 
    72557102         if (1==0 || Yi+1>YA1_S3dt-1 || Yj-1<0){}  
    72567103         else  
    7257            //    #pragma omp atomic 
     7104         #pragma omp atomic 
    72587105         YG1_va( Yi+1, Yj-1, 1, YTemps-1) += Ytbeta[17]; 
    72597106         if (1==0 || Yi+1>YA1_S3dt-1){}  
    72607107         else  
    7261            // #pragma omp atomic 
     7108         #pragma omp atomic 
    72627109         YG1_va( Yi+1, Yj, 1, YTemps-1) += Ytbeta[18]; 
    7263          Ydbg_nanf ("B: ", YNBI_zdf_avt_avmu_avmv, "zdf_avt_avmu_avmv"); 
    7264  
    7265 } 
    7266 } 
    7267 //#pragma omp parallel for num_threads(16) 
    7268 for(YY=YA3_S3dt-1; YY>=0; --YY){  
     7110 
     7111} 
     7112} 
     7113#pragma omp parallel for num_threads(16) 
     7114for(YY=YA3_S3dt-1; YY>=0; --YY){ 
    72697115  Yk=YY; 
    72707116for(Yj=YA2_S3dt-1; Yj>=0; --Yj) 
    … …  
    72847130         Ybn2_rn2(Yi, Yj, Yk)->backward( Yting[0], Yting[1], Yting[2], Yting[3] ); 
    72857131         Yvsmatt (YNBS_bn2_rn2, YNBI_bn2_rn2, YMAX_JAC_NBI, &YG1Y_bn2_rn2, Yjac[0], Ytbeta); 
    7286          //      #pragma omp atomic 
     7132         #pragma omp atomic 
    72877133         YG1_tb( Yi, Yj, Yk, YTemps-1) += Ytbeta[0]; 
    72887134         if (1==0 || Yk-1<0){}  
    72897135         else  
    7290            //    #pragma omp atomic 
     7136         #pragma omp atomic 
    72917137         YG1_tb( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[1]; 
    7292          //      #pragma omp atomic 
     7138         #pragma omp atomic 
    72937139         YG1_sb( Yi, Yj, Yk, YTemps-1) += Ytbeta[2]; 
    72947140         if (1==0 || Yk-1<0){}  
    72957141         else  
    7296            //    #pragma omp atomic 
     7142         #pragma omp atomic 
    72977143         YG1_sb( Yi, Yj, Yk-1, YTemps-1) += Ytbeta[3]; 
    7298          Ydbg_nanf ("B: ", YNBI_bn2_rn2, "bn2_rn2"); 
    72997144 
    73007145} 
    … …  
    73247169         YG1_sbcgyre_utau_vtau( Yi, Yj, YTemps) += Ytbeta[2]; 
    73257170         YG2_sbcgyre_utau_vtau( Yi, Yj, YTemps) += Ytbeta[3]; 
    7326          Ydbg_nanf ("B: ", YNBI_sbcgyre_wndm, "sbcgyre_wndm"); 
    73277171         Ysbcgyre_utau_vtau(Yi, Yj)->backward(); 
    73287172         Yting[0]=YS1_sbcgyre_emp( Yi, Yj, YTemps); 
    … …  
    73317175         Yvsmatt (YNBS_sbcgyre_emps, YNBI_sbcgyre_emps, YMAX_JAC_NBI, &YG1Y_sbcgyre_emps, Yjac[0], Ytbeta); 
    73327176         YG1_sbcgyre_emp( Yi, Yj, YTemps) += Ytbeta[0]; 
    7333          Ydbg_nanf ("B: ", YNBI_sbcgyre_emps, "sbcgyre_emps"); 
    73347177         Yting[0]=YS1_sbcgyre_emp1( Yi, Yj, YTemps); 
    73357178         Yting[1]=YS1_sbcgyre_zsumemp_zsurf( 31, 21, YTemps); 
    … …  
    73417184         YG1_sbcgyre_zsumemp_zsurf( 31, 21, YTemps) += Ytbeta[1]; 
    73427185         YG2_sbcgyre_zsumemp_zsurf( 31, 21, YTemps) += Ytbeta[2]; 
    7343          Ydbg_nanf ("B: ", YNBI_sbcgyre_emp, "sbcgyre_emp"); 
    73447186 
    73457187} 
    … …  
    73807222         YG2_sbcgyre_zsumemp_zsurf( 31, Yj-1, YTemps) += Ytbeta[3]; 
    73817223         YG1_sbcgyre_emp1( Yi, Yj, YTemps) += Ytbeta[4]; 
    7382          Ydbg_nanf ("B: ", YNBI_sbcgyre_zsumemp_zsurf, "sbcgyre_zsumemp_zsurf"); 
    73837224         Ysbcgyre_emp1(Yi, Yj)->backward(); 
    73847225         Yting[0]=YS1_tb( Yi, Yj, 0, YTemps-1); 
    … …  
    73877228         Yvsmatt (YNBS_sbcgyre_qsr_qns, YNBI_sbcgyre_qsr_qns, YMAX_JAC_NBI, &YG1Y_sbcgyre_qsr_qns, Yjac[0], Ytbeta); 
    73887229         YG1_tb( Yi, Yj, 0, YTemps-1) += Ytbeta[0]; 
    7389          Ydbg_nanf ("B: ", YNBI_sbcgyre_qsr_qns, "sbcgyre_qsr_qns"); 
    73907230 
    73917231} 

  • trunk/src/cost.dat

    r37 r41  
    22 0.0000000000000000e+00 
    33 0.0000000000000000e+00 
     4 2.9733593674719423e+03 
     5 5.4561410776679122e+02 
     6 2.4475352252509134e+02 
     7 2.6358555318068841e+02 
     8 2.1527711430855038e+02 
     9 2.0204126463990215e+02 
     10 2.1416662011590833e+02 
     11 2.0207208737645496e+02 
     12 2.0211049147089304e+02 
     13 2.0235703823123197e+02 
     14 2.0019969413511828e+02 
     15 2.2286623273509434e+02 
     16 1.9878349444274735e+02 
     17 2.0070292696001030e+02 
     18 2.0014258993979595e+02 
     19 1.9972801663599338e+02 
     20 1.9293537153893590e+02 
     21 2.0814944891652880e+02 
     22 1.8213155756443817e+02 
     23 1.7950672225933866e+02 
     24 1.8080233872884580e+02 
     25 1.7876027971923725e+02 
     26 2.0259497236690703e+02 
     27 1.7881103040469617e+02 
     28 1.7886996041817454e+02 
     29 1.7907554124113310e+02 
     30 1.8045371934564810e+02 
     31 1.7569372841627006e+02 
     32 1.9113648527038990e+02 
     33 1.7516955251891332e+02 
     34 1.9550146370923514e+02 
     35 1.7455342251319595e+02 
     36 2.1077987977826763e+02 
     37 1.7473743190883036e+02 
     38 1.7476659805882767e+02 
     39 1.7491607437730161e+02 
     40 1.7534570209262407e+02 
     41 1.7146968928574532e+02 
     42 1.7918664114219453e+02 
     43 1.7311809756377912e+02 
     44 1.7320818032615955e+02 
     45 1.7322190117805096e+02 
     46 1.7328346667692077e+02 
     47 1.7390681350015083e+02 
     48 1.7091883954920155e+02 
     49 1.7074257389516592e+02 
     50 1.7423942813610816e+02 
     51 1.7069031205380722e+02 
     52 1.7425456057653693e+02 
     53 1.7425631257707835e+02 
     54 1.7425753342598696e+02 

  • trunk/src/solsor_flag.h

    r39 r41  
    1212        3       from    solsor_flag             1  i-1  j    k  t 
    1313        4       from    solsor_flag             1  i    j-1  k  t 
    14 } */ 
     14}  
     15output 1 : resmax 
     16oupput 2 : flag (0 -> continue; 1 -> stop) 
     17 
     18*/ 
    1519 
    1620//jn (current iteration) = Yk+1 
    … …  
    1822if (Yk<nmin) 
    1923{ 
    20 YS1=0; 
     24  YS1=0; 
     25 YS2=0; 
    2126} 
    22  
    23 if (Yi==0) { 
    24  
     27 else { 
     28   if (Yi==0 & Yj==0) YS1=0 
     29      
     30 } //Yk>=nmin 
Note: See TracChangeset for help on using the changeset viewer.